Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 921,692 T→G 19.9% V449G (GTG→GGG)  macB → macrolide ABC transporter peremase/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913921,6920TG19.9% 94.6 / 2.1 52V449G (GTG→GGG) macBmacrolide ABC transporter peremase/ATPase
Reads supporting (aligned to +/- strand):  ref base T (15/24);  new base G (9/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTACTTTACA  >  NC_000913/921473‑921911
                                                                                                                                                                                                                             |                                                                                                                                                                                                                             
tttttAATTTTTATGGCCTGACCTTCAGTGAAGGAAACCCCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCgt                                                                                                                                                                                                                        <  2:632091/226‑1 (MQ=255)
    taaggtttttGGCATGCCCTTAAGGGGAGGAAACACCTTTAATCAGGGTCAGCTGAACGGTCGT‑‑GCCCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCt                                                                                                                                                                                                                    <  2:33800/218‑1 (MQ=255)
       gggcTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCCCGAGAGTAATACTCGCCGCCAGGTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGc                                                                                                                                                                                                                <  1:584916/224‑1 (MQ=255)
          cTATGGCATGGCCTTCAGTGAAGGAACCACCTTTAATCAGGAAGAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTAcc                                                                                                                                                                                                            <  1:13564/228‑1 (MQ=255)
           tATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGAAAAACAGTCGATGTTTGTTAGCTGTAAAGTGCTGCGGGTCTGGCTAACt                                                                                                                                                                                                           >  1:54571/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGTGTCGGGCTACCTTACAg                                                                                                                                                                                                      >  1:544518/1‑227 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCGTTGGTGTGGCGGCAGAAAAACAGTCGAGGTTTTGTAGCAGTAAAAGGATGCGGGACTGGCGAACTTACaaaac                                                                                                                                                                                                   >  2:45525/1‑224 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGATGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:106295/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGGGTCTGGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:85211/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGAAAAACAGTCGATGTTTGTTAGCAGTAAAGGGCTGCTGGTCTGGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:133574/1‑228 (MQ=255)
                              aaGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGAACGATGTCCggg                                                                                                                                                                                         >  1:444021/1‑227 (MQ=255)
                                      acCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGGGTCTGGCTACCTTAAAGCACGATGTCCGGGCGGGTTATg                                                                                                                                                                                >  1:144453/1‑228 (MQ=255)
                                                 ggacccgcTTAACGGTCGT‑‑GCCCAGGTCGTGGTTCCCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCg                                                                                                                                                                       <  1:189903/220‑1 (MQ=255)
                                                     cagcTCAAGGGTTGG‑‑GGGCCGGTCGGGGGTCTCGACAGTAACACGTGCCGCCAGCTTTTCCCCCAGAAAGCAGATGGGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 <  2:106295/228‑1 (MQ=255)
                                                       gcTGAACGGTCGT‑‑GCGCAGGGCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTaa                                                                                                                                                               <  1:327777/228‑1 (MQ=255)
                                                           aaCGGCCGT‑‑GCCCAGGTCGTGGTTCGCGGCAGTAATACTCGCCGCCGCTTTTTCCCCCAAAAAGCGAATGTGGTTGGGGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACt                                                                                                                                                             <  2:444021/226‑1 (MQ=255)
                                                                  gT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGCAAAACAGTCGCTGTTTGGTAGACGTAAAGGGCTGCGTGGCTGGACACCTTAAAGCACGCTGTCCGGGCGGGTTATGGGCCAGCCGGGGCTTAACCCCATAACtt                                                                                                                                                    >  2:414433/1‑227 (MQ=255)
                                                                      gcgcAGGTCGTGGTTCTCGACAGTAAGACTCGCCGCCTGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTc                                                                                                                                                  <  1:179099/228‑1 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 >  2:223446/1‑207 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 <  1:223446/207‑1 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGAGCTGCGTGTCTGGCTACCTTACAGCACGATGTCAGGGCGAGTTATGGGCCAGTCGTGGCTTAACTACATTACTGTCAggggg                                                                                                                                            >  1:303402/1‑226 (MQ=255)
                                                                                    tctcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGGGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGACAGGGTGAaagaag                                                                                                                                    >  1:288778/1‑227 (MQ=255)
                                                                                      tcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGGTGTTTGGTAGCAGTAAAGGGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCTGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGGGACAGAAGGAt                                                                                                                                 >  1:321847/1‑228 (MQ=255)
                                                                                                  cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGCGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGTGTCTGGGTACCTTACAGCAAGATGTCCGGGCGAGTTAAGGGCCAGTGGGGGCTTAAATCCATTACTGTCAGGGGGAAAGAAGGATTTTGCAGAGcca                                                                                                                     >  2:474802/1‑227 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTTGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGc                                                                                                                                                                             <  2:266931/170‑1 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGc                                                                                                                                                                             >  1:266931/1‑170 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTGGTGGCTTAACTCCCTTATTGTCAGGGGGAAAGAAGGGTTTGCAGGCGCCgggg                                                                                                                  >  2:418863/1‑228 (MQ=255)
                                                                                                       cgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTACTGTCAGGGGGGAAGAAGGATTTGACAGCGCCgcggcg                                                                                                                >  1:214260/1‑228 (MQ=255)
                                                                                                                  cccccATAAAGGCTGTGGTGTTGGGTAAGTTATTTTGGGCGGAAATATTCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACt                                                                                                     <  1:474802/228‑1 (MQ=255)
                                                                                                                     ctaaaagtcAGGTGGGGTTGGGGGGGTGGTTCTGTGCGGCAATATGCGCGCCAGGGTCAGGGGTGTGGTGGAAGAAAAACAGGCGATGTTTGGTGGCGGTAAAGTTCTGCGTGTCTGGCTACCTTACAGCACGAGTTCCGGGCGGGTTATGGGGCAGTCGTGGCTTAACGCCATTACGGTCAGGGTGAAAGAAGGATTGGACAGCGCCGAGGCGGAACAGCAACTCa                                                                                                   <  2:187437/219‑1 (MQ=255)
                                                                                                                               gttgtgGTTGGGGGGTTGATTCTGGTCGGAAATAGGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACtt                                                                                        <  2:73413/226‑1 (MQ=255)
                                                                                                                                aTGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATTTTTGGTAGCAGTAAAGAGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGGTAGTTATGGGCCAGTCGGGGCTTAACTCCATTACTGGCAGGGTTGAAGAATGATTTGACAGGGCCCAGGCGGACAAGCAACTCACGGGTTTCCttt                                                                                       >  1:498777/1‑228 (MQ=255)
                                                                                                                                aTGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGGTTGGTAGCATTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGCCAGGGTGAAGTAAGGATTTAACAGCGCCGAGGCGGAAACGCAACTCACGCGTTTGCttt                                                                                       >  1:301724/1‑228 (MQ=255)
                                                                                                                                          cGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCAAGGCGGACCAGCACCTCAAGCGTTTCATTTTACTGCGCCa                                                                             >  1:11944/1‑228 (MQ=255)
                                                                                                                                                       gTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAACAAAGTCGATGTTTGGCAGCAGTAAAGGGCGGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTCCTGTCAGGGGGACAGACGGGTTTTGACGGGCCGCGGCGGACAAGGCACTCACGCGTTTACTTTCACTCCGACACGTAAAAGAGGAt                                                                >  2:431435/1‑228 (MQ=255)
                                                                                                                                                       gTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTTCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGa                                                                 <  1:400922/227‑1 (MQ=255)
                                                                                                                                                       gTCGGCAATACGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGGTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGa                                                                 <  1:45525/227‑1 (MQ=255)
                                                                                                                                                            cAATATGCCCGGCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAACACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCtt                                                           <  1:47824/228‑1 (MQ=255)
                                                                                                                                                                      gCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTACTGTCCGGGGGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGCACGAAGGATTTCTTTACCTGGAAc                                                 >  1:159242/1‑228 (MQ=255)
                                                                                                                                                                      gCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGCAGCAGTAAAGGGCTGCGTGTCTGGCTAACTTACAGCACAATGTCCGGGGGAGGTATGGGACAGTCGTGGCTTAAATCCATTACTGTACGGGGGACAGAAGGCTTTTACAGCGGCGATGCGGCACAGGCACCAGCGCGATAAATTTGACTGGGCGACGCAAAGCAGGAGTTCTTATCCTgggaa                                                 >  2:451489/1‑225 (MQ=255)
                                                                                                                                                                             tttttGGTGGGGCGGAAGAAAAACATTTGATGGTTGGTAGCAGTAAGGTTCTGCGTGTCTGGCTTCCTTACAGCACGATGTCCGGGCGGGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACg                                          <  2:288778/225‑1 (MQ=255)
                                                                                                                                                                                       tgCGGAAGAAAAACCTTCGATTTTTGGTAGGAGTAAAGTGCTGGGTGTCTGGCTTCCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCtt                                   <  2:303402/224‑1 (MQ=255)
                                                                                                                                                                                       ggCGGAAGAAAAACAGTCGATGTTTTGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGGTATGGGCCAGCAGGGGCTTAACTCCATTACTGTCCGGGCGAAAGAAGGCTTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCCCTGGCCCACGGAAAGAAGGATTTCTTTACCGTGAAAATTGCCGGCGGCTTGaa                                >  2:44362/1‑228 (MQ=255)
                                                                                                                                                                                                 aaaccgttcaTGTTTGGTAGGAGGAAAGGGGTGGTGTGTGGGCTGACTTACAGCATGATGTCCGGGCGAGTTATGGGCCAGGCGTGGCTTAACTCCATTACTGGCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  2:85211/219‑1 (MQ=255)
                                                                                                                                                                                                 aaaCAGTCGAGTTTTGGTAGCAGTAAAGTGCTGCGTTTCTGGCTACCTTACAGCACGATGGCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  2:214260/228‑1 (MQ=255)
                                                                                                                                                                                                  aaCAGTCGATGTTTGGTAGCAGGAAAGTGCTGCGTGTCTGGCTTCCTTCCAGCACGATGGCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  1:500128/227‑1 (MQ=255)
                                                                                                                                                                                                     aGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAg                   <  1:364651/227‑1 (MQ=255)
                                                                                                                                                                                                            gTTTGGTTGCAGGAAAGGGCTGCGTGTCTGGCTTCCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGGCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACg           <  2:321847/228‑1 (MQ=255)
                                                                                                                                                                                                            gTTTGGTAGCAGTAAAGTGCTGCGTGGCGTGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACGCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACg           <  1:431435/228‑1 (MQ=255)
                                                                                                                                                                                                              tgtgTAGGAGTGTAGTGTTGGGTGTCTGGCTCCCTTGACGCACGAGTGCCGGGGGAGTTTTGGGCCCTTCGTGGCTGTACTCCATTGATGGCCGGGTGAAAGAAGGATTTGGCAGGGGCGAGGGCGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGACAGAAGGAATTCATTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTa         <  2:54571/225‑1 (MQ=255)
                                                                                                                                                                                                                     cAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGGCTTGAAAACTGTTGAAAAAACCACACGTACTTTaaa  >  2:306593/1‑226 (MQ=255)
                                                                                                                                                                                                                     cAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGAGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGCCTTGCAAACTGTGGAAAAGACCACACGGACTTTAc   >  1:306534/1‑227 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                             
ATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTACTTTACA  >  NC_000913/921473‑921911

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: