Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,409,558 G→T 6.5% T16T (ACG→ACT) ‡ dbpA → ATP‑dependent RNA helicase, specific for 23S rRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,409,5580GT 6.5% 249.9 / 5.4 85T16T (ACG→ACT) ‡dbpAATP‑dependent RNA helicase, specific for 23S rRNA
Reads supporting (aligned to +/- strand):  ref base G (41/38);  new base T (3/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCCACCACCCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGA  >  NC_000913/1409409‑1409707
                                                                                                                                                     |                                                                                                                                                     
cccACCACCCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACg                                                                                                                                                       >  2:1358227/1‑150 (MQ=255)
cccACCACCCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACg                                                                                                                                                       >  1:1379566/1‑150 (MQ=255)
   accaccCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAAc                                                                                                                                                    >  2:1426770/1‑150 (MQ=255)
      accCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCtt                                                                                                                                                 >  1:209658/1‑150 (MQ=255)
      accCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCtt                                                                                                                                                 <  2:417796/150‑1 (MQ=255)
       cccGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTa                                                                                                                                                <  2:42910/150‑1 (MQ=255)
          gATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATg                                                                                                                                             >  1:806183/1‑150 (MQ=255)
             gaagaaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCTCGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGt                                                                                                                                          <  1:1010921/150‑1 (MQ=255)
             gaagaaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGt                                                                                                                                          <  1:479762/150‑1 (MQ=255)
             gaagaaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGt                                                                                                                                          <  2:736/150‑1 (MQ=255)
               agaaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTg                                                                                                                                        <  2:453586/150‑1 (MQ=255)
                gaaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTgg                                                                                                                                       >  2:256415/1‑150 (MQ=255)
                 aaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTggg                                                                                                                                      <  2:1703736/150‑1 (MQ=255)
                 aaCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTggg                                                                                                                                      <  1:1188826/150‑1 (MQ=255)
                 aaCCGATTAAATTGTCGCATCGTGAGCGAAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTggg                                                                                                                                      <  2:1432563/150‑1 (MQ=255)
                  aCCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGt                                                                                                                                     >  2:28860/1‑150 (MQ=255)
                   ccGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCATTAAGGTTCGTGAGTTgggcg                                                                                                                                    <  1:204936/150‑3 (MQ=255)
                    cGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTATGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTGTTGCCTCCCGCACAACTCACGAACCTTAAAGAGTTGGGTTa                                                                                                                                   >  2:861843/1‑150 (MQ=255)
                     gATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCACCCAACTCATTAAGGTTCGTGAGTTgggcggg                                                                                                                                  <  1:936989/150‑5 (MQ=255)
                     gATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATCTTTTGCCTCCCGCCCAACTCACGAACCTTAATTAGTTGGCTTAt                                                                                                                                  >  1:252356/1‑150 (MQ=255)
                        tAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGCCCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTa                                                                                                                               >  2:1222295/1‑150 (MQ=255)
                          aaTTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAAc                                                                                                                             >  1:1595108/1‑150 (MQ=255)
                           aTTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAAcc                                                                                                                            <  1:266235/150‑1 (MQ=255)
                            ttGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCTTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCa                                                                                                                           <  2:1242091/150‑1 (MQ=255)
                              gTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATg                                                                                                                         >  1:471515/1‑150 (MQ=255)
                              gTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGAGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATg                                                                                                                         >  2:515873/1‑150 (MQ=255)
                               tCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGa                                                                                                                        >  2:148094/1‑150 (MQ=255)
                                cGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTAGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGAc                                                                                                                       <  1:300765/150‑1 (MQ=255)
                                           cGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATTAGTTGGGTTATTTAACCATTACGCCGGTGCAgg                                                                                                            <  1:1002540/150‑1 (MQ=255)
                                             tAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGgcc                                                                                                          >  1:916271/1‑150 (MQ=255)
                                                  cTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGTGTTGGGTTATTTAACCATGACGCCGGTGCAGGCcgccgc                                                                                                     >  1:944967/1‑150 (MQ=255)
                                                     cGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCt                                                                                                  >  1:1110553/1‑150 (MQ=255)
                                                     cGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCt                                                                                                  <  1:927362/150‑1 (MQ=255)
                                                     cGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCt                                                                                                  <  1:1243068/150‑1 (MQ=255)
                                                       aGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTc                                                                                                >  2:1184470/1‑150 (MQ=255)
                                                        ggttACGCTCCTGCTAACTTTTCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGATCCGTACTGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTGcc                                                                                               <  1:916722/147‑1 (MQ=255)
                                                          aTACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCgg                                                                                             >  1:211851/1‑150 (MQ=255)
                                                             cgcgCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGa                                                                                          >  1:1139310/1‑150 (MQ=255)
                                                             cgcgCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGa                                                                                          >  2:214348/1‑150 (MQ=255)
                                                                 ccTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTATTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCAAGACGCACGTGCAGGCCGCCGCGCTTCCGGCGATCCt                                                                                      >  1:83038/1‑150 (MQ=255)
                                                                     cTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGcc                                                                                  <  1:736404/150‑1 (MQ=255)
                                                                     cTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGcc                                                                                  <  2:1164687/150‑1 (MQ=255)
                                                                     cTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATCTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGcc                                                                                  <  2:98901/150‑1 (MQ=255)
                                                                     cTAACTTTGCGTCGATGACCACGAGAATAGATTCTGACCGCGTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGcc                                                                                  <  1:81120/150‑1 (MQ=255)
                                                                       aaCTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCgg                                                                                >  2:1296323/1‑150 (MQ=255)
                                                                       aaCTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCgg                                                                                >  1:949321/1‑150 (MQ=255)
                                                                         atttGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTAACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGaa                                                                              >  2:1188101/2‑150 (MQ=255)
                                                                          tttGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGaaa                                                                             <  1:377501/150‑1 (MQ=255)
                                                                              cGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGAt                                                                         >  2:1645948/1‑150 (MQ=255)
                                                                                 cGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGtt                                                                      >  1:363837/1‑150 (MQ=255)
                                                                                  gATGACCACGAGAATAGATTGAGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTc                                                                     <  1:515873/150‑1 (MQ=255)
                                                                                   aTGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTcg                                                                    >  1:470643/1‑150 (MQ=255)
                                                                                       ccACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTTATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTg                                                                >  1:1506981/1‑150 (MQ=255)
                                                                                       ccACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTg                                                                >  2:954473/1‑150 (MQ=255)
                                                                                          cGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAg                                                             >  2:972513/1‑150 (MQ=255)
                                                                                           gagaAGAGATTGTGACCGCTTTTTCTACCGTGACTGTTTTGCCTCAGGGCAAACTCACTAACCTTAATGTGTTGTGTTTTTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAgg                                                            <  2:140649/150‑1 (MQ=255)
                                                                                            agaATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGc                                                           <  2:211851/150‑1 (MQ=255)
                                                                                             gaATAGATTGGGACCGCTTTTTCTACCGTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCg                                                          <  2:83038/150‑1 (MQ=255)
                                                                                               aTAGATTGTGACCGCTTTTTCTTCCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGGGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGaa                                                        <  1:1322464/150‑1 (MQ=255)
                                                                                               aTAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGaa                                                        <  2:806183/150‑1 (MQ=255)
                                                                                               aTAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGaa                                                        >  1:409852/1‑150 (MQ=255)
                                                                                                   aTTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAAcc                                                    <  2:1568819/150‑1 (MQ=255)
                                                                                                      gtgACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGc                                                 >  1:314976/1‑150 (MQ=255)
                                                                                                         aCCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGc                                              >  1:73917/1‑150 (MQ=255)
                                                                                                         aCCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCAAGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGc                                              >  2:385798/1‑150 (MQ=255)
                                                                                                          ccGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCg                                             >  1:1417081/1‑150 (MQ=255)
                                                                                                                    tACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAAcggc                                   >  1:801124/1‑150 (MQ=255)
                                                                                                                     aCCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAAcggcg                                  >  1:317799/1‑150 (MQ=255)
                                                                                                                           aaTGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCtttt                            <  2:916271/150‑1 (MQ=255)
                                                                                                                            aTGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTg                           <  2:1295258/150‑1 (MQ=255)
                                                                                                                               ttttGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGcc                        <  1:1222295/150‑1 (MQ=255)
                                                                                                                                tttGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTGGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCt                       <  1:256415/150‑1 (MQ=255)
                                                                                                                                   gCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCgg                    >  2:67994/1‑150 (MQ=255)
                                                                                                                                     cTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCt                  <  2:1379566/150‑1 (MQ=255)
                                                                                                                                      tCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTGCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCtt                 <  1:57858/150‑1 (MQ=255)
                                                                                                                                       cccGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTg                >  1:573297/1‑150 (MQ=255)
                                                                                                                                        taacccAACTCATTAAGGTTCGTGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGt               >  1:1261079/4‑150 (MQ=255)
                                                                                                                                        ccGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGt               >  1:1675945/1‑150 (MQ=255)
                                                                                                                                         cGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGtt              <  1:981961/150‑1 (MQ=255)
                                                                                                                                         cGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGtt              <  2:1139310/150‑1 (MQ=255)
                                                                                                                                         aacccAACTCATTAAGGTTCGTGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGtt              >  1:5736/3‑150 (MQ=255)
                                                                                                                                          acccAACTCATTAAGGTTCGTGAGTTGGGTTATTTAACCATGATGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTa             >  2:1020835/2‑150 (MQ=255)
                                                                                                                                            ccAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACa           >  1:475747/1‑150 (MQ=255)
                                                                                                                                             cAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAg          <  1:1296323/150‑1 (MQ=255)
                                                                                                                                             cAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAg          <  2:949321/150‑1 (MQ=255)
                                                                                                                                                cTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCaa       >  1:1631707/1‑150 (MQ=255)
                                                                                                                                                  cACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAAt     <  1:1358227/150‑1 (MQ=255)
                                                                                                                                                     gAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGa  >  2:1642308/1‑150 (MQ=255)
                                                                                                                                                     gAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGa  >  2:1498188/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCACCACCCGATGAAGAACCGATTAAATTGTCGCATCGTGAGCGTAGATCTGCGAGGATACGCGCCTGCTAACTTTGCGTCGATGACCACGAGAATAGATTGTGACCGCTTTTTCTACCCTGAATGTTTTGCCTCCCGCCCAACTCACGAACCTTAATGAGTTGGGTTATTTAACCATGACGCCGGTGCAGGCCGCCGCGCTTCCGGCGATCCTTGCCGGAAAAGATGTTCGCGTGCAGGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGA  >  NC_000913/1409409‑1409707

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