Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000964 165,826 +C intergenic (+1/‑4) trnI‑Asn → / → trnI‑Thr tRNA‑Asn/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000964165,8261.C97.5% 145.4 / ‑2.6 40intergenic (+1/‑4)trnI‑Asn/trnI‑ThrtRNA‑Asn/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (13/26);  total (14/26)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

ACACTCTCAATTTGTTTGGTGGCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGC‑A‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTA‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAG  >  NC_000964/165577‑166073
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
acacTCTCAATTTGTTTGGTGGCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGccat                                                                                                                                                                                                                                                         <  2:322001/251‑4 (MQ=255)
 cacTCTCAATTTGTTTGGTGGCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGccatt                                                                                                                                                                                                                                                        >  1:239948/1‑247 (MQ=255)
       cAATTTGTTTGGTGGCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGGGCCA‑TTGCtt                                                                                                                                                                                                                                                   <  2:170982/249‑1 (MQ=255)
                 ggtggCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTcagc                                                                                                                                                                                                                                       >  1:226097/1‑251 (MQ=255)
                    ggCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAgg                                                                                                                                                                                                                                    <  1:75642/251‑1 (MQ=255)
                                 agGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCa                                                                                                                                                                                                                       <  1:321115/251‑1 (MQ=255)
                                      acacCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGt                                                                                                                                                                                                                   >  1:138094/1‑250 (MQ=255)
                                        acCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAg                                                                                                                                                                                                                <  2:69648/251‑1 (MQ=255)
                                           cGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGaa                                                                                                                                                                                                             <  1:338043/251‑1 (MQ=255)
                                             ttCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGa                                                                                                                                                                                                              >  1:53767/1‑248 (MQ=255)
                                                         acacGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTGCCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc                                                                                                                                                                                               <  2:273455/251‑1 (MQ=255)
                                                                     aGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTgg                                                                                                                                                                                   <  2:265547/251‑1 (MQ=255)
                                                                                   gATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGTTCTGTAGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCAAAGCTCAGCAGGTCGAGCACTTCCCTTGAAAGGAAGAGGTCATCGGTTCGAGCCCGATTGGAAGCTTAAATggat                                                                                                                                                                     >  2:173709/1‑248 (MQ=37)
                                                                                      ggTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGtatta                                                                                                                                                                   <  1:238829/250‑1 (MQ=255)
                                                                                            cGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCaa                                                                                                                                                            <  2:203643/251‑1 (MQ=255)
                                                                                             gggggTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAg                                                                                                                                                           <  1:82499/251‑1 (MQ=255)
                                                                                                 gTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGttt                                                                                                                                                       >  1:40218/1‑251 (MQ=255)
                                                                                                           gtgAGAGTAGGACGCCTCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑AGTCCGCAGTAGCTCAGGGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAgg                                                                                                                                              <  2:20455/250‑1 (MQ=38)
                                                                                                                 gTAGGACGCCGCCAAGCAGATTTATTAATTCTTCCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGc                                                                                                                                       <  2:140253/251‑1 (MQ=255)
                                                                                                                 gTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGc                                                                                                                                       <  1:123443/251‑1 (MQ=255)
                                                                                                                    ggACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCttt                                                                                                                                    <  1:124970/251‑1 (MQ=255)
                                                                                                                     gACGCCGCCCAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTAGTAGGTTCGAATCCTACCTGCGGAGC‑CCTTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCtttt                                                                                                                                   >  2:274259/1‑251 (MQ=38)
                                                                                                                        gccgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGTAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCttttttt                                                                                                                                <  1:284614/251‑1 (MQ=255)
                                                                                                                              aaGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCg                                                                                                                          >  1:152068/1‑251 (MQ=255)
                                                                                                                               aGCAGATTTATTAATTCTTCCATGTTTTTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTTCCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTAATATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGa                                                                                                                         <  1:176034/251‑1 (MQ=38)
                                                                                                                                      ttattaATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGGAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTAGGAAGCATAAATGTATTATTACAAAAGTTTCTCTTAAGGAGAAAGATTTTTTTATTGCGATATACgg                                                                                                                  >  2:22828/1‑251 (MQ=37)
                                                                                                                                                              aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt                                                                                                                                                                                                        >  2:149286/1‑141 (MQ=21)
                                                                                                                                                              aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt                                                                                                                                                                                                        <  1:149286/141‑1 (MQ=21)
                                                                                                                                                              aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCttt                                                                                                                                    >  2:170523/1‑209 (MQ=255)
                                                                                                                                                              aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCttt                                                                                                                                    <  1:170523/209‑1 (MQ=255)
                                                                                                                                                                 aaGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTa                                                                                                                                                                            <  1:18887/166‑1 (MQ=255)
                                                                                                                                                                 aaGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTa                                                                                                                                                                            >  2:18887/1‑166 (MQ=255)
                                                                                                                                                                  aGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACTCCACCtt                                                                                      <  2:315332/251‑1 (MQ=255)
                                                                                                                                                                                  ccGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGt                                                                        <  1:188781/251‑1 (MQ=255)
                                                                                                                                                                                     cAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATTGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTcgc                                                                     <  1:141892/251‑1 (MQ=255)
                                                                                                                                                                                                     gagCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGCATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGAATTATTACCAAGGTTTCTCCTAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTc                                                     >  2:86778/1‑251 (MQ=255)
                                                                                                                                                                                                              gCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTcc                                            <  2:226097/251‑1 (MQ=255)
                                                                                                                                                                                                                          tCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTAAc                                >  1:212797/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggTCGAAGGTTCGAATCCTACCTACGGTGCCA‑TTGCTTCCATAGCTCATCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGGATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAAGCCCGTCTGCCGCTCCATTTATCCTAACgg                              <  2:166912/251‑1 (MQ=38)
                                                                                                                                                                                                                                         aaTCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTAACGGGATATTGTAaatt                 <  1:274259/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaTCCTACCTGCGGAGCCA‑TTGCTTCCATAGCTCAGCAGGTAGAGAACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGATTGGAAGCTTAAATGTATTATTACCAAGGATTATCCTAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGCAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTAACGGGATCTTGTAaatt                 >  2:262524/1‑251 (MQ=38)
                                                                                                                                                                                                                                          aTCCTCCCTGCGGAGCCA‑TTGCCTCCATAGCTCAGCAGGTAGAGCTCTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACCCCACCTTGCCAAGGTGGGGGTCTCGGGTTCGAATCCCTTCTTCCGCTCCATTTATCCTAACGGGATATTGTAaattg                <  1:261345/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gcc‑A‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg  <  1:284918/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
ACACTCTCAATTTGTTTGGTGGCGATAGCGAAGAGGTCACACCCGTTCCCATACCGAACACGGAAGTTAAGCTCTTCAGCGCCGATGGTAGTCGGGGGTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGC‑A‑TTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTA‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAG  >  NC_000964/165577‑166073

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑