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‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTACTTTACA  >  NC_000913/921473‑921911
                                                                                                                                                                                                                             |                                                                                                                                                                                                                             
tttttAATTTTTATGGCCTGACCTTCAGTGAAGGAAACCCCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCgt                                                                                                                                                                                                                        <  2:632091/226‑1 (MQ=255)
    taaggtttttGGCATGCCCTTAAGGGGAGGAAACACCTTTAATCAGGGTCAGCTGAACGGTCGT‑‑GCCCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCt                                                                                                                                                                                                                    <  2:33800/218‑1 (MQ=255)
       gggcTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCCCGAGAGTAATACTCGCCGCCAGGTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGc                                                                                                                                                                                                                <  1:584916/224‑1 (MQ=255)
          cTATGGCATGGCCTTCAGTGAAGGAACCACCTTTAATCAGGAAGAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTAcc                                                                                                                                                                                                            <  1:13564/228‑1 (MQ=255)
           tATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGAAAAACAGTCGATGTTTGTTAGCTGTAAAGTGCTGCGGGTCT‑GGCTAACt                                                                                                                                                                                                           >  1:54571/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGTGTCG‑GGCTACCTTACAg                                                                                                                                                                                                      >  1:544518/1‑227 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCGTTGGTGTGGCGGCAGAAAAACAGTCGAGGTTTTGTAGCAGTAAAAGGATGCGGGACT‑GGCGAACTTACaaaac                                                                                                                                                                                                   >  2:45525/1‑224 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGATGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:106295/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGGGTCT‑GGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:85211/1‑228 (MQ=255)
                   gACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGAAAAACAGTCGATGTTTGTTAGCAGTAAAGGGCTGCTGGTCT‑GGCTACCTTACAGCAc                                                                                                                                                                                                   >  1:133574/1‑228 (MQ=255)
                              aaGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGAACGATGTCCggg                                                                                                                                                                                         >  1:444021/1‑227 (MQ=255)
                                      acCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGGGTCT‑GGCTACCTTAAAGCACGATGTCCGGGCGGGTTATg                                                                                                                                                                                >  1:144453/1‑228 (MQ=255)
                                                 ggacccgcTTAACGGTCGT‑‑GCCCAGGTCGTGGTTCCCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCg                                                                                                                                                                       <  1:189903/220‑1 (MQ=255)
                                                     cagcTCAAGGGTTGG‑‑GGGCCGGTCGGGGGTCTCGACAGTAACACGTGCCGCCAGCTTTTCCCCCAGAAAGCAGATGGGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 <  2:106295/228‑1 (MQ=255)
                                                       gcTGAACGGTCGT‑‑GCGCAGGGCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTaa                                                                                                                                                               <  1:327777/228‑1 (MQ=255)
                                                           aaCGGCCGT‑‑GCCCAGGTCGTGGTTCGCGGCAGTAATACTCGCCGCCGCTTTTTCCCCCAAAAAGCGAATGTGGTTGGGGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACt                                                                                                                                                             <  2:444021/226‑1 (MQ=255)
                                                                  gT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGCAAAACAGTCGCTGTTTGGTAGACGTAAAGGGCTGCGTGGCT‑GGACACCTTAAAGCACGCTGTCCGGGCGGGTTATGGGCCAGCCGGGGCTTAACCCCATAACtt                                                                                                                                                    >  2:414433/1‑227 (MQ=255)
                                                                      gcgcAGGTCGTGGTTCTCGACAGTAAGACTCGCCGCCTGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTc                                                                                                                                                  <  1:179099/228‑1 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 >  2:223446/1‑207 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCtt                                                                                                                                                                 <  1:223446/207‑1 (MQ=255)
                                                                            gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGAGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCAGGGCGAGTTATGGGCCAGTCGTGGCTTAACTACATTACTGTCAggggg                                                                                                                                            >  1:303402/1‑226 (MQ=255)
                                                                                    tctcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGGGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGACAGGGTGAaaga‑ag                                                                                                                                    >  1:288778/1‑227 (MQ=255)
                                                                                      tcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGGTGTTTGGTAGCAGTAAAGGGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCTGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGGGACAGA‑AGGAt                                                                                                                                 >  1:321847/1‑228 (MQ=255)
                                                                                                  cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGCGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGGGCTGCGTGTCT‑GGGTACCTTACAGCAAGATGTCCGGGCGAGTTAAGGGCCAGTGGGGGCTTAAATCCATTACTGTCAGGGGGAAAGA‑AGGATTTTGCAGAGcca                                                                                                                     >  2:474802/1‑227 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTTGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGc                                                                                                                                                                             <  2:266931/170‑1 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGc                                                                                                                                                                             >  1:266931/1‑170 (MQ=255)
                                                                                                     gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTGGTGGCTTAACTCCCTTATTGTCAGGGGGAAAGA‑AGGGTTTGCAGGCGCCgggg                                                                                                                  >  2:418863/1‑228 (MQ=255)
                                                                                                       cgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTACTGTCAGGGGGGAAGA‑AGGATTTGACAGCGCCgcggcg                                                                                                                >  1:214260/1‑228 (MQ=255)
                                                                                                                  cccccATAAAGGCTGTGGTGTTGGGTAAGTTATTTTGGGCGGAAATATTCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACt                                                                                                     <  1:474802/228‑1 (MQ=255)
                                                                                                                     ctaaaagtcAGGTGGGGTTGGGGGGGTGGTTCTGTGCGGCAATATGCGCGCCAGGGTCAGGGGTGTGGTGGAAGAAAAACAGGCGATGTTTGGTGGCGGTAAAGTTCTGCGTGTCT‑GGCTACCTTACAGCACGAGTTCCGGGCGGGTTATGGGGCAGTCGTGGCTTAACGCCATTACGGTCAGGGTGAAAGA‑AGGATTGGACAGCGCCGAGGCGGAACAGCAACTCa                                                                                                   <  2:187437/219‑1 (MQ=255)
                                                                                                                               gttgtgGTTGGGGGGTTGATTCTGGTCGGAAATAGGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACtt                                                                                        <  2:73413/226‑1 (MQ=255)
                                                                                                                                aTGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATTTTTGGTAGCAGTAAAGAGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGGTAGTTATGGGCCAGTCGGGGCTTAACTCCATTACTGGCAGGGTTGAAGA‑ATGATTTGACAGGGCCCAGGCGGACAAGCAACTCACGGGTTTCCttt                                                                                       >  1:498777/1‑228 (MQ=255)
                                                                                                                                aTGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGGTTGGTAGCATTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGCCAGGGTGAAGTA‑AGGATTTAACAGCGCCGAGGCGGAAACGCAACTCACGCGTTTGCttt                                                                                       >  1:301724/1‑228 (MQ=255)
                                                                                                                                          cGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCAAGGCGGACCAGCACCTCAAGCGTTTCATTTTACTGCGCCa                                                                             >  1:11944/1‑228 (MQ=255)
                                                                                                                                                       gTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAACAAAGTCGATGTTTGGCAGCAGTAAAGGGCGGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTCCTGTCAGGGGGACAGACGGGTTTTGAC‑GGGCCGCGGCGGACAAGGCACTCACGCGTTTACTTTCACTCCGACACGTAAAAGAGGAt                                                                >  2:431435/1‑228 (MQ=255)
                                                                                                                                                       gTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTTCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGa                                                                 <  1:400922/227‑1 (MQ=255)
                                                                                                                                                       gTCGGCAATACGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGGTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGa                                                                 <  1:45525/227‑1 (MQ=255)
                                                                                                                                                            cAATATGCCCGGCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAACACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCtt                                                           <  1:47824/228‑1 (MQ=255)
                                                                                                                                                                      gCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCCTTACTGTCCGGGGGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGCACGAAGGATTTCTTTACCTGGAAc                                                 >  1:159242/1‑228 (MQ=255)
                                                                                                                                                                      gCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGCAGCAGTAAAGGGCTGCGTGTCT‑GGCTAACTTACAGCACAATGTCCGGGGGAGGTATGGGACAGTCGTGGCTTAAATCCATTACTGTACGGGGGACAGA‑AGGCTTTTACAGCGGCGATGCGGCACAGGCACCAGCGCGATAAATTTGACTGGGCGACGCAAAGCAGGAGTTCTTATCCTgggaa                                                 >  2:451489/1‑225 (MQ=255)
                                                                                                                                                                             tttttGGTGGGGCGGAAGAAAAACATTTGATGGTTGGTAGCAGTAAGGTTCTGCGTGTCT‑GGCTTCCTTACAGCACGATGTCCGGGCGGGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACg                                          <  2:288778/225‑1 (MQ=255)
                                                                                                                                                                                       tgCGGAAGAAAAACCTTCGATTTTTGGTAGGAGTAAAGTGCTGGGTGTCT‑GGCTTCCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCtt                                   <  2:303402/224‑1 (MQ=255)
                                                                                                                                                                                       ggCGGAAGAAAAACAGTCGATGTTTTGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGGTATGGGCCAGCAGGGGCTTAACTCCATTACTGTCCGGGCGAAAGA‑AGGCTTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCCCTGGCCCACGGAAAGAAGGATTTCTTTACCGTGAAAATTGCCGGCGGCTTGaa                                >  2:44362/1‑228 (MQ=255)
                                                                                                                                                                                                 aaaccgttcaTGTTTGGTAGGAGGAAAGGGGTG‑GTGTGTGGGCTGACTTACAGCATGATGTCCGGGCGAGTTATGGGCCAGGCGTGGCTTAACTCCATTACTGGCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  2:85211/219‑1 (MQ=255)
                                                                                                                                                                                                 aaaCAGTCGAGTTTTGGTAGCAGTAAAGTGCTGCGTTTCT‑GGCTACCTTACAGCACGATGGCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  2:214260/228‑1 (MQ=255)
                                                                                                                                                                                                  aaCAGTCGATGTTTGGTAGCAGGAAAGTGCTGCGTGTCT‑GGCTTCCTTCCAGCACGATGGCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGaa                      <  1:500128/227‑1 (MQ=255)
                                                                                                                                                                                                     aGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAg                   <  1:364651/227‑1 (MQ=255)
                                                                                                                                                                                                            gTTTGGTTGCAGGAAAGGGCTGCGTGTCT‑GGCTTCCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGGCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACg           <  2:321847/228‑1 (MQ=255)
                                                                                                                                                                                                            gTTTGGTAGCAGTAAAGTGCTGCGTGGCG‑TGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACGCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACg           <  1:431435/228‑1 (MQ=255)
                                                                                                                                                                                                              tgtgTAGGAGTGTAGTGTTGGGTGTCT‑GGCTCCCTTGACGCACGAGTGCCGGGGGAGTTTTGGGCCCTTCGTGGCTGTACTCCATTGATGGCCGGGTGAAAGA‑AGGATTTGGCAGGGGCGAGGGCGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGACAGAAGGAATTCATTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTa         <  2:54571/225‑1 (MQ=255)
                                                                                                                                                                                                                     cAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGGCTTGAAAACTGTTGAAAAAACCACACGTACTTTaaa  >  2:306593/1‑226 (MQ=255)
                                                                                                                                                                                                                     cAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGAGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGCCTTGCAAACTGTGGAAAAGACCACACGGACTTTAc   >  1:306534/1‑227 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                             
ATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGT‑‑GCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCT‑GGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGA‑AGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAAGGATTTCTTTACCTGGAACATGGACGGCGTCTTGAAAACTGTTGAAAAGACCACACGTACTTTACA  >  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: ‑