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.2% 185.9 / NA 55I97N (ATC→AAC) ygiNputative quinol monooxygenase YgiN
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (24/29);  minor base C (1/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGT  >  NC_000913/3173189‑3173668
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
acggcaacgcacggtccaggatcacgtacatngagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTc                                                                                                                                                                                                                                       <  2:259834‑M2/14‑1 (MQ=255)
  ggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTg                                                                                                                                                                                                                                     >  2:183827‑M2/236‑251 (MQ=255)
     aacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgatcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGc                                                                                                                                                                                                                                    >  1:211399‑M2/233‑249 (MQ=255)
          acggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgcaggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcataaGAATAACCGTATTCTGCAGCCAgg                                                                                                                                                                                                                             >  1:21186‑M2/228‑251 (MQ=255)
               cctttatcacgtacatggagtattttaagtaggcgttttccgtaaattcatggagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgcgccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAtt                                                                                                                                                                                                                         <  2:376454‑M2/28‑1 (MQ=255)
               ccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAttt                                                                                                                                                                                                                        <  2:358280‑M2/29‑1 (MQ=255)
                             atggagtagtttaagtaggagttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgcttcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattacctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgcgccgacggcgataacgacatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGtt                                                                                                                                                                                                          >  2:290687‑M2/209‑251 (MQ=255)
                                          agtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggatatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATcc                                                                                                                                                                                             <  1:61673‑M2/56‑1 (MQ=255)
                                             aggcgttttccgtaaattcatgtagctcaaggcgctctgcaatatcgctcattaattctgattcctcaacttattcgccagcctgcaaaggggcaatattgccggggatactgcattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGa                                                                                                                                                                                          <  2:318193‑M2/59‑1 (MQ=255)
                                                      ccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAAt                                                                                                                                                                                 <  1:435568‑M2/68‑1 (MQ=255)
                                                          aaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGtt                                                                                                                                                                              <  1:416605‑M2/71‑1 (MQ=255)
                                                          aaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                                                             >  1:293757‑M2/180‑251 (MQ=255)
                                                              tcaggtagctcaaggcgctctgccatatcgatcattaagtctgattccgcagcttattcgccagccttcaaaggggcaatattgccggagatactacatgatctgccgcgtcgagtcacaggagatcaacgtgaatagcgcctccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCTGGGATTGAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCttt                                                                                                                                                                          <  1:224130‑M2/75‑1 (MQ=255)
                                                               catgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTcc                                                                                                                                                                        <  1:164757‑M2/77‑1 (MQ=255)
                                                               catgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttagtcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                                                                         <  2:3371‑M2/76‑1 (MQ=255)
                                                                  gaagcgcaaggcgttctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCttt                                                                                                                                                                     <  2:480940‑M2/80‑1 (MQ=255)
                                                                         aaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgcgtcacatgagatcaacgtgaatagctccgccggacgtgctccgacggcgctaaaggcatcaGCATAACCGTAGTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCNTGACCGGGCAATGTTCTTTCCTTTAAACCTc                                                                                                                                                              >  2:257990‑M2/165‑251 (MQ=255)
                                                                          aggcgcactgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCAGCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCa                                                                                                                                                             <  2:421794‑M2/88‑1 (MQ=255)
                                                                           ggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCa                                                                                                                                                             >  1:84054‑M2/163‑250 (MQ=255)
                                                                              gctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtct                                                                                                                                                         <  1:485478‑M2/92‑1 (MQ=255)
                                                                                        atcgctcatgaattctgattccccaacttattcgccagccttcaaaggggcaatattgccggggatactacattatcggccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggtgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCAGCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCACTCTCCGCCATGTc                                                                                                                                               <  2:38579‑M2/102‑1 (MQ=255)
                                                                                                  aattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTctct                                                                                                                                     >  1:385669‑M2/140‑251 (MQ=255)
                                                                                                           tcctccacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgcgtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAAcc                                                                                                                              >  1:151715‑M2/131‑249 (MQ=255)
                                                                                                           tcctcaccttattcgccagccttcgaaggggcaatagtgccggagataatacattatctgcgtcgtcgggtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAAGAACCGGAGTCGGCAGCCAGGGATTTAATCCGGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTGCCTTTAAACCTCAATCTCCGCCAGGTCTCCTATCTCTTGCATCCa                                                                                                                             <  2:289187‑M2/120‑1 (MQ=255)
                                                                                                                aacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcg                                                                                                                       >  1:264379‑M2/126‑251 (MQ=255)
                                                                                                                 acttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgatacgacggcgaaaaaggcaacaGAATACCCGAATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCACTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACATGTTgcgg                                                                                                                      >  2:479895‑M2/125‑251 (MQ=255)
                                                                                                                     agtcgccagcgctcaaagtggcaatattgccggagatactccattaactgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGAt                                                                                                                  <  1:267404‑M2/131‑1 (MQ=255)
                                                                                                                          ccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTcc                                                                                                             <  1:178769‑M2/136‑1 (MQ=255)
                                                                                                                                 tcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCg                                                                                                        <  1:266474‑M2/141‑1 (MQ=255)
                                                                                                                                     aggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttctt                                                                                                  <  1:525495‑M2/147‑1 (MQ=255)
                                                                                                                                                 gccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGCGCGTTTCTTCGCCAGCAGCat                                                                                      >  1:213310‑M2/93‑251 (MQ=255)
                                                                                                                                                 gccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCat                                                                                      >  1:219635‑M2/93‑251 (MQ=255)
                                                                                                                                                   cggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCatat                                                                                    >  1:420341‑M2/91‑251 (MQ=255)
                                                                                                                                                     gagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATcc                                                                                  <  1:9449‑M2/163‑1 (MQ=255)
                                                                                                                                                          actacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTGGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatcat                                                                             >  1:277411‑M2/84‑251 (MQ=255)
                                                                                                                                                             acattcccggccgcgtcgattcaccggagagcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGTGATTAAATCCGGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTGCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCTGCGAGCTTCCGAGCGTTTCTGCGCCAGCAGCATATCCATCATcgc                                                                          <  2:110295‑M2/171‑1 (MQ=255)
                                                                                                                                                             acattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgc                                                                          <  1:18485‑M2/171‑1 (MQ=255)
                                                                                                                                                             acattatctgccgcgtcgagccccgggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                                                                            <  2:109041‑M2/169‑1 (MQ=255)
                                                                                                                                                              cattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                                                                         >  2:401935‑M2/80‑251 (MQ=255)
                                                                                                                                                                      gccgcgtccagtcacaggagatcaacgtgaagagcgcctccggacgtgctccgacggcgataaaggcttcaGAATATCCTGATTCTGCAGCCAGGTATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCGTGCAACCATTTGCTGCGATCTTCCGAGCGTTTCTTCGCCTTCAGCATAGCCATCATCGCGTCAGTAc                                                                  <  2:230069‑M2/179‑1 (MQ=255)
                                                                                                                                                                          cgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTg                                                              >  1:288582‑M2/68‑250 (MQ=255)
                                                                                                                                                                          cgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTg                                                              >  1:233460‑M2/68‑250 (MQ=255)
                                                                                                                                                                           gtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGAt                                                            >  1:251288‑M2/67‑251 (MQ=255)
                                                                                                                                                                                  cacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtc                                                     >  1:225789‑M2/60‑251 (MQ=255)
                                                                                                                                                                                  cacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtc                                                     <  1:401935‑M2/192‑1 (MQ=255)
                                                                                                                                                                                       gagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTcatcat                                                >  2:399304‑M2/55‑251 (MQ=255)
                                                                                                                                                                                         gatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGACTAACCGTATTCTGCAACCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACAGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCg                                              >  1:57383‑M2/53‑251 (MQ=255)
                                                                                                                                                                                                   aatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACt                                    <  1:399304‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                        cgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCAcc                               <  1:5166‑M2/214‑1 (MQ=255)
                                                                                                                                                                                                        cgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCAcc                               <  1:339482‑M2/214‑1 (MQ=255)
                                                                                                                                                                                                               ggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGc                        >  1:222270‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                   gtgccccggcggcgataaaggcatcaGAATAACCGTATTCTGCAGCCATGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCAGCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAgt                    <  2:213310‑M2/225‑1 (MQ=255)
                                                                                                                                                                                                                         cgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCgg              <  1:409395‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaggcatcaGCATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCCTGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGAGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTg    >  1:78004‑M2/11‑251 (MQ=255)
                                                                                                                                                                                                                                     aggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGCCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGt  >  1:502054‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
CCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGT  >  NC_000913/3173189‑3173668

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