Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,879,936 T→G 22.3% T62P (ACC→CCC)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,879,9360TG22.3% 127.8 / 6.9 68T62P (ACC→CCC) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (18/33);  new base G (15/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACC  >  NC_000913/3879699‑3880165
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
tatatttCCCCGCGCGGTGGATACCCAGTCTAGGTTTGTGGGACAGGGCTTTGCTCACCGAACCACCACCGCCGTGCAGACCGCGGGATCTTTTATGGGTGTTAGCGCCAAATTTGCCGCCTGCGTGCAGAACGGTCATGATCACTCCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTg                                                                                                                                                                                                                                      <  2:320698/235‑1 (MQ=255)
           tCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGGGCAGAACGGTCATGATCACTTCCGCCGCCGATCCGCCCCCTTCCGGGAGAATACCGGTCGGACTCCCGCGCCCGACATCCTGTACAGAGACAGAGGTATAGAGGAGGATGGGGACGAAAAta                                                                                                                                                                                                                             >  2:309357/1‑236 (MQ=255)
           tCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGGGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCTCCTTCCGGGTGAATACCGGTCGGAATCCCGCGCCCGCCATCCTGTAACAGGCAAGGGTTATCGGCGTGGAGGGGGAGGAAAAtttt                                                                                                                                                                                                                           >  2:698734/1‑238 (MQ=255)
                             tCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGCGGGTAATTTTTTTCCACGGACCCCCGAg                                                                                                                                                                                                         >  1:229065/1‑239 (MQ=255)
                             tCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGGGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGCCAGAGTTATCGGCGTGCATGTGGACGAGAATTTTTTTACAGGGACCcgcgc                                                                                                                                                                                                          >  1:1107007/1‑237 (MQ=255)
                                       tGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGGGAAAGGGGACGGTAATTTCTTTACAGTGACCCGCGAGCGCTtcgtc                                                                                                                                                                                                >  2:125018/1‑239 (MQ=255)
                                          gtcggggCGTTTCTGACCCAACCCCCAACGCGTGGCAGACCGCCGACAATTTTTAAGGAGTTATCACAAAATTTACCCACTGGGTGAAGAAGGGTCAGGATCATTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGa                                                                                                                                                                                              <  2:870482/234‑1 (MQ=255)
                                          gACAGGGCATTTACTACCGAAACCCAAACGCCGTGCAGACCGCCGGACACTTTATAGGGTTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGa                                                                                                                                                                                              <  1:703526/238‑1 (MQ=255)
                                                                       cGCCGTGCAGCCCGCGGGACACTTTATAGGAGTTACCGTCAAATTTCCCGCCTGCGTGCAGAACGGTAATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCAtg                                                                                                                                                                 <  1:962813/239‑1 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGCCAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCCTCGGCGTTAGAGGTATCTACAACCCCGAATaaccagtgggg                                                                                                                                                            >  2:213861/1‑229 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGGCGATAATTTCTTTTCAGTGACCCGCCAGCGCTTCGTCGATAGAGGTATTAACCACCCCCAAAACCaattgggg                                                                                                                                                            >  2:168495/1‑232 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGACGGTAATTTCTTTaaa                                                                                                                                                                                                                     >  1:181600/1‑180 (MQ=255)
                                                                                         actcTTTAAAGGGTTTACCGTAAAGTTGCCCGCCGGCGTGCAGAACGGTCATGATCACTCCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGcc                                                                                                                                                <  1:125018/235‑1 (MQ=255)
                                                                                             ttttaagggAGTTTCCGTAAATTTTCCCGCCTGCGTGCAGAACGGTCATGATCCCTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCa                                                                                                                                            <  1:413180/232‑1 (MQ=255)
                                                                                             tttATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAAAGGTGACGATAATTTCTTTTCAGTGACCCGCGAGCCCTTCGTCGGTGGCGGTTTATACCACCCCCAATACCACGTGGGGGAGACCGGTCCCGGCa                                                                                                                                            >  2:761159/1‑238 (MQ=255)
                                                                                                       gTTTTCGTCATTTTTCCCGCTGGCGTGCGGAGCGGTCAGTGTCCCTCCCGCCGCCTGTACGCCCTGTTACGGGTGAACGCCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCg                                                                                                                                   <  2:691464/238‑1 (MQ=255)
                                                                                                        ttATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCTCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCCTCGGAGCTAGCGTTAACTAACACCTCGAAAAACATGTGGTGCCGACCGGGGGCGGCAACCGTGACGca                                                                                                                                 >  2:793951/1‑238 (MQ=255)
                                                                                                        ttATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGCCGGTAATTTCTTTACAGTGCCCCCCGGGCGCTTCGCGGTTAGCGGTATTTTACACCCCGAATACCCTTTGGTGAACACCGGAGCCGTCACCCggtttcc                                                                                                                                  >  2:64312/1‑232 (MQ=255)
                                                                                                         taccGTCAAATTTCCCGGCTGTGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  1:1012611/236‑1 (MQ=255)
                                                                                                          ttCGTCAAATTTTCCCGCGTGCTGGCAGAACGGTCAGGATCCCTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATGCCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  2:46096/238‑1 (MQ=255)
                                                                                                          aTCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATTTGGTGCAGACCGGGGCCGTCATCCGGGTCGcccg                                                                                                                               >  1:92112/1‑237 (MQ=255)
                                                                                                          aTCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCGCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTATAGTGACAGAGTTATCGGCGTGCATGGTGACGATAATTTCTTTACa                                                                                                                                                                                                                     >  1:472115/1‑153 (MQ=255)
                                                                                                                         cGCCTGCGGGCAGAACGGGCAAGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCggg                                                                                                                 <  1:213861/239‑1 (MQ=255)
                                                                                                                                           tcctGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCggg                                                                                                                 <  1:391191/218‑1 (MQ=255)
                                                                                                                                                aTCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGCTAGCGTTATCTACCACCCCGAATACCAGTGGGTGCCGACCGGGGCCGCCCTCCGTGTCGCCGCAATACATAACCCGGCGCTTACGCACCGCACCCAGc                                                                                            >  2:339936/1‑237 (MQ=255)
                                                                                                                                                    ctccccccgccgATACCCCCTATTCCGGGTGAATACCTGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTTACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTc                                                                                      <  2:184758/233‑1 (MQ=255)
                                                                                                                                                     ttCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGGCGCTAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGAATGCGTTATCTTCCACCTCGAAAACCATGTGGGGCAGACCGGGGCCGTCACCCCTCTCCCCCGTTTACATACCCGGGCGCTTCCGCCCCCCAACCAAcccatccg                                                                                     >  2:568964/1‑234 (MQ=255)
                                                                                                                                                           cgccgATTCGCCCTCTTCCGGGGTAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACtt                                                                               <  2:92112/239‑1 (MQ=255)
                                                                                                                                                             ccgATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGAt                                                                                                                                                                                                                                 <  1:244374/91‑1 (MQ=255)
                                                                                                                                                             ccgATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGACGAt                                                                                                                                                                                                                                 >  2:244374/1‑91 (MQ=255)
                                                                                                                                                                 tACGCCCGCTGCCGGGTGAATACGGGTCGGAATGCCGCGCACGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTCACGATAATTTCTTTACAGTGGCCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATGCCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATa                                                                          <  1:64696/238‑1 (MQ=255)
                                                                                                                                                                     cccTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTTCCGTCATCCGTGTCGCCGATATACATACCCGCGCGCTTTCGCACCCCATCCAGCCCTTTTCGGACTTTCATACTtga                                                                     >  1:479352/1‑239 (MQ=255)
                                                                                                                                                                      cctcctCCGGGGTAATCCCGGTCGGAATGCCGCGACCGTCATCCTTGACAGGAGCAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGGCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTggag                                                                    <  2:698661/239‑1 (MQ=255)
                                                                                                                                                                          ttCCGGGTTAATCCCCTTCGGAAGTGCGGGCCCGGACTCCTTTACAGGGACAGCGTTATCGGCGTTAATGGTGACGATAATCTCTGTACAGTGCCCCTCGAGGCCTTCGTCGATAGCGTTAGCTACCGCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGCCTTTGATACTGGAGGAGt                                                                <  2:836142/239‑1 (MQ=255)
                                                                                                                                                                          gccaggttGAATTCCGGTCGGAATGCCCCCGCCGTCAACCTGTACAGAGACAGAGTTATCGGCTGGAATGGTGACGATAATTTCTTTACAGGGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTggagga                                                                  <  2:339940/230‑1 (MQ=255)
                                                                                                                                                                             cGGGTTAATCCCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGGTATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCAt                                                             <  1:568964/239‑1 (MQ=255)
                                                                                                                                                                                gTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGat                                                                                                                              <  1:248296/171‑1 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGCGACAGCGTTATCGGCGTGAACGGGGACGATAATTTCTTTACAAGGACCCCCCAGCGCTTCGTCGATATCGTTATCTACCACCTCGTATACCATGTGGTGCCGAGCGGGGACGCCTTCCGTGTCGCCGATTTACATACCCGGTCGCCTACTGCCCGCCCCCACCCCTTTCACGAACTTTAAACTCGAGGAGtacaaaaaaat                                                      >  2:370449/1‑229 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGGATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCCTCCATAACGTTATCTCCCACCTCGAGAAACATGTGGTGCAGACCGGCGCCGTCATCCGGTGCGCCGACATCCAAATCCTGGCGCTTCGCCCCCGCACCCATCCCTTTCGCAACTTTGTAACTGGAGGGGCCACAAGaaat                                                      >  2:266624/1‑237 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAAtt                                                      <  1:168495/239‑1 (MQ=255)
                                                                                                                                                                                     tACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTc                                                     <  1:579498/239‑1 (MQ=255)
                                                                                                                                                                                         gggtggAATGCCGCGCCCGTCATCCTTGTAAGGGAACGAGTTATCGGCGTGAATGGGTACGATAATTTCCTTCCAGTGCCCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACa                                                 <  2:714792/235‑1 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATa                                                                                                                                                                      <  1:168085/121‑1 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCCCATCCAGCCCTTTCCGGACTTTGATCCTGGAGGAGTCCAAAGAATTCGACAt                                                >  1:412792/1‑239 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCCGGCGCTTACGCACCGCATCCCGCCCTTTCAGGACTTTGATATTGGCGGAGCCATAAGaa                                                        >  1:339940/1‑231 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATa                                                                                                                                                                      >  2:168085/1‑121 (MQ=255)
                                                                                                                                                                                              gAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGGATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCCTTTGGGGCAGACCGGTGCCCCCCTCCCTGTCGCCGCGATACATACCCCGGCGCTTACGCACCCCCTCCCCCCCCTTCCCGACCTTTATACTTGGGGGGTCAAAAGAAtttaaacacaaa                                            >  2:1011147/1‑229 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                               <  2:213618/59‑1 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                               >  1:213618/1‑59 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGGGGAGTCCTAAGAATTCGACATCaaag                                           >  1:714792/1‑237 (MQ=255)
                                                                                                                                                                                                     gcgcCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGGGGCCTAAGAATTCGACATCAACGttttc                                      >  1:908296/1‑236 (MQ=255)
                                                                                                                                                                                                      ccccctttacccTGTACAGAGACAGAGTTTTCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTATGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCg                                    <  2:412792/229‑1 (MQ=255)
                                                                                                                                                                                                             cATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCAtt                              <  1:20228/238‑1 (MQ=255)
                                                                                                                                                                                                                 cTGTACAGAGACAGAGTTATCGGTGTGAATGGTGACGATAATTTCTTTACCGTGACCCGCGGGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATAc                         <  1:309357/239‑1 (MQ=255)
                                                                                                                                                                                                                  ctttACAGAGACAGAGTTATCGGTGTTAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCCACCCCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATAc                         <  1:698734/236‑1 (MQ=255)
                                                                                                                                                                                                                   tGTACAGAGACAGAGTTACCGGCGTTAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACt                        <  1:726933/239‑1 (MQ=255)
                                                                                                                                                                                                                     tACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                               <  1:269645/38‑1 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                >  1:500456/1‑128 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  2:500456/128‑1 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTAAGCACCGCATCCAGCCCTTTCAGGACCTTGATACTGGAGGAGGCATAAGAATTCGACATCAACGTTTTTTGCTCATTTTTACTTTGGTTa                 >  1:411566/1‑239 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCCTTATACAAACCCGGGCGCTTACGACCCGCCTCCACCCCTTTCAGGGCTTTGTTACTGGGGGAGGCATAAGAATTCGACCTCAACCTTTCTCCCTCATTTTTCTCTGGGTTa                 >  2:521872/1‑239 (MQ=255)
                                                                                                                                                                                                                             acagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTTCGCACCGCATCCAGCCCTTTCAGGACTTTGATACTTGAGGAGGCATAAGAATTCGACATCAACGATTCCCGCCCCTTTATACTTGGGTTaa                >  1:535529/1‑238 (MQ=255)
                                                                                                                                                                                                                                gagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCCGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTTTCGCTCATTTATACTTGGGTTAATCCg            >  1:999033/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGGCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCg            <  1:498538/237‑1 (MQ=255)
                                                                                                                                                                                                                                     gtCGGCGTTAATGGTTACGATAATTTCTTTACAGTGGCCCGCGAGCGGTTCGTCGATAGCGTTATCTACCGCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTAtt       <  1:691132/238‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGTTACTGGAGGAGGCATAAGAATTCGACATCACCGTTTCCCGCCCATTTTTACTTGGGTTAATCCGTTAtttttc   >  1:586584/1‑235 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGGGCCCCGCGAGCGCTTCGGCGATAGCGTTATCTACCCCCTCGAATACCATGTGGTGCAGACCGGGGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTAcc  <  2:999033/238‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGGGACCCGCGAGCGCTTCGTCGATCGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGTTTCCGTCATCCGTGTCCCCGATATACATCCCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACATTGTAACTGGGGGTTTATAAAGAATTCGACAGCCCCCTTTTTCTCCCATTTTACCTTGTGGTAATACGTttttttcc   >  2:691132/1‑235 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
GAATTTTACCCTCGCGCTGGATAACCAGCTCCAGTTTTTGCGACAGGGCGTTTACTACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACC  >  NC_000913/3879699‑3880165

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