Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,670 A→T 7.3% 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 7.3% 440.1 / 22.4 159intergenic (+18/‑415)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (66/81);  new base T (4/8);  total (70/89)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 170 total aligned reads displayed.

GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781521‑781818
                                                                                                                                                     |                                                                                                                                                    
 acctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                                     <  1:625761‑M1/95‑1 (MQ=255)
  cctcccttacttggtgggggtcggcggttcgatcccgtcatcacccaccacttgGGGTCGTTAGCTCTTTTGGTAGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                    >  1:4032224‑M1/55‑150 (MQ=255)
   tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                   <  2:1425429/150‑1 (MQ=255)
   ctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTCGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAg                                                                                                                                                   <  1:1796674‑M1/97‑1 (MQ=255)
   ctcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                   <  1:1796644‑M1/97‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAATCGCCCTAAAgg                                                                                                                                                  <  1:783791‑M1/98‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                  <  1:3043615‑M1/98‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                  <  1:776908‑M1/98‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAgg                                                                                                                                                  <  2:3088676‑M1/98‑1 (MQ=255)
    tcccttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                  <  2:1864322‑M1/98‑1 (MQ=255)
    cATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                  <  2:2798097/150‑1 (MQ=255)
     aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATGAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                 <  2:89305/150‑1 (MQ=255)
     aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                 <  2:2112925/150‑1 (MQ=255)
     aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                 <  2:75867/150‑1 (MQ=255)
      tCCTGCACGACCCTCCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                <  2:1450016/150‑1 (MQ=255)
        ttacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                              <  1:170972‑M1/102‑1 (MQ=255)
        cTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCt                                                                                                                                              >  1:1991168/1‑150 (MQ=255)
         tacaaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtt                                                                                                                                             >  2:1473159‑M1/48‑150 (MQ=255)
         tGCACGACCCACCATCCTTAATTATTAAGGCAGCATAATCCCGCACGGGGTCGTTAGCTCTGTTGGTTGAGCAGTTGACTTTTAATCAATGGGTCGCCGTTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtt                                                                                                                                             <  1:739452/150‑1 (MQ=18)
         tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                             <  1:2185479/150‑1 (MQ=255)
          gCACCACCCACCAACCTGAATGATTAAGGCAGCATAACCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTAAATCATTTGGTCGCCGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCCTTTAGGGCGCttt                                                                                                                                            <  1:3902350/150‑1 (MQ=18)
           caaggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                           >  1:2220508‑M1/46‑150 (MQ=255)
           cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                           >  2:3937524/1‑150 (MQ=255)
           cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                           >  2:929692/1‑150 (MQ=255)
           cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                           <  2:1991168/150‑1 (MQ=255)
            ccggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCACAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACCGGCGCttttt                                                                                                                                          >  2:3216561‑M1/45‑150 (MQ=255)
            caggagggggtcggcggttcgatcccgtcatcacccaccacttcGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATGGGTCGCAGGTTCGACTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                          >  1:709584‑M1/45‑150 (MQ=255)
            aCGACCCACCATCCTGACTGATTAAGGCAGCATACTCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                          <  1:541645/150‑1 (MQ=18)
             cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                         <  1:2187832/150‑1 (MQ=255)
              gACCCACCATCCTGAATTTTTATGTCAGCATAATCCCGCAAGGGGTCGTTATCTCAGTTGGTAGCGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                        <  1:3810821/150‑2 (MQ=17)
                  caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                    >  1:2006144/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:2560049/150‑1 (MQ=255)
                         cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAt                                                                                                                             <  1:1913728/150‑1 (MQ=255)
                         cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAt                                                                                                                             <  1:234283/150‑1 (MQ=255)
                          tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                            >  1:3084776/1‑150 (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)
                                    aGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                  <  2:1350427/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)
                                           aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCt                                                                                                           <  1:3299993/150‑1 (MQ=255)
                                             aaTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATAGGCGATACTAAAAGATTCGAACCTgc                                                                                                         >  1:29160/1‑150 (MQ=255)
                                              ccaccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAAGCGCCCGAAAGGCGCTTTTTTTCAATCTGCGATAATCAAAGATTCGAACCTgca                                                                                                        <  1:675491‑M1/140‑1 (MQ=255)
                                               tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                       <  1:2526240/150‑1 (MQ=39)
                                               caccacttcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTTCTATCTGCGATACTCAAATATTCGAACCTgcag                                                                                                       <  1:3388267‑M1/141‑1 (MQ=255)
                                               caccacttcGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAATAGCCCAAAAGGCGCTTTTTTTATATATGCGAAAGTCAAAATTTCCAACCGgcag                                                                                                       <  1:2994221‑M1/141‑1 (MQ=255)
                                                      aGGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                <  1:3346513/150‑1 (MQ=18)
                                                        gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                              <  1:1134429/150‑1 (MQ=17)
                                                           tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                           <  2:3315571/150‑1 (MQ=32)
                                                           tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCTGCAGGTttgagtt                                                                                           <  2:1756685/150‑1 (MQ=17)
                                                           tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                           <  1:1917615/150‑1 (MQ=17)
                                                             gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                         >  1:2149039/1‑150 (MQ=32)
                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAg                                                                                   >  2:2220012/1‑150 (MQ=17)
                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGATCGCAg                                                                                   >  2:2336167/1‑150 (MQ=17)
                                                                    cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                  <  1:1958845/150‑1 (MQ=32)
                                                                    cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                  <  1:2682952/150‑1 (MQ=32)
                                                                       ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                               >  2:1558787/1‑150 (MQ=32)
                                                                            agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                          <  2:1864027/150‑1 (MQ=32)
                                                                              agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                        <  2:3026323/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:2808310/1‑150 (MQ=32)
                                                                                                 aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTCTCTGGGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                     >  1:657098/1‑150 (MQ=17)
                                                                                                          aGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc                                            >  1:2668591/1‑150 (MQ=32)
                                                                                                             ttCGACTCCTGCTCGCCCCACCCCGGTCACCACGCGCCCTACAGGCGCTTTTTGGCTCTCTGCGATCCTCAACGATTCGACCCTGCAGCTGGTGTGAGTGGAGGGCCGCGACACTCCGGAGCCGCGCGCGGCGACGGCCCGAAGGgcgag                                         <  2:812116/150‑1 (MQ=2)
                                                                                                             ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         >  1:1030020/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)
                                                                                                                               caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATAGTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                       >  2:2476110/1‑150 (MQ=17)
                                                                                                                                accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                      >  2:3307784/1‑150 (MQ=32)
                                                                                                                                accaATGTAAAAAAGCGCCCTAAAGGCGCTGTTTTGCTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                      >  2:2769708/1‑150 (MQ=17)
                                                                                                                                       caaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATATGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCCAAACAACGGAGCCGCACGCGGCGACGGCCCCAAGGGCGAGCCAAGCGAGTCATCCAGCcccaccca               >  1:3264670/2‑150 (MQ=11)
                                                                                                                                       aaaaaaaGCGCCCTAACGGCTCTTTTTTGCTTTCTGCCATACTCACAGACTCGAACCTGCAGCGGTTTTGCTTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca               <  1:3375469/149‑1 (MQ=11)
                                                                                                                                        aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac              >  2:3792869/1‑150 (MQ=32)
                                                                                                                                            aaGCGCCCTAACGGCGCTTTTTTGCTCTCTGCGATACTATAAGATTCGAACGTGCTGCAGTTTTGAGTTGAGCGCAGCTATACAACTGTGGCGCTGGCTGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt          <  2:2027184/150‑1 (MQ=11)
                                                                                                                                                  ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa    >  2:2509870/1‑150 (MQ=32)
                                                                                                                                                    tAAAGGCTTTTTTTTGCTTTCGGCTATACTCAAATTTTGGAACGTGCAGCTGTTTTGTTTTGCGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa  <  1:1558787/150‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                                                    
GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/781521‑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: