Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,668 C→T 6.0% intergenic (+16/‑417) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6680CT 6.0% 644.4 / 23.6 233intergenic (+16/‑417)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (122/97);  new base T (4/10);  total (126/107)
Fisher's exact test for biased strand distribution p-value = 5.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781519‑781817
                                                                                                                                                      |                                                                                                                                                     
tCGAGTCATCCTGCACGACCCACCATCCTGAATGATAAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGccc                                                                                                                                                       >  2:275921/1‑150 (MQ=255)
  gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTa                                                                                                                                                     >  1:5470101/1‑150 (MQ=255)
     taaTCCTGCACGACCCAACATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAg                                                                                                                                                  <  2:5072270/148‑1 (MQ=255)
     tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAg                                                                                                                                                  >  1:1319748/1‑150 (MQ=255)
        tCCTGCACGACCCACCATCCTGAATGAGTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGgcg                                                                                                                                               >  2:3016209/1‑150 (MQ=255)
        tCCTGCAAGACCCACCATCCTGATTTATTAAGTCATCATAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGgcg                                                                                                                                               <  2:3037766/150‑1 (MQ=17)
           tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTCAAGGCGCtt                                                                                                                                            >  2:3389496/1‑150 (MQ=255)
           tGCAAGACCCACCATCCTTATTGATTTAGGCAGCATAAGCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtt                                                                                                                                            <  2:4496273/150‑1 (MQ=17)
            gCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttt                                                                                                                                           <  1:1197788/150‑1 (MQ=255)
            gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttt                                                                                                                                           >  1:4495136/1‑150 (MQ=255)
             cccGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttt                                                                                                                                          <  1:3117922/148‑1 (MQ=255)
             cACGACCGACCATCCTTTTTTATTTATGCAGCCTAATCCCGGTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtttt                                                                                                                                          <  2:5344934/150‑1 (MQ=14)
             cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttt                                                                                                                                          <  2:1936042/150‑1 (MQ=255)
              aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttttt                                                                                                                                         >  1:4341404/1‑150 (MQ=255)
               cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttttt                                                                                                                                        <  2:1786100/150‑1 (MQ=255)
                  cccACCATCCTGAATGATTAAGGCAGCATACCCCCGCACGGGGTCGTTATCTCAGTGGTTAGAGCAGTTCACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAA‑GCGCCTTTAGGGCGCTTTTTTTCt                                                                                                                                     <  1:1147844/150‑1 (MQ=17)
                  cacaccaTCCTGTTTGATTAATGCATCTTAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCt                                                                                                                                     <  2:2062591/148‑1 (MQ=17)
                   ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAA‑GCGCCTTTAGGGCGCTTTTTTacat                                                                                                                                    >  2:5242796‑M2/1‑146 (MQ=255)
                    caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  1:4633906/1‑150 (MQ=255)
                    caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  2:4069015/1‑150 (MQ=255)
                      ccaTCCTGAATGATTCAGGCAGCATACTCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                 <  2:3285419/150‑1 (MQ=255)
                        cgcccgAAAGATTAAGGCAGCATACTCCCGCAAGGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:418786/145‑1 (MQ=17)
                        aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:1266848/150‑1 (MQ=255)
                          ccggcaagATTTAGGCCGCCTACGCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                             <  2:3804495/143‑1 (MQ=11)
                           cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATAAGCGAt                                                                                                                            >  2:4834408/1‑150 (MQ=255)
                            tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                           >  2:2633753/1‑150 (MQ=255)
                             gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATAc                                                                                                                          >  2:4570804/1‑150 (MQ=255)
                                tGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCa                                                                                                                       <  2:2265849/150‑1 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:4426116/150‑1 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:782611/150‑1 (MQ=255)
                                        gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                               >  2:1520730/1‑150 (MQ=255)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                             >  2:1341104/1‑150 (MQ=255)
                                            caTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACCAAcc                                                                                                           >  2:4113771/1‑150 (MQ=255)
                                             aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtttttt                                                                                                                                        <  2:3465197/120‑1 (MQ=255)
                                              tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                         <  1:1198174/150‑1 (MQ=255)
                                                 tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      >  2:3960605/1‑150 (MQ=39)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  2:4018495/150‑1 (MQ=33)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:662897/150‑1 (MQ=33)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:4678950/150‑1 (MQ=33)
                                                          gggTCGTTTGTCAAGTTGGTAGAGCTGTTGAATTTTAAACAATTGGCCGACGGTTCGAAGCCTGCACGACCCACCAATTTAAAAAA‑GCGCCCAAACGCCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  2:4632020/150‑1 (MQ=11)
                                                          gggTCGGTAGCTGAGTTTGTTGTTCTGTTGTCTTTTAATCTATTGGTCGCATGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  1:1413928/150‑1 (MQ=11)
                                                            gTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                           <  1:5209880/150‑1 (MQ=17)
                                                               gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                        >  1:2000925/1‑150 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            <  1:5292019/129‑1 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            >  2:5292019/1‑129 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                       >  1:4995909/1‑150 (MQ=32)
                                                                 tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                                      <  1:4289101/150‑1 (MQ=32)
                                                                   gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                    <  2:1319748/150‑1 (MQ=32)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   <  1:818955/150‑1 (MQ=32)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   >  1:5063542/1‑150 (MQ=32)
                                                                     tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                  <  1:2633628/150‑1 (MQ=32)
                                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCAAGGATTCGAACCTGCAGCAGGTGTGTGTTGAGCGCAGCGa                                                                               >  2:599863/1‑150 (MQ=255)
                                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                               >  2:3279174/1‑150 (MQ=17)
                                                                          tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                             <  1:4743045/150‑1 (MQ=32)
                                                                            gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                           >  1:3438665/1‑150 (MQ=17)
                                                                             tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                          <  2:4542379/150‑1 (MQ=32)
                                                                                 gaaGTTTAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCTAAACAACGGa                                                                      <  2:1869540/148‑1 (MQ=255)
                                                                                   aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    >  1:5101210/1‑150 (MQ=32)
                                                                                    gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                   >  1:3053346/1‑150 (MQ=32)
                                                                                     ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCAGAAGCAGTTAGGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                  >  1:2219301/1‑150 (MQ=11)
                                                                                      tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCCAACCTGCAGCAGGTTTGAGAGGAGCGCAGAGAAACAAAGGAGCCGc                                                                 >  1:1696659/1‑150 (MQ=11)
                                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                             >  2:1988189/1‑150 (MQ=32)
                                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACGCAAAGATACGATCCTGCACCAGTTTTCAGATGATCGCAGCGAAGCGCACGAGCCGCTcgc                                                             >  1:3612563/1‑150 (MQ=11)
                                                                                               aTCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGAAGTTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:3233045/150‑1 (MQ=11)
                                                                                               aTCAATTGGCCGCAGGTTCGAACCCTGCCCGACCCACCACTGTAAAAAC‑GCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:2360506/150‑1 (MQ=11)
                                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGATCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  2:2296445/1‑150 (MQ=17)
                                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  1:404248/1‑150 (MQ=17)
                                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                     >  2:1472973/1‑150 (MQ=32)
                                                                                                    ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   >  2:381594/1‑150 (MQ=17)
                                                                                                    ttGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   <  2:4120873/150‑1 (MQ=17)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTTAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               <  2:1110520/150‑1 (MQ=12)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:419337/1‑150 (MQ=32)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGCAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:4223706/1‑150 (MQ=17)
                                                                                                         cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATGTGCGATACTCAAATATTCGTACCTGCACCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg                                              >  1:4378854/1‑150 (MQ=11)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                        >  1:3682750/1‑150 (MQ=32)
                                                                                                                tCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTAGCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                       >  2:2846183/1‑150 (MQ=17)
                                                                                                                  gAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                     >  1:3017862/1‑150 (MQ=32)
                                                                                                                    aTCCTGCACGACCCACCAATTTAAAAAT‑GCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                   <  2:754114/150‑1 (MQ=255)
                                                                                                                       cTGCACGACCAACCAATGTAAAAAA‑GCGCCCTAATGGCGCTTTTTTGATATATGCGATACTAAAAGATTCTAACCTGCAGCTGGTTTTTGTTGAGCGCAGCTAAACAACGGATCCGATCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:4137630/150‑1 (MQ=11)
                                                                                                                        tGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                               >  2:5110497/1‑150 (MQ=32)
                                                                                                                           aCGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgcaggattcgaacctgcgaccaagtgagtacaagtcaactgctctaccaactgggctaacgacccgacgtggtgggtgatgacgggatcgaac                            <  2:4857347‑M2/150‑109 (MQ=255)
                                                                                                                           aCGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  2:4287498/150‑1 (MQ=32)
                                                                                                                           aCGACCCA‑CAATTTAAAAAACGCCCCCTACGGGGGTTTTTTTGCTTTATGCTATACTCAAAGATTCCACCCTGCTGCTGTTTTGTGTTGCGGGCTCCCAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  1:4954983/150‑1 (MQ=2)
                                                                                                                            cgccccACCACTGTAAAAAA‑CCGCCCTAAAGGCGCTTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCTGTTTTTAGTTGTGCGCAGCGAAACAACGGAGCCGCTCGCTGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                           <  2:1022029/147‑1 (MQ=11)
                                                                                                                             gACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTCTGAAAGATTAGAACCTGCAGCAGGTTTGAGTTGAGCGGAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCGGAAGGGCGAGCGAAGCGAGTCATc                          >  1:1843078/1‑150 (MQ=11)
                                                                                                                              aCCCACCACTGAAAACAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGACACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATcc                         <  2:4445137/150‑1 (MQ=11)
                                                                                                                                  accaATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAATCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:2005787/1‑150 (MQ=17)
                                                                                                                                  accaATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:5650/1‑150 (MQ=17)
                                                                                                                                    caATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2094423/1‑150 (MQ=32)
                                                                                                                                    caATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2535318/1‑150 (MQ=32)
                                                                                                                                     atagcaaaaaa‑GCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                  >  1:2693030/6‑150 (MQ=255)
                                                                                                                                     aaTGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGAAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTTCACGa                  >  2:3527676/1‑150 (MQ=12)
                                                                                                                                         tAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca              >  2:2871559/1‑150 (MQ=32)
                                                                                                                                          aaaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAACCTGCACGACccac             >  1:4071702/1‑150 (MQ=17)
                                                                                                                                           aaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  2:665264/1‑150 (MQ=32)
                                                                                                                                           aaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:155835/1‑150 (MQ=32)
                                                                                                                                            aaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:4992355/150‑1 (MQ=17)
                                                                                                                                            aaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:1634060/150‑1 (MQ=32)
                                                                                                                                             aaa‑GCGCCCTAATGGCTATTTTTTGCTATCTGCGATACTCAAATATTCTAACCTGATGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa          <  1:4116458/150‑1 (MQ=11)
                                                                                                                                               atgcgcCCTTAACTCGATTTTTTGCTTTCTGCGATCCTCAACGACACGAAACTGCAGCTGGTTTGAGTTGTTGAGAGCGTAACTACGTAGGCTCTCGCTGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         <  2:5456368/148‑1 (MQ=2)
                                                                                                                                                      |                                                                                                                                                     
TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781519‑781817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑