Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,173,431 T→A I97N (ATC→AAC)  ygiN → putative quinol monooxygenase YgiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,173,4310TA98.0% 186.7 / NA 51I97N (ATC→AAC) ygiNputative quinol monooxygenase YgiN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (21/29);  minor base C (1/0);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTT  >  NC_000913/3173191‑3173669
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
ggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcaggtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTg                                                                                                                                                                                                                                      <  2:352330‑M2/16‑1 (MQ=255)
        acggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATACCCGTATTCTGCAGCCAgg                                                                                                                                                                                                                              >  2:330222‑M2/228‑251 (MQ=255)
           gtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAt                                                                                                                                                                                                                           <  1:75956‑M2/27‑1 (MQ=255)
                           atggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGtt                                                                                                                                                                                                           <  2:217245‑M2/43‑1 (MQ=255)
                               agtcgtttaagtaggcgttttccgtaaattcatttagctcaaggcgctctgccatatcgctcattaattctgatgcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGGTTGcc                                                                                                                                                                                                       <  2:288675‑M2/47‑1 (MQ=255)
                               agtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGcc                                                                                                                                                                                                       >  2:161120‑M2/205‑251 (MQ=255)
                                           aggcgttttccggaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTg                                                                                                                                                                                            <  2:256356‑M2/58‑1 (MQ=255)
                                            ggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAc                                                                                                                                                                                          >  2:109618‑M2/192‑251 (MQ=255)
                                            ggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAc                                                                                                                                                                                          <  1:456739‑M2/60‑1 (MQ=255)
                                                     cgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATg                                                                                                                                                                                 >  1:320934‑M2/183‑251 (MQ=255)
                                                          atgtctgtagcgtaaggtgctctgccatctcgctcattaattctgattcttccaattaatcgtcagccttcaaaggggcaatattgccggagatacttcgttatctgccgcggcgagtaacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGTATTTAATCCTGCCTTGTGTGTCCGTCCATCCTGACCGGGCAATGTTCtt                                                                                                                                                                            <  2:375893‑M2/74‑1 (MQ=255)
                                                                gtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCttt                                                                                                                                                                      >  2:125807‑M2/172‑251 (MQ=255)
                                                                gtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                                                                       <  1:330222‑M2/79‑1 (MQ=255)
                                                                      caaggcgctctgccatatcgctcattaattctgattcctcatcttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                                                                <  2:8522‑M2/86‑1 (MQ=255)
                                                                         ggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggatgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCaa                                                                                                                                                             <  2:238589‑M2/89‑1 (MQ=255)
                                                                          gcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagctactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgatcaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAt                                                                                                                                                            >  2:145195‑M2/162‑251 (MQ=255)
                                                                           cgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtc                                                                                                                                                           <  1:302348‑M2/91‑1 (MQ=255)
                                                                           cgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtc                                                                                                                                                           <  1:224332‑M2/91‑1 (MQ=255)
                                                                            gctctgccataccgctcattaattctgattcctcaactaattcgccagccttcaaaggggcaatagtgccggagatgctacattatcggccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacgtctataaaggcatcaGAATAACCGGATTCTGCAGCCAGGGATTTAATCCTGCCTAGTTTGCCCGGCCATCCTGACCGGTCAATGGTCTTTCCCGGAAACCTCAAtct                                                                                                                                                          <  2:394059‑M2/92‑1 (MQ=255)
                                                                                           tcattaattctgattcctcaacttattcgccagccttcaaaggtgcaatatttccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCCATCTCCGCCATGTCGCCtt                                                                                                                                           >  2:179880‑M2/145‑251 (MQ=255)
                                                                                            cattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACAGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCttt                                                                                                                                          >  2:132532‑M2/144‑251 (MQ=255)
                                                                                                              aacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcg                                                                                                                        <  2:396876‑M2/126‑1 (MQ=255)
                                                                                                               acttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcgg                                                                                                                       <  1:140‑M2/127‑1 (MQ=255)
                                                                                                                   attcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGAt                                                                                                                   >  1:351156‑M2/121‑251 (MQ=255)
                                                                                                                       gccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTc                                                                                                               >  2:180626‑M2/117‑251 (MQ=255)
                                                                                                                                caaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGtt                                                                                                       <  1:368152‑M2/143‑1 (MQ=255)
                                                                                                                                caaaggggcaatattgccagagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                                                                      >  1:217129‑M2/108‑251 (MQ=255)
                                                                                                                                    ggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttcttc                                                                                                  <  1:285607‑M2/148‑1 (MQ=255)
                                                                                                                                     gggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                                                                 >  1:390146‑M2/103‑251 (MQ=255)
                                                                                                                                           tattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                                                                            >  1:262597‑M2/97‑250 (MQ=255)
                                                                                                                                               gccggagatactacattatctgccgcgtcgagtcccaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTATCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTGGCGGCGATCGTCCGAGCGTTTCTTCGCTAgcagca                                                                                        <  1:216433‑M2/158‑1 (MQ=255)
                                                                                                                                               gccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCat                                                                                       >  1:301084‑M2/93‑251 (MQ=255)
                                                                                                                                                ccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCata                                                                                      <  1:132532‑M2/160‑1 (MQ=255)
                                                                                                                                                    agatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccacgggacgtgctccgacggcgataaaggcatcaGAATAACCGTATGCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCAATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTAAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCACTCTTCCGAGCGTTTCTTCGCAAGCAGCATATcc                                                                                   >  1:177467‑M2/88‑250 (MQ=255)
                                                                                                                                                       tactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatc                                                                                >  1:457412‑M2/85‑250 (MQ=255)
                                                                                                                                                        actacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatcat                                                                              <  1:90113‑M2/168‑1 (MQ=255)
                                                                                                                                                              ttatctgacgcgtcgagtcacaggagatcaacgtgactagcgccgccggccgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCACTCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                        >  2:24843‑M2/78‑251 (MQ=255)
                                                                                                                                                              gtagctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCACCGCGTc                                                                        <  2:128115‑M2/174‑1 (MQ=255)
                                                                                                                                                               tatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                        <  2:286854‑M2/174‑1 (MQ=255)
                                                                                                                                                                    gccgcgtcgtgtcacaggagatcaacgtgaatagcgccgccggacgggctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAg                                                                      <  2:39858‑M2/176‑1 (MQ=255)
                                                                                                                                                                           cgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATc                                                            <  2:453742‑M2/186‑1 (MQ=255)
                                                                                                                                                                                 acaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTCCGCTGATCGTCttcc                                                     >  2:245041‑M2/59‑250 (MQ=255)
                                                                                                                                                                                 acaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat                                                    <  1:414705‑M2/193‑1 (MQ=255)
                                                                                                                                                                                          caacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATa                                            <  2:128248‑M2/202‑1 (MQ=255)
                                                                                                                                                                                           aacgtgaatagcgccgccggacgtgctccgacggcgataaaggcagcaGAATAACCGTCTTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAACCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAg                                           <  2:439286‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                               acgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGCCGGCg                        >  1:44112‑M2/29‑250 (MQ=255)
                                                                                                                                                                                                                    ctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGtt                  >  1:96816‑M2/24‑251 (MQ=255)
                                                                                                                                                                                                                       cgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCgg               <  1:108341‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                         acggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGAt             >  1:283348‑M2/19‑251 (MQ=255)
                                                                                                                                                                                                                                   aggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGt   >  1:127281‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                    ggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGtt  <  2:351156‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
GTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTT  >  NC_000913/3173191‑3173669

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