Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 5.1% intergenic (+20/‑413) 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,6720AG 5.1% 695.8 / 21.6 235intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (123/96);  new base G (10/2);  total (133/102)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 245 total aligned reads displayed.

GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCG  >  NC_000913/781523‑781821
                                                                                                                                                     |                                                                                                                                                     
 tCATCCTGCACGACCCACCATCCTTAATTATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCTGTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCCCTAAAG                                                                                                                                                      <  1:869226‑M1/150‑1 (MQ=255)
 tCATCCTGCACGACCCACCATCCTGTATGATTAAGGCAGCATAATCCCGATAGGGGTCGTTTGCTCTGTTGGTTTAGCATTTGTCTTTTAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAG                                                                                                                                                      <  1:859657‑M1/150‑1 (MQ=255)
  cATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGG                                                                                                                                                     >  2:3823218‑M1/1‑150 (MQ=255)
    tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAATCGCCCTAAAGGCG                                                                                                                                                   >  1:133093‑M1/1‑150 (MQ=255)
    ccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  2:1881023‑M1/100‑1 (MQ=255)
    ccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  2:607877‑M1/100‑1 (MQ=255)
    ccttacaaggaggggggcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  2:3387684‑M1/100‑1 (MQ=255)
     cttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                                  >  2:116049‑M1/50‑150 (MQ=255)
     cttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                                  >  2:2420897‑M1/50‑150 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  1:2248657‑M1/102‑1 (MQ=255)
        acaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                               <  2:514263‑M1/104‑1 (MQ=255)
          aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                             <  2:3423140‑M1/150‑1 (MQ=255)
          aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                             >  1:2265760‑M1/1‑150 (MQ=255)
            gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTG                                                                                                                                           >  1:3976726‑M1/1‑150 (MQ=255)
                caccaTCCTTAATGATTAAGGAAGCTTAAGCCCCGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAT                                                                                                                                       <  1:3073243‑M1/150‑1 (MQ=255)
                  ccaTCCTGAATGATTAAGGCATCATAAGGCCGCTTGGGGTCTTTAGCTCAGTGGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                     <  1:2090834/150‑1 (MQ=17)
                     tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTGCTATCTGCGa                                                                                                                                 <  1:122045/150‑1 (MQ=255)
                      ccTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                                 <  1:1332306/150‑1 (MQ=255)
                        tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                               <  1:3410550/150‑1 (MQ=255)
                        tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                               >  2:256512/1‑150 (MQ=255)
                             gATTAAGGCAGCATAATCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaa                                                                                                                          <  2:4179531/150‑1 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCGGCGATACTCaac                                                                                                                         >  2:2621204/1‑149 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                         >  2:2812544/1‑150 (MQ=255)
                                   ggCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                    <  2:144751/150‑1 (MQ=255)
                                               ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCGATCTGCGATACTCAAAGATTCGAACCTgcagcc                                                                                                        >  2:1314419/1‑149 (MQ=38)
                                                 gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                      >  2:4062271/1‑150 (MQ=35)
                                                 gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                      >  2:3269927/1‑150 (MQ=35)
                                                    aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                   <  1:2812544/150‑1 (MQ=33)
                                                      gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                                 <  2:1625983/150‑1 (MQ=32)
                                                          cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg                                                                                             <  1:3823218/150‑1 (MQ=32)
                                                           gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                            >  1:190506/1‑150 (MQ=32)
                                                            gtAGCTCAGTGGGTAGCGCAGTTGACTTGTAATCTATGGGGCGCAGGTTCGAATCCTGCACGACCCACCAAGGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                           <  2:2184005/149‑1 (MQ=11)
                                                              aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                                                                         <  2:3109115/150‑1 (MQ=32)
                                                                cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                       >  2:3343200/1‑150 (MQ=32)
                                                                cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                       <  1:2939113/150‑1 (MQ=32)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                      <  2:2234059/150‑1 (MQ=32)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                      <  1:1130399/150‑1 (MQ=32)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                      >  1:2684490/1‑150 (MQ=17)
                                                                      tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                 <  1:100831/150‑1 (MQ=32)
                                                                       ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTTCAGCAGGTTTGAGATGAGCGCAGCGAaac                                                                                >  2:1934413/1‑150 (MQ=17)
                                                                       ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                >  2:1905051/1‑150 (MQ=32)
                                                                       ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                >  2:1936002/1‑150 (MQ=17)
                                                                         tAGAGCAGTTGACTTTTAATCACTCGGTCGCAGGTTCGAGTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgc                                                                              <  2:3976726‑M1/150‑58 (MQ=255)
                                                                          agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                             <  2:3256224/150‑1 (MQ=32)
                                                                              cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                                                                         >  1:3369767/1‑150 (MQ=32)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                       >  2:1419961/1‑150 (MQ=32)
                                                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                      >  1:3520118/1‑150 (MQ=32)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                    >  1:797510/1‑150 (MQ=32)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                    >  2:3593972/1‑150 (MQ=32)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                    >  1:250784/1‑150 (MQ=32)
                                                                                    aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGCTACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                                                                   >  1:3116736/1‑150 (MQ=17)
                                                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACACAAAGATACGAACCTGAACCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                 >  1:4148115/1‑150 (MQ=11)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           >  1:1068126/1‑150 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           >  1:1438062/1‑150 (MQ=32)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                          <  1:738606/101‑1 (MQ=32)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                          >  2:738606/1‑101 (MQ=32)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGATATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                         <  2:3699353/150‑1 (MQ=17)
                                                                                                       cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                <  2:2276884/150‑1 (MQ=17)
                                                                                                       cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                <  2:4009702/150‑1 (MQ=32)
                                                                                                         ggTTCGACTCCTGCCCGACCCACCAATGTAAAAAAGGGCCCTAAAGGCGCTTTTTTGCTATCTGCTATACTCATAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGTAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                              <  2:2359318/150‑1 (MQ=11)
                                                                                                           ttCGAATCCTGCACGACCCACCATTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                            <  2:3057072/150‑1 (MQ=17)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                            <  2:2067346/150‑1 (MQ=32)
                                                                                                                 tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgc                                      >  2:3484297‑M1/1‑53 (MQ=255)
                                                                                                                 tCCTGCACGACCCACCAATGTAAAAAAGCGCCATAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                      <  2:1237932/150‑1 (MQ=17)
                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcat                                    >  1:1151334‑M1/1‑51 (MQ=255)
                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                    <  2:250784/150‑1 (MQ=32)
                                                                                                                    tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatc                                   <  1:409660‑M1/150‑101 (MQ=255)
                                                                                                                     gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcacctctatgtcggtggaattaggcgtaaaatgacgcatcc                                  >  1:547393‑M1/1‑49 (MQ=255)
                                                                                                                      cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGATTGAGATGAGCGCAGCGAAACAACGCACCCGCACGCGGCGACGGCCCGAAGGGCGACCGAAGCCAgac                                 >  2:2712995/1‑148 (MQ=11)
                                                                                                                       aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttttaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctg                                <  2:133093‑M1/150‑104 (MQ=255)
                                                                                                                         gACCCACCCATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcacctctctggcggtggacttaggcgtaaaatgacgcatcctgca                              >  2:8948‑M1/1‑45 (MQ=255)
                                                                                                                         gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCATc                              >  1:2896384/1‑150 (MQ=17)
                                                                                                                            ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAATGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTg                           >  2:1213749/1‑150 (MQ=17)
                                                                                                                            ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGATCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGGCATCCTg                           >  1:1263686/1‑150 (MQ=17)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacatta                         <  2:465311‑M1/150‑111 (MQ=255)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacatta                         >  2:876263‑M1/1‑40 (MQ=255)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactcgatgtcggtggaattaggcgtaaaatgacgcatcctgcacatta                         >  2:2393936‑M1/1‑40 (MQ=255)
                                                                                                                                caATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattagg                       <  1:1409741‑M1/150‑113 (MQ=255)
                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattagg                       <  1:1912375‑M1/150‑113 (MQ=255)
                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattagg                       <  1:1180044‑M1/150‑113 (MQ=255)
                                                                                                                                    gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgta                   >  2:2677644‑M1/1‑34 (MQ=255)
                                                                                                                                     tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatagtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgatgcatcctgcacattaggcgtaa                  >  1:3931232‑M1/1‑33 (MQ=255)
                                                                                                                                     caaaaaagCGCCTTTAGGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaa                  >  2:3323412‑M1/2‑33 (MQ=255)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATtaaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat                 >  1:92976‑M1/1‑32 (MQ=255)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat                 <  2:1390104‑M1/150‑119 (MQ=255)
                                                                                                                                        aaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca               >  1:2988487/1‑150 (MQ=12)
                                                                                                                                          aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga             <  2:1701345‑M1/150‑123 (MQ=255)
                                                                                                                                          aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga             <  2:805274‑M1/150‑123 (MQ=255)
                                                                                                                                           aGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgag            <  2:200180‑M1/150‑124 (MQ=255)
                                                                                                                                           aGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgag            <  2:3197824‑M1/150‑124 (MQ=255)
                                                                                                                                            gcgcCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagt           >  2:2078280‑M1/1‑26 (MQ=255)
                                                                                                                                                ccTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgact       >  2:140512‑M1/1‑22 (MQ=255)
                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattatgcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgactt      >  1:3716493‑M1/1‑21 (MQ=255)
                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATATGCGATACACAAATATTCGAAACTGCAGCAGGTTTTAGTTTAGCGCAGCGAAACAACGGAGCCGCTCGCCGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCAGCACGACCCACCAAAGTaaaaa      >  1:429835/1‑150 (MQ=11)
                                                                                                                                                 cTAAAGGCGCTATTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGCTTTGAGTTGAGCGCAGCGAAACAAAGCAGCCGCTCGCGGCGACGGCACGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa      >  2:1065951/1‑150 (MQ=11)
                                                                                                                                                  tAAAGGCGCTTTTTTGCTATATGCTATACTCAAATATTCGATCCTGCAGCAGGTTTGTGTTGAGCGCTGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     <  1:2312397/150‑1 (MQ=11)
                                                                                                                                                    aaGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttc   >  2:1605590‑M1/1‑18 (MQ=255)
                                                                                                                                                    aaGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttc   >  2:2668401‑M1/1‑18 (MQ=255)
                                                                                                                                                     aGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttcc  <  1:1735970‑M1/150‑134 (MQ=255)
                                                                                                                                                     aGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgacttttcc  <  1:2124‑M1/150‑134 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCG  >  NC_000913/781523‑781821

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: