Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,670 A→T 5.1% intergenic (+18/‑415) 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,6700AT 5.1% 722.8 / 27.6 237intergenic (+18/‑415)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (124/97);  new base T (10/2);  total (135/102)
Fisher's exact test for biased strand distribution p-value = 7.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

Only 100 of 244 total aligned reads displayed.

AGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781522‑781818
                                                                                                                                                    |                                                                                                                                                    
acctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                                     <  2:1566098‑M1/95‑1 (MQ=255)
acctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                                     <  2:3104501‑M1/95‑1 (MQ=255)
acctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                                     <  2:1006350‑M1/95‑1 (MQ=255)
aGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAA                                                                                                                                                     >  1:1917591‑M1/1‑150 (MQ=255)
 cctcccttacatggagggggtcggcggttcgatcccgtcatcacccaccacttgGGGTCTTTAGCTCAGTTGGTAGAGCTGTTGACTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                    >  1:1402431‑M1/55‑150 (MQ=255)
 cctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                    >  2:1266894‑M1/55‑150 (MQ=255)
  tCATCCTGCACGACCCACCATCCTTAATTATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCTGTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCCCTAAAG                                                                                                                                                   <  1:869226‑M1/150‑1 (MQ=255)
  tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAG                                                                                                                                                   <  2:3468492‑M1/150‑1 (MQ=255)
     tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAATCGCCCTAAAGGCG                                                                                                                                                >  1:133093‑M1/1‑150 (MQ=255)
     tCCTGCACGACCCACCAACCTGAATGAGAAACGCAGCATAAACCCTCAAGGGTTCGTAAGCTCAGTTGGTAGAGCAGTTGACATTTAATCAATAGGTCGCAGGATCGAAACCTTCACGACCCACCAATGTAAAAAAGCGCCCTATAGgcg                                                                                                                                                >  2:2816753/1‑150 (MQ=17)
     ccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                <  2:1881023‑M1/100‑1 (MQ=255)
      cttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                               >  2:2420897‑M1/50‑150 (MQ=255)
       ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                              <  1:2248657‑M1/102‑1 (MQ=255)
           aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                          >  1:2265760‑M1/1‑150 (MQ=255)
             gACCCACCATCCTTAATGATTACGGCAGCATAATCCCGCTTGGGGTGGTTTGCTCTGTTGTTAGTGCTGTTGACTTTTAATCTTTTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                        <  1:711946/150‑2 (MQ=17)
             gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTG                                                                                                                                        >  1:3976726‑M1/1‑150 (MQ=255)
              aCCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGAAGGTTCTAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGC                                                                                                                                       <  2:3286498‑M1/150‑1 (MQ=255)
                 caccaTCCTTAATGATTAAGGAAGCTTAAGCCCCGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAT                                                                                                                                    <  1:3073243‑M1/150‑1 (MQ=255)
                  accaTCCTGTATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCCGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATc                                                                                                                                   <  1:1873091/150‑1 (MQ=255)
                   acaTCCTGAATGATTAATGCAGCATATTCCCGCATGGGGTCTTTAGCTAAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                  <  1:2298426/149‑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)
                               aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCGGCGATACTCaac                                                                                                                      >  2:2621204/1‑149 (MQ=255)
                               aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATgcaggcatcaaaaa                                                                                                                      >  2:1581212‑M1/1‑136 (MQ=255)
                                 tAAGGCAGCATACTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGa                                                                                                                    <  2:1936550/150‑1 (MQ=255)
                                        gcaTAATCCCGCCAGGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                             <  1:813411/150‑1 (MQ=255)
                                                  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                   >  2:752396/1‑150 (MQ=35)
                                                     aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                <  1:2812544/150‑1 (MQ=33)
                                                       gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                              <  2:1625983/150‑1 (MQ=32)
                                                       gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                              <  2:1181263/150‑1 (MQ=32)
                                                           cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg                                                                                          <  1:3823218/150‑1 (MQ=32)
                                                            gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTACTCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                         <  2:3914020/150‑1 (MQ=17)
                                                             gtAGCTCAGTGGGTAGCGCAGTTGACTTGTAATCTATGGGGCGCAGGTTCGAATCCTGCACGACCCACCAAGGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                        <  2:2184005/149‑1 (MQ=11)
                                                                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                    >  2:3343200/1‑150 (MQ=32)
                                                                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                    <  1:2939113/150‑1 (MQ=32)
                                                                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                   <  1:1130399/150‑1 (MQ=32)
                                                                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                   <  2:2234059/150‑1 (MQ=32)
                                                                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                   >  1:3495449/1‑150 (MQ=32)
                                                                       tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                              <  1:100831/150‑1 (MQ=32)
                                                                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                             >  2:1905051/1‑150 (MQ=32)
                                                                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                             >  2:1936002/1‑150 (MQ=17)
                                                                              tcAGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                       <  2:649537/149‑1 (MQ=17)
                                                                              gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                       >  1:4011027/1‑150 (MQ=17)
                                                                              gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAACTCAAAGATTAGAAACTGCAGCAGGTTTGATATTAGCGCAGAGAAACAACAGa                                                                       >  1:2022156/1‑150 (MQ=11)
                                                                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                 >  1:250784/1‑150 (MQ=32)
                                                                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                 >  1:797510/1‑150 (MQ=32)
                                                                                     aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGCTACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                                                                >  1:3116736/1‑150 (MQ=17)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAACCAACGGAGCCGCTcgc                                                              >  2:457177/1‑150 (MQ=17)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  1:335600/1‑150 (MQ=17)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACACAAAGATACGAACCTGAACCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  1:4148115/1‑150 (MQ=11)
                                                                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                        >  1:1438062/1‑150 (MQ=32)
                                                                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                        >  1:1068126/1‑150 (MQ=32)
                                                                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGAAACACAAATATTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGAGACAACGGAGCCGCTAGCGTCGAc                                                        >  2:4116163/1‑150 (MQ=11)
                                                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                      <  2:1327029/150‑1 (MQ=32)
                                                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGATATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                      <  2:3699353/150‑1 (MQ=17)
                                                                                                     tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTAGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                >  2:3387034/1‑150 (MQ=17)
                                                                                                        cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                             <  2:2276884/150‑1 (MQ=17)
                                                                                                          ggTTCGACTCCTGCCCGACCCACCAATGTAAAAAAGGGCCCTAAAGGCGCTTTTTTGCTATCTGCTATACTCATAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGTAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                           <  2:2359318/150‑1 (MQ=11)
                                                                                                            ttCGACTCCTGCACGCCCCACCCATGTAAAAACGCGCCCTCCCGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         <  2:1917632/150‑1 (MQ=11)
                                                                                                            ttCGAATCCTGCACGACCCACCATTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         <  2:3057072/150‑1 (MQ=17)
                                                                                                            ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         <  2:2067346/150‑1 (MQ=32)
                                                                                                            tgcGAATCCTGCACGACCCACCTATTTAAAAACGCGCCCTAAAGGCGATTTTTTGCTATATGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCTGCGACGGCCCGAAGGgcgag                                         <  2:2477464/148‑1 (MQ=11)
                                                                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAACTCAAAGATAAGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                     >  2:1409741/1‑150 (MQ=255)
                                                                                                                 aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacg                                    >  2:326919‑M1/1‑54 (MQ=255)
                                                                                                                  tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgc                                   >  2:3228281‑M1/1‑53 (MQ=255)
                                                                                                                    cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcat                                 >  1:1151334‑M1/1‑51 (MQ=255)
                                                                                                                      gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcacctctatgtcggtggaattaggcgtaaaatgacgcatcc                               >  1:547393‑M1/1‑49 (MQ=255)
                                                                                                                       cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGATTGAGATGAGCGCAGCGAAACAACGCACCCGCACGCGGCGACGGCCCGAAGGGCGACCGAAGCCAgac                              >  2:2712995/1‑148 (MQ=11)
                                                                                                                        aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCATTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctg                             <  2:133099‑M1/150‑104 (MQ=255)
                                                                                                                          gACCCACCCATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcacctctctggcggtggacttaggcgtaaaatgacgcatcctgca                           >  2:8948‑M1/1‑45 (MQ=255)
                                                                                                                          gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCATc                           >  1:2896384/1‑150 (MQ=17)
                                                                                                                             ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacat                        >  2:318861‑M1/1‑42 (MQ=255)
                                                                                                                             ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACGTGCAGCAGGTTTGAGTTGAGCGCAGCGATACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                        >  2:1213800/1‑150 (MQ=17)
                                                                                                                             ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGATCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGGCATCCTg                        >  1:1263686/1‑150 (MQ=17)
                                                                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacctta                      >  2:2393802‑M1/1‑40 (MQ=255)
                                                                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacatta                      >  2:876263‑M1/1‑40 (MQ=255)
                                                                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacatta                      <  2:465347‑M1/150‑111 (MQ=255)
                                                                                                                                 caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattagg                    <  1:1180044‑M1/150‑113 (MQ=255)
                                                                                                                                 caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattagg                    <  1:1908355‑M1/150‑113 (MQ=255)
                                                                                                                                  actgtAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggc                   <  1:2093773‑M1/148‑114 (MQ=255)
                                                                                                                                     gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgta                >  2:2677644‑M1/1‑34 (MQ=255)
                                                                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatagtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgatgcatcctgcacattaggcgtaa               >  1:3931232‑M1/1‑33 (MQ=255)
                                                                                                                                       aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat              >  2:3527394‑M1/1‑32 (MQ=255)
                                                                                                                                       aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATtaaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat              >  1:92976‑M1/1‑32 (MQ=255)
                                                                                                                                       aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaat              <  2:1390143‑M1/150‑119 (MQ=255)
                                                                                                                                         aaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca            >  1:2988487/1‑150 (MQ=12)
                                                                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattc            >  2:4003504‑M1/1‑30 (MQ=255)
                                                                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattc            >  1:2568710‑M1/1‑30 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga          <  2:3402053‑M1/150‑123 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga          <  2:1701345‑M1/150‑123 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga          <  2:805274‑M1/150‑123 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga          <  2:571005‑M1/150‑123 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcga          <  2:4011027‑M1/150‑123 (MQ=255)
                                                                                                                                             gcgcCCTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagt        >  2:2078280‑M1/1‑26 (MQ=255)
                                                                                                                                               gcCCTAATGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtga      <  2:1433421‑M1/150‑127 (MQ=255)
                                                                                                                                                cccTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattaggcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgac     <  2:110308‑M1/150‑128 (MQ=255)
                                                                                                                                                  cTAAAGGCGCTTTTTTGCTATtcaggcatcctcaatttcactttgtaaacctgatgacatcgtcagagcttactgtgcaagcaactctatgtcggtggaattatgcgtaaaatgacgcatcctgcacattaggcgtaattcgagtgactt   >  1:3716493‑M1/1‑21 (MQ=255)
                                                                                                                                                  cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTAGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa   >  1:429789/1‑150 (MQ=17)
                                                                                                                                                    |                                                                                                                                                    
AGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781522‑781818

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: