Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 921,243 G→T 6.3% M299I (ATG→ATT macB → macrolide ABC transporter peremase/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913921,2430GT 6.3% 270.9 / 2.8 78M299I (ATG→ATTmacBmacrolide ABC transporter peremase/ATPase
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base T (3/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CGCAATCCTCCCGCCATTGAAAAAGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTA  >  NC_000913/921013‑921478
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
cgcAATCCTCCCGCCATTGAAAAAGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCCGGCCGCGGAAAAAAATCGGTCCTTTCTTACCATGCTGGGGGTTTTTTTCGGGTATTGGGTGGTGGTTTTCCTTTGTGGGGGTGGGTGCGGCGCCCCAAAAAAGGtggtggg                                                                                                                                                                                                                                      >  2:826545/1‑238 (MQ=255)
 gcAATCCCGCCGCCATTGAAAAAGGGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCg                                                                                                                                                                                                                                     <  1:6911/239‑1 (MQ=255)
       cccccGGCATTGAAAAACGGAATGCGACGTGCGGGATGGAACCTGTGGTCAAGACGGTGTCTGGGTGGCGGAAGGTTGTCAGCGGTTTTATCGAGGCGCTGACGATGGCATGGCAGGCGCTGGCAGCGAATAAAATGCGTTCTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGATTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATAtt                                                                                                                                                                                                                               <  2:882882/237‑1 (MQ=255)
         tccGCCAGTGAAAAAGTGGATGTTACTGGCGGGACGGAACCTTTTGGCAACACGGTGTCTGGCTGGCGGGAGTTTGTCAGCGGTTTTAACGAGGCGCTGGCGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCg                                                                                                                                                                                                                             <  2:860385/238‑1 (MQ=255)
         cccGCCATTGAAAAAGTGAATGTTACTGGCGGGACGGAACCTTTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCg                                                                                                                                                                                                                             <  1:761740/239‑1 (MQ=255)
                ttGAAAAAGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTAtt                                                                                                                                                                                                                      <  1:947902/239‑1 (MQ=255)
                    aaaaGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGt                                                                                                                                                                                                                   <  1:1034382/238‑1 (MQ=255)
                    aaaaGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGt                                                                                                                                                                                                                   <  1:826545/238‑1 (MQ=255)
                    aaaaGGGAATGTTACTGGCGGGACGGAACCTGTTGGCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGt                                                                                                                                                                                                                   <  2:1251383/238‑1 (MQ=255)
                     aaaGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGAGGATATTCGTTCTATTGGTAc                                                                                                                                                                                                                 <  1:202186/239‑1 (MQ=255)
                     aaaGTGAATGTGACCGGCGGGGCGGAGCCTGTGTGCAACTCGGTGTGTTGCTGGCGGCAGTTGGTCAGCGGTTTTAACGAGGCGCAGACGATGGCATGGCGGGATCTGGCAGCGAATAAAATGCGTATTTTACTGACCATGGTGGGGATTATTATCGGTATTGCGTCGGTGGTGTCCATTGTCGTGGTGGGTGACGCCGCCAAGCAAATGGTGCTGGCGGATATACGTTCTATTGGTAc                                                                                                                                                                                                                 <  2:1011840/239‑1 (MQ=255)
                           aaTGTTTCTGGGGGGACGGAAACTGTTTTTAACACGGTTGCTGGCTGGGGGCAATTTTGCAGCGGGTTTAACGAGGCGCTGACGATGGGATGGCGGGCGCTGGGAGCGAATAAAATGCGTTCTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTTGGTTCCATTTGCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATa                                                                                                                                                                                                            <  2:1325707/238‑1 (MQ=255)
                           aaTGTTACTGGCGGGACGGACCCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATa                                                                                                                                                                                                            <  2:1220171/238‑1 (MQ=255)
                             tGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGTGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGGATTGCGTCGGGGGTTTCCATTGTCGTGGTGGGGGACGCCGCCAAACAAATGGTGCTGGCGGAGATTCGTTCTATTGGGACGGATACTa                                                                                                                                                                                                         >  1:251777/1‑239 (MQ=255)
                                     gcggggcgggACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGGCAGCGGGTTTTACGAGGCGCTGACGCTGGCATGGCGGGCGGTGTGCGGGAAAAAAAAGCGTACCTTTCTTGCCATTCTGGGGGTTATTTATCGTAATTGGTGCGTGGTTTTTCTTTCCGGTGGTGGGGGACGCCGCCGAACAAATGTTGCTGTTGGGTTTTGGTTTTATTGGTCCgcaacttatttgt                                                                                                                                                                                                    >  2:476724/1‑225 (MQ=255)
                                                   gttgtCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGGGCTGGCGGGTATTCGTTCTATTGGTACGGATAATATTTATGTCTTTCCCGGGAAAGa                                                                                                                                                                                   >  1:897989/1‑239 (MQ=255)
                                                        cAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGAATTGAGTCAGTGGTTTCCCTTGTCGTGGTGGGGGACGCCGCCAAACACATGGGGCGGGCGGATATTCGTTTTATTGGTAAGAATACTATTGATGGCGATCCCGGGAAAGATTTTg                                                                                                                                                                              >  1:425554/1‑239 (MQ=255)
                                                              ggTGTCTGGCTGGCGGCAGATTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGGCGGGCACGCATAAAATAAGTTCTTTTCTGACCGTGGTGGGGGTTATTATCGGTCATGAGTAGGGGGTTTTCATTGACGCGGTGGGTGATGCAGCCAAACAAATGGCAGTTGTGGATATTAGCTCTTTTTGTTCGGATACTATAGATTCGCATCCCGGGAAAGATTTTTGctgcg                                                                                                                                                                        >  2:1134300/1‑235 (MQ=255)
                                                                ggtctgtggggcggcAGTTTGTCAGCGGTTTTAACTAGGCGCTGACGATGGCTTGACGGGCGCTGGCAGCGAATAAAATGCGTACTTGTCTGACCATTCTGGGGATTATTAGCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTAGCCCGGGAAAGATTTTTGCGATGa                                                                                                                                                                       <  2:919863/229‑1 (MQ=255)
                                                                 gtCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGGACTTTACTAACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCCAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATTCCGGGAAAGATTTTGGCgcgcacg                                                                                                                                                                     >  1:1039623/1‑233 (MQ=255)
                                                                  tCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACg                                                                                                                                                                     <  1:466658/238‑1 (MQ=255)
                                                                   ctggctggCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGGACTTTACTGACCATGCCGGGGATTATTATCGGTACTGCGCTGGTGGGTTCCATTGTCGTGGTGGGGGAGCGCGCCAAACATATTGTGCCTGCGGATTTTTCGTCTATTGGCACGGATACCATTTTTGTTTTTTCCGGGATAGAtttttgggactgggtt                                                                                                                                                                   >  2:549433/1‑227 (MQ=255)
                                                                    tggctggcGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATc                                                                                                                                                                  <  1:1111047/239‑1 (MQ=255)
                                                                    tggctggcGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATc                                                                                                                                                                  <  1:645935/239‑1 (MQ=255)
                                                                          gcggcAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGCATAAAATGCGTACTTTACTGACCATGCTGGGGGTTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGGGGGTGACGCCGCCAAACAAAAGGGGGTGGCGGAGTTTTGTTCTGTTGGTACGAACACTTTTTAATTCTATCCCGGGGAAGATTTTTGCGCTGTAGATCCggcat                                                                                                                                                            >  2:127420/1‑235 (MQ=255)
                                                                                      cAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGgcg                                                                                                                                                 <  1:425903/238‑1 (MQ=255)
                                                                                                               atggcatggcGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCa                                                                                                                       >  1:541320/1‑238 (MQ=255)
                                                                                                                tggcatggcGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCAt                                                                                                                      <  1:1383240/238‑1 (MQ=255)
                                                                                                                      ggcgggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATTGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCaaaaa                                                                                                               >  1:585544/1‑239 (MQ=255)
                                                                                                                      ggcgggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTTGCGATGACGATCCGCAATATCAACCGGGGCTGAAGTACGACGACTTAATCGCCATCCaaaa                                                                                                                >  1:504293/1‑238 (MQ=255)
                                                                                                                      ggcgggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCaaaaa                                                                                                               <  2:1202492/239‑1 (MQ=255)
                                                                                                                      ggcgggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGGGGGTGACGCCGCCCAACAAATGGTGCTGGCGGATATTCGTTCTATTTGTACGAAAACTATTTATGTCTATACCGGGAAAAATTTTTGGGGTGACAGTCCCCCATTACACCAGGCGCCGAAAGACGACCACCTTAATCGcctcacaaaa                                                                                                                >  2:1130185/1‑231 (MQ=255)
                                                                                                                                 ggCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTggg                                                                                                     <  1:225338/239‑1 (MQ=255)
                                                                                                                                        gAATAGAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGc                                                                                                 <  2:1040598/237‑1 (MQ=255)
                                                                                                                                        gAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGc                                                                                                 <  1:1130185/237‑1 (MQ=255)
                                                                                                                                            aaaaGGCGTTCTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCTATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATAGATGTCTAACCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCAAAAAACAACCGTGGGTCGCCTCAGc                                                                                           <  1:1419131/239‑1 (MQ=255)
                                                                                                                                                    tACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCt                                                                                    >  1:107254/1‑238 (MQ=255)
                                                                                                                                                       tttaaggACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTGTCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATGCTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCg                                                                                <  2:107254/233‑1 (MQ=255)
                                                                                                                                                              aCCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGGTGTTTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATAGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCggtgcgca                                                                         >  2:1352166/1‑232 (MQ=255)
                                                                                                                                                              aCCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCa                                                                         >  1:151485/1‑239 (MQ=255)
                                                                                                                                                                 aTGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCaaa                                                                       >  1:623678/1‑238 (MQ=255)
                                                                                                                                                                 aTGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCaaa                                                                       >  1:439184/1‑238 (MQ=255)
                                                                                                                                                                 aTGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCaaa                                                                       >  1:777322/1‑238 (MQ=255)
                                                                                                                                                                   gCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCaaaa                                                                      >  2:1056722/1‑237 (MQ=255)
                                                                                                                                                                   gCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAc                                                                     >  2:385740/1‑238 (MQ=255)
                                                                                                                                                                        ggATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTAtt                                                                                                                                                                                                       <  1:233284/104‑1 (MQ=255)
                                                                                                                                                                        ggATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTAtt                                                                                                                                                                                                       >  2:233284/1‑104 (MQ=255)
                                                                                                                                                                          attattATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGGCGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  2:581422/238‑1 (MQ=255)
                                                                                                                                                                          attattATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  1:686530/238‑1 (MQ=255)
                                                                                                                                                                          attattATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  1:476724/238‑1 (MQ=255)
                                                                                                                                                                          attattATCGGTATTGCGTCGGTGATTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCAGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  1:549433/238‑1 (MQ=255)
                                                                                                                                                                          attattATCGGTATTGAGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  1:1126982/238‑1 (MQ=255)
                                                                                                                                                                          attattATCGGGATTGGGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcc                                                              <  2:1389957/238‑1 (MQ=255)
                                                                                                                                                                                     ggTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGGTATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTTGAGGACGACGACCTAATCGCCCTCCCAAAACAACCGTGGGTCGCCCCCCCCACCCCTGCCGTCCCGCAAAACCTGCCACTGCtttttt                                                     >  2:1059329/1‑235 (MQ=255)
                                                                                                                                                                                     ggTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTtata                                                     >  1:242711/1‑238 (MQ=255)
                                                                                                                                                                                     ggTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATaa                                                    <  1:1277805/239‑1 (MQ=255)
                                                                                                                                                                                          tGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAatg                                               >  2:469474/1‑239 (MQ=255)
                                                                                                                                                                                          tGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTAGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAatg                                               >  1:1020982/1‑239 (MQ=255)
                                                                                                                                                                                          tGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCCATTGGTACGAATCCTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGGACGACGACTTAATTGCCATCCAAAAACAACCGTGGGTCGCCTCATCCCCCCCTGCCGTCTCGCAAACACTGCGCCCGCGTTATCACCatg                                               >  2:849635/1‑239 (MQ=255)
                                                                                                                                                                                              tCGGTGGTTTCCATTGTCGTGGCGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAACAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCACAACCTGCGCCTGCGTTATAACAatgttga                                           >  1:1135161/1‑239 (MQ=255)
                                                                                                                                                                                                 gtggtTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATTGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGAGGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATcagc                                                                                                                                                      <  2:104571/130‑1 (MQ=255)
                                                                                                                                                                                                 gtggtTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATTGTGCTGGCGGATATTCGTTCTATTGGTACGAATAATATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATTCGCAATATcagc                                                                                                                                                      >  1:104571/1‑130 (MQ=255)
                                                                                                                                                                                                 gtggtTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcctg                                                            >  1:362083/1‑219 (MQ=255)
                                                                                                                                                                                                 gtggtTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAAcctgcgcctg                                                            <  2:362083/219‑1 (MQ=255)
                                                                                                                                                                                                 gtggtTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAatgttgatgt                                        >  2:1105840/1‑239 (MQ=255)
                                                                                                                                                                                                           aTTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATTTTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAg                               <  1:1084174/239‑1 (MQ=255)
                                                                                                                                                                                                           aTTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAg                               <  1:780642/239‑1 (MQ=255)
                                                                                                                                                                                                                  tggtggGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAAt                        >  2:1299777/1‑239 (MQ=255)
                                                                                                                                                                                                                    gtgggtgACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAAc                                                   <  1:1161289/210‑1 (MQ=255)
                                                                                                                                                                                                                    gtgggtgACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAAc                                                   >  2:1161289/1‑210 (MQ=255)
                                                                                                                                                                                                                    gtgggtgACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATg                       >  2:1071231/1‑238 (MQ=255)
                                                                                                                                                                                                                         tgACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGa                 >  1:1095710/1‑239 (MQ=255)
                                                                                                                                                                                                                                gccAAACAAATTGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGa          <  2:585544/239‑1 (MQ=255)
                                                                                                                                                                                                                                     acaaaTGGTGCTGGCGGATATTCGTTCTATTGGTACGaa                                                                                                                                                                                                              <  1:652334/39‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaTGGTGCTGGCGGATATTCGTTCTATTAGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTAtttt   <  1:1134300/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGGTGCTGGCGGATCTTCGTTCTATTGGCACGAATAATATTGATGCCTATCCCGGGATAGATCTTGGCGATGAAGATCCGCAATAGCAGCAGGCGCTTGAGTACGACGACTTAAACGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGATGTCAGTGCCAATGGCGTGAGCGGCGATTATTTTa  <  2:657803/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTa  <  2:439184/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTa  <  1:127420/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
CGCAATCCTCCCGCCATTGAAAAAGTGAATGTTACTGGCGGGACGGAACCTGTTGTCAACACGGTGTCTGGCTGGCGGCAGTTTGTCAGCGGTTTTAACGAGGCGCTGACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACAAATGGTGCTGGCGGATATTCGTTCTATTGGTACGAATACTATTGATGTCTATCCCGGGAAAGATTTTGGCGATGACGATCCGCAATATCAGCAGGCGCTGAAGTACGACGACTTAATCGCCATCCAAAAACAACCGTGGGTCGCCTCAGCCACACCTGCCGTCTCGCAAAACCTGCGCCTGCGTTATAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTA  >  NC_000913/921013‑921478

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