Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,879,062 A→T 15.7% M353K (ATG→AAG)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,879,0620AT15.7% 169.0 / 5.6 74M353K (ATG→AAG) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (2/10);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878840‑3879293
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
gAGCCCCCTCCCACCAGGTGTGATCCGGAAGGGTCGGGATCAGGTCCCGGGCTGTCTTCCAGTTTCCCCGGCGGCCCCTTTAGGCCGCGGGAACCTTTCCGGCGGGTATTTTCCCGCGCGCGACCCCCCTCTTCGGGGGCCCGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCTTCGTTTGGGTTTGCCAGCAGATTTTCTGCCAGAAGTTCGTTCATCTGCTGTTCAAccg                                                                                                                                                                                                                          <  2:171776/238‑1 (MQ=255)
         tCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCCGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTCATCCCCGCGCCGGCCCGCGCCCCGGGCCCGGGCAGTATAGCTGATTTTGGCCAAACCTCGCTTCGCGTCGGTTGCGGTTTTTTGCAGGCACTATGCCACCGCTGCGTTCTTCTTCTGCTCTAccgccggcatca                                                                                                                                                                                                                 >  2:146631/1‑232 (MQ=255)
           cgccgggaaCAGTCGGGAACGCGCCGGACGGCGTCCCGGGCATCCTGCAGGTTGCCCGGGCCGCCCGCCAAGGGCGAGCGCCACTTGCCGGGGGGTCTTTTCCGGCGCGGGCCGCGCCGCTTCACGGGCACGGGCACCATCGATATTTTTGCCACCCCCGTTTTTCGCGCGGGTTGGGTTTTCAGGCAGTTATTCCGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCt                                                                                                                                                                                                              <  2:1367909/231‑1 (MQ=255)
             ccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTTACCCGCGCGGCGCGCCCCCTCACGGGGACGGGGGGGACCCGGAAATTTTCCAACCACCATTTTCGCGGCGGTTGGGTGTTTTAGCCGGGATTCTGCCAGCCGTTCGTTTATCTTCTGTTCATccgcccccttctcccc                                                                                                                                                                                                             >  2:596705/1‑227 (MQ=255)
                  tACAGTTCGGAAAGCGCCGGATCGCGTTCCGGGCATTCTGCCAGTTGGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCagaa                                                                                                                                                                                                         <  1:308699/237‑1 (MQ=255)
                   gcAGTTCGGGAACCGCGGGTGCGCGTTCCGGGCATTCGGCCAGTTTCCCCGGAGGGCCCGCTAAGTCGGGCGCTCCTTTACGGGGGGTCATTTCCCGCGCGCGACGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCGTGGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAgaagaa                                                                                                                                                                                                      <  2:11523/238‑1 (MQ=255)
                     gcttggaaacgcggcgGTTCGCGTCCCGGGCATTCTCCGATTTGCCCCGCCAGGCCCCTAAAGTCGGGCGCCCCTGTGGGGCGGGTCATTCCCCGCGCGCGCCGCGCCGCTTCACGGGCACGGGCGGCATCGATAATTTTCCCAACCACGTTTTTCGTGTCGGTTGGGTTTTCCAGCGGGTTTTCTCCAAGCAGTTCGTTCATCTGTTGTTAAACCGCCGATTTCACCTCAGAAGAAAc                                                                                                                                                                                                    <  2:636172/229‑1 (MQ=255)
                           gAACGCGCCGGACCGCGTCCCTGGCAGTCGGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCCCTTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTCCCAACCACGATTTCCGCGTCGGTTGGGTTTTCCAGCAGGTTTTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGtt                                                                                                                                                                                               <  2:396778/238‑1 (MQ=255)
                           gAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTTTCGCGCGCGCCGCGCCGCCTCCCGGCCCCGGGCAGCATCGGCAATTTTTCCACACCCGACTTTTTCGTCGGGGTGGGTTTCCCGCCGGGATTTTGTCTCCAGCTACTTTTTTTTCTGTTCCAACGCCGCCTTCTCCTCCCTAGAAAACAGttt                                                                                                                                                                                              >  2:1336818/1‑239 (MQ=255)
                                        cgcgTTCCGGGCAGTCTGCGAGTTCGCCCGGCGGGCCGGCTAAGCGGAGCGCCCCTTTCCGGCGGGTCTTTTCACGGGCGCGCCGCGCCGCTTCCCGGGCGCGGGGAGCATCGATAATTTTGCCAACCACGATTTTCGGGTGGGTTGGTTTTCCAGCAGTTATTCTGCCAGCATTTCTTTCTTCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTgg                                                                                                                                                                                 <  2:625087/238‑1 (MQ=255)
                                          cgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCACTTTTCGCGCGCGCCGCCCCGCTTCACCGGGACGGGGGGCAACGATAATTTTGCCAACCACCACTTTTGCGGCGGTTGGGGTTTTCCGCAGGGGTTTTTCCCCACGGTCTTTCATCTGCTGTTCAAACCCCCCGTTTTCCTCACAAGAAACCCACTTGTCttttttttgggt                                                                                                                                                                               >  2:939279/1‑230 (MQ=255)
                                            tcccgggCAGTTCGCCAGTTTGCCCGGCCGGCCCGCGAAGCCGAGCGCCCCTTTCGGGGGGGTATTTTCACGCGCGCGCCGCGCGGTTCCCCGGGAGCGGGCAGCATCGATAATTTTGCAAACCACGTTTTCCGGGTGGGTTGGTTTTTCCAGCAGGTATTCTCCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGggagg                                                                                                                                                                             <  2:628268/234‑1 (MQ=255)
                                                                              gCTAGGTCGAGCGCACCTTTACGGCGGGTCTTTTCCCGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGAc                                                                                                                                           <  1:492065/239‑1 (MQ=255)
                                                                                      cgggccctcttttCCGGGGGGTCATTTCCCGCGCGCGCGGCGCCGCTGCCGGGGCACGGGCAGCCTCGATATTTTTGCCACCCCCTTTTTGCGTTTTGGTTGGTTTTTCCGGCAGGTATTCGCCCACCAGTTCGTTCATCTGTTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATc                                                                                                                                    <  2:628023/229‑1 (MQ=255)
                                                                                                 tgagggggggTCATTGCCCGCGCGCGCGGCGCCGCTTCACGGGCAGGGGCAGCATCGATAATTTTCCCCCCCACGTTTTCCGCGTCGGTTGGGTTTCCCAGCAGGTTTTCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCAcg                                                                                                                         <  2:900402/232‑1 (MQ=255)
                                                                                                     ggcggGTCATTCCCCGCGCGCGCCGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTTCCCAACCACGATTTCCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:146091/239‑1 (MQ=255)
                                                                                                      gcggGTCTTTTCACCCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:1326172/238‑1 (MQ=255)
                                                                                                      gcggGTCATGTCAGGCGCCGGCCGCGCGGCTACAGGGGCACGGACAGCGTGAATATTTTTCCAACCCACGTTTTTCGCGTCGTTTGGTTTGTCCAGCAGGTTTCCGCCCAGCAGTTCGTTCATCTGTTTTCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  2:1297212/238‑1 (MQ=255)
                                                                                                                    gcgcgggCCGCCCCGCTTCACGGGCCGGGGCAGAGTCGTATATTTTCCCAACCACTATTTTCGCGCTGTTTGGGTTTCCAAGCAGGTATTCTGCAAGCATTTGTTTAATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgcg                                                                                                      <  2:742935/239‑1 (MQ=255)
                                                                                                                     cgggcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGGTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGCAACCAGTTTGTCTTTGGTCTGGGCGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGACTTCACGCGCATCCTAACCGgcggcg                                                                                                      >  2:872418/4‑238 (MQ=255)
                                                                                                                     cgcgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGATAATTTCGGGTCCGGCACTTTTACGGAAACGACCGCAATCAGGCCTTCACGCGCATAGTCACCGGTTgc                                                                                                       >  2:488126/1‑237 (MQ=255)
                                                                                                                       cgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGGGAATTTCGGGTCCGGCACTTTCACGGAAACGCCCGCAATCTGGCCTTCAGGCGCTTCGTCCCCGCTGGCCCtc                                                                                                   >  2:61391/1‑238 (MQ=255)
                                                                                                                         cgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTg                                                                                                   >  1:322136/1‑237 (MQ=255)
                                                                                                                                        cgggcacgggcaGATTGGATAATTTTCCAACCCACGATTTTCGCGTGGGTTGTGTTTCCCAGCGGGTTTTCGGCCAGCATTTCTTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:336593/239‑1 (MQ=255)
                                                                                                                                         gggcacgggaaGCATCGATAATTTTGCCACCCACGATTTTCGCGTCGGTGGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:1198522/238‑1 (MQ=255)
                                                                                                                                             acgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:422835/1‑238 (MQ=255)
                                                                                                                                             acgggcaGCATCGATAATTTTCCCAACCCCGATTTTCGCGTCGTTTGGGTTTTCCAGCAGTTTTTCGGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  1:440348/239‑1 (MQ=255)
                                                                                                                                             acggccaGCATCTATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTTTTCTGCCAGCATTTCGTTCATCTGCTGTTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  2:1096902/239‑1 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:783619/1‑237 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:68176/1‑237 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:1064874/1‑237 (MQ=255)
                                                                                                                                                  gagtatgaaTAATGTTGCTATCCACGATTTTCGGTGGGGTTGGGTTTTCCATGGGTGATTCCGTGCACAGTTTGTACATCGGATGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttc                                                                         <  1:445884/230‑1 (MQ=255)
                                                                                                                                                   agcaTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttctt                                                                       >  2:1214009/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTTGCCTTCTTTg                                                                     >  2:245373/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:893336/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:508849/1‑239 (MQ=255)
                                                                                                                                                      aTCGATAATTTTGCCAACCGCGATTTTCGCGTCGTTTGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCGGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCttt                                                                      <  1:1359505/237‑1 (MQ=255)
                                                                                                                                                      aTCGATAATTTGGCCACCCACGATTTGCGCGTCGGTTGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGt                                                                    <  1:532057/239‑1 (MQ=255)
                                                                                                                                                       tagATAATTTTCCCACCCACGATTTTCGCGCCGGTTGTGTTTCCCAGCGGGTATTCTGCCAGCATTTCGTTCATCGGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttct                                                                        <  2:946703/232‑1 (MQ=255)
                                                                                                                                                          aTAATTTTGCCAACCACGCTTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTTACCGCAGTTCGTTCATCTGCTGTTCAACCGCCGCTTTCACCTCAGAAGGAACCAGTTTGTCTTTGGTCTGGGAGCGGCATTTCGGGTCCGGCACTTTCACGGCCTCGACCGCAATCCGGCATTCACGCGCATCGTCACCGGGGGCGTTGACTTTCGCTTTTTTGCTGCAGCATTCTTTGTCCAt                                                                >  2:747158/1‑239 (MQ=255)
                                                                                                                                                              ttttGCCAACCACTATTTGCGGGTGGTTTGGGTTTTCCGGAAGGTATTCTGCCAGCAGTTCGTGCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:740951/239‑1 (MQ=255)
                                                                                                                                                               tttgtccAACGACTGTTTTCGCGTCGGTTGTGTTTGCCAGAAGGTATTCTCCCAGCATTTCTTTCAGCTGTTGTTAAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:488126/234‑1 (MQ=255)
                                                                                                                                                                ttGCCAACCACATTTTTCGGGCCGTTGGGTTTTTCAGCCTGGTTTTCTGCCACGGTTTGTTTTACCTCTGTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGc                                                          <  2:488704/239‑1 (MQ=255)
                                                                                                                                                                 tGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                                                              >  1:946703/1‑234 (MQ=255)
                                                                                                                                                                    ccAACCCCGTTTTCCGCGTGGTTTGGGTTTTCCAGCAGGTATGCGGCGAGCGGTTCGTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:410838/239‑1 (MQ=255)
                                                                                                                                                                    ccAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCCGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTATTTGTCCATGTAGGCGt                                                        >  1:1210082/1‑238 (MQ=255)
                                                                                                                                                                    ccAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATCCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:928743/239‑1 (MQ=255)
                                                                                                                                                                    caaaCCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCTTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:122586/237‑1 (MQ=255)
                                                                                                                                                                     cccccccgtttttTGCGGCGGGTGGGTTTTTCCAGCAGGTTTTCCGCCGGCTGTTCCTTCCTTTGCTTTCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:1228690/230‑1 (MQ=255)
                                                                                                                                                                     cAACCACTATTTTCGCGTGGGTTGGGTTTTCCAGCAGGTTTTCGGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:960852/238‑1 (MQ=255)
                                                                                                                                                                      acccACGTTTTGCCCGGCGGTTGGGTTTTCCAGCAGTTTTCCACCAAGATGTTCGTTAATCGGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCa                                                     <  1:331492/237‑1 (MQ=255)
                                                                                                                                                                                   gcggcgGTGGGTTTTCCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCGCCTGTTCAGCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTGTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGGCCATGTAGGCGTTCAGGGTACGGGTCa                                         <  1:1116401/238‑1 (MQ=255)
                                                                                                                                                                                   cgttggtttGTGTTTTCCAGCAGTTATCCCGCCAGGAGTCCGTCCATCTGGTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCAt                                        <  1:376224/232‑1 (MQ=255)
                                                                                                                                                                                           tGGGTTTTCCAGCAGGTTTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCa                                 <  1:596705/238‑1 (MQ=255)
                                                                                                                                                                                             ggTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACgg                              >  1:1313841/1‑239 (MQ=255)
                                                                                                                                                                                                ttttCCAGCGGTTATTCTGCCAGCTGTTCGTTTATCTGCTGTACAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:897731/238‑1 (MQ=255)
                                                                                                                                                                                                tttcCCAGCAGTTTTTCGGCCGGTAGTTCGTACATTGGCTGATCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:1064205/238‑1 (MQ=255)
                                                                                                                                                                                                 ttcccAGCGGGTTTTCTCCCCGCAGTTCTTTCTTCTCTTGTTCAACCGCCGATTTCCCCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:188007/236‑1 (MQ=255)
                                                                                                                                                                                                  ttcccagaggttATCCTGCCGGCAGTGCGTTCTTCTGCGGTTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:1357132/228‑1 (MQ=255)
                                                                                                                                                                                                          gggTTTTCTGCCAGCCGTTCTTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGcc                     <  1:1244337/236‑1 (MQ=255)
                                                                                                                                                                                                           ggTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  >  2:842048/1‑239 (MQ=255)
                                                                                                                                                                                                            gTATTCTGCCAGCAGTTCTTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  <  1:1057828/238‑1 (MQ=255)
                                                                                                                                                                                                              aTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:971927/1‑239 (MQ=255)
                                                                                                                                                                                                              aTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:127132/1‑239 (MQ=255)
                                                                                                                                                                                                                  tGCCAGAATTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAc           <  1:1185375/239‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcc          >  1:1303228/1‑238 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgcc       >  2:426679/1‑239 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGGGAGTAccgc        >  2:791272/1‑238 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGCAGCCTTCTTTGTCCATGTAGGCGTTCATGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcagc        >  2:1357594/1‑238 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:784689/1‑239 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:77958/1‑239 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:206676/1‑239 (MQ=255)
                                                                                                                                                                                                                        cagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:536294/1‑238 (MQ=255)
                                                                                                                                                                                                                           ttcgttcATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCGGCCAGGTGAGTACCGCCGTccc  >  1:962044/1‑237 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878840‑3879293

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