Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 7.0% 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 7.0% 443.4 / 19.7 159intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (65/77);  new base G (4/7);  total (70/89)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

Only 100 of 167 total aligned reads displayed.

GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCG  >  NC_000913/781523‑781821
                                                                                                                                                     |                                                                                                                                                     
gTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTATTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                       <  1:654580/150‑1 (MQ=255)
cctcccttacttggtgggggtcggcggttcgatcccgtcatcacccaccacttgGGGTCGTTAGCTCTTTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                       >  1:4032224‑M1/55‑150 (MQ=255)
cctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                       <  1:2727626‑M1/96‑1 (MQ=255)
cctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                       <  1:1981218‑M1/96‑1 (MQ=255)
 tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                      <  2:1425429/150‑1 (MQ=255)
 ctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTCGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAg                                                                                                                                                      <  1:1796674‑M1/97‑1 (MQ=255)
  tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                     <  1:776908‑M1/98‑1 (MQ=255)
  tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAgg                                                                                                                                                     <  2:3088676‑M1/98‑1 (MQ=255)
  cATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                     <  2:2798097/150‑1 (MQ=255)
  cATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                     >  1:3046438/1‑150 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATGAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  2:89305/150‑1 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  2:3181051/150‑1 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  2:75867/150‑1 (MQ=255)
   aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                    <  2:2112925/150‑1 (MQ=255)
    tCCTGCACGACCCTCCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  2:1450016/150‑1 (MQ=255)
    ccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                   <  2:1514703‑M1/100‑1 (MQ=255)
      ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  1:170972‑M1/102‑1 (MQ=255)
      ttacaaggagggggtcggcgggtcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGGTGGTAGAGCAGGTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                 <  2:1792200‑M1/102‑1 (MQ=255)
      cTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCt                                                                                                                                                 >  1:1991168/1‑150 (MQ=255)
       tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                                <  1:2185479/150‑1 (MQ=255)
        gCACCACCCACCAACCTGAATGATTAAGGCAGCATAACCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTAAATCATTTGGTCGCCGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCCTTTAGGGCGCttt                                                                                                                                               <  1:3902350/150‑1 (MQ=18)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                              >  2:929692/1‑150 (MQ=255)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                              >  2:3937524/1‑150 (MQ=255)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                              <  2:1991168/150‑1 (MQ=255)
          aaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                             >  1:303543‑M1/45‑150 (MQ=255)
           cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                            <  1:2187832/150‑1 (MQ=255)
            gACCCACCATCCTGAATTTTTATGTCAGCATAATCCCGCAAGGGGTCGTTATCTCAGTTGGTAGCGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                           <  1:3810821/150‑2 (MQ=17)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCtct                                                                                                                                       >  1:2006205/1‑148 (MQ=255)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                       >  1:1789092/1‑150 (MQ=255)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTGCTAt                                                                                                                                       >  1:1789083/1‑150 (MQ=255)
                  ccttCCAGAACGATAACGGCTGCATAACCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                     <  2:3669213/147‑1 (MQ=14)
                  ccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                     <  2:2633758/150‑1 (MQ=255)
                     tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                                                                                                                                  <  1:2559932/150‑1 (MQ=255)
                     cggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGGTCAGTTGGTAGAGCAGTTGACTTTTAATGAATTTGTCGCAGGTTGGAATCCTGCACGACCCACCAATGTAAATAAGCGCCCTAAAGGGGGTTTTTTTATatcttcg                                                                                                                                  <  1:1876438‑M1/117‑4 (MQ=255)
                       cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAt                                                                                                                                <  1:234283/150‑1 (MQ=255)
                           aTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGCTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTc                                                                                                                            <  2:1915716/150‑1 (MQ=255)
                             gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaa                                                                                                                          <  1:4104881/150‑1 (MQ=255)
                             gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaa                                                                                                                          >  1:383007/1‑150 (MQ=255)
                              aatAAGGCAGCATAAGCCCGCAAGGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                         <  1:1243911/148‑1 (MQ=255)
                               ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAg                                                                                                                        >  1:3849703/1‑150 (MQ=255)
                               atcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCGAAg                                                                                                                        <  2:317388‑M1/127‑1 (MQ=255)
                                  aGGCAGCATAATCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                     <  1:2458994/150‑1 (MQ=255)
                                  aGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                     <  2:2143417/150‑1 (MQ=255)
                                    gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                                   >  1:1267899/1‑150 (MQ=255)
                                     cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                                  <  1:1840589/150‑1 (MQ=255)
                                      cgcaTTATCCCTATAGGGGTCTTTAGCTCAGTTGGTAGATCTTTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                                 <  1:3595415/149‑1 (MQ=11)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                                 <  1:3595389/150‑1 (MQ=255)
                                        tcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTTCTATATGCGATACTCAAAGATTCGAAcc                                                                                                               <  1:3900464‑M1/136‑1 (MQ=255)
                                        caTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAcc                                                                                                               >  2:3968659/1‑150 (MQ=255)
                                         aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCt                                                                                                              >  1:2741739/1‑150 (MQ=255)
                                            ccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAAGCGCCCGAAAGGCGCTTTTTTTCAATCTGCGATAATCAAAGATTCGAACCTgca                                                                                                           <  1:675491‑M1/140‑1 (MQ=255)
                                            aTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                                           <  2:1365161/150‑1 (MQ=40)
                                             caccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTTCTATCTGCGATACTCAAATATTCGAACCTgcag                                                                                                          <  1:3388267‑M1/141‑1 (MQ=255)
                                             caccacttcGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAATAGCCCAAAAGGCGCTTTTTTTATATATGCGAAAGTCAAAATTTCCAACCGgcag                                                                                                          <  1:2994221‑M1/141‑1 (MQ=255)
                                                    aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                   <  1:3265995/150‑1 (MQ=33)
                                                      gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                                 <  1:1134429/150‑1 (MQ=17)
                                                         tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                              <  1:716030/150‑1 (MQ=17)
                                                         tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                              <  1:1917615/150‑1 (MQ=17)
                                                         tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                              <  2:1111456/150‑1 (MQ=17)
                                                           gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGCTACAAAAATATTAGAAACTGCAGCAGGTttgatatca                                                                                            >  1:2664092/1‑145 (MQ=11)
                                                            ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                           <  1:756399/150‑1 (MQ=32)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAg                                                                                      >  2:2220012/1‑150 (MQ=17)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGATCGCAg                                                                                      >  2:2336167/1‑150 (MQ=17)
                                                                   aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                    >  2:3375038/1‑150 (MQ=32)
                                                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                                  >  2:1558787/1‑150 (MQ=32)
                                                                          agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                             <  2:1864027/150‑1 (MQ=32)
                                                                                 ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                      >  2:537447/1‑150 (MQ=32)
                                                                                     cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                                  >  1:3650805/1‑150 (MQ=32)
                                                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                 >  2:1830447/1‑150 (MQ=32)
                                                                                      ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                 >  1:1029101/1‑150 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                           >  2:328606/1‑150 (MQ=17)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGAAACACAAATATACGAGCCAGCAGCAGCATTGAGATGAGCTCAGCGAAAACAAGGAGCCGCTCGCGGCGAc                                                           >  1:1079859/1‑150 (MQ=11)
                                                                                             cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg                                                          >  1:2808406/1‑150 (MQ=32)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc                                               >  1:2668591/1‑150 (MQ=32)
                                                                                                          gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGTGCGATACACAAATATACGAACATGCGGCAGGATTGAGTTGAGAGAACCGAAACAACTGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                             >  1:2802391/1‑150 (MQ=11)
                                                                                                           ttCGACTCCTGCTCGCCCCACCCCGGTCACCACGCGCCCTACAGGCGCTTTTTGGCTCTCTGCGATCCTCAACGATTCGACCCTGCAGCTGGTGTGAGTGGAGGGCCGCGACACTCCGGAGCCGCGCGCGGCGACGGCCCGAAGGgcgag                                            <  2:812116/150‑1 (MQ=2)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                            >  1:1030020/1‑150 (MQ=32)
                                                                                                                aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                       >  1:669745/1‑150 (MQ=32)
                                                                                                                aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                       >  1:1093013/1‑150 (MQ=32)
                                                                                                                 tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                      >  2:2892466/1‑150 (MQ=32)
                                                                                                                 tCCGGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                      <  2:1917987/150‑1 (MQ=255)
                                                                                                                 accTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                      >  1:3454439/2‑150 (MQ=17)
                                                                                                                   cTGCACGACCCTCCAATGTAAAAAAGCGCCCTTAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGTTTTGAGTTGTGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                    <  2:2074298/150‑1 (MQ=11)
                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                    <  2:272325/150‑1 (MQ=32)
                                                                                                                    tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                   <  2:3659653/150‑1 (MQ=32)
                                                                                                                       aCGACCCTCCACTGTCTTATAGCGCCGTATAGGCGCTTTTTTGCTATCTGCGATTCTCAAAGTTTCGAACATGCAGCAGGTTTGAGTTGAGCGCAGCGTAACTACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                                <  2:1699205/150‑1 (MQ=11)
                                                                                                                            ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACACCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                           >  2:2052895/1‑150 (MQ=17)
                                                                                                                             caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATAGTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                          >  2:2476110/1‑150 (MQ=17)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                         >  2:3307784/1‑150 (MQ=32)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                         >  2:3307788/1‑150 (MQ=32)
                                                                                                                               ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATGCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGACGCGAGTCCTCCGGCAc                        >  2:2733670/1‑150 (MQ=11)
                                                                                                                                     caaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATATGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCCAAACAACGGAGCCGCACGCGGCGACGGCCCCAAGGGCGAGCCAAGCGAGTCATCCAGCcccaccca                  >  1:3264670/2‑150 (MQ=11)
                                                                                                                                     aaaaaaaGCGCCCTAACGGCTCTTTTTTGCTTTCTGCCATACTCACAGACTCGAACCTGCAGCGGTTTTGCTTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                  <  1:3375469/149‑1 (MQ=11)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                 >  2:3792869/1‑150 (MQ=32)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                 >  1:477442/1‑150 (MQ=32)
                                                                                                                                       aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                >  1:2163830/1‑150 (MQ=32)
                                                                                                                                         aaaGCCCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa              <  2:3018401/150‑1 (MQ=17)
                                                                                                                                                ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa       >  2:2509870/1‑150 (MQ=32)
                                                                                                                                                  tAACGGCGCTTTTTTGCTTTCTGCTATACTAAAATTTTCGACCGTGCAGCTGGTTTGAGTTGTGCGCAGGGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     <  1:3317940/150‑1 (MQ=11)
                                                                                                                                                  tAAAGGCTTTTTTTTGCTTTCGGCTATACTCAAATTTTGGAACGTGCAGCTGTTTTGTTTTGCGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa     <  1:1558787/150‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
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: