Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,583 A→G 5.3% intergenic (‑426/+17) rrsG ← / ← clpB 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,5830AG 5.3% 311.8 / 5.7 94intergenic (‑426/+17)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base A (46/43);  new base G (3/2);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

TCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAC  >  NC_000913/2731438‑2731730
                                                                                                                                                 |                                                                                                                                                   
tCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAgg                                                                                                                                                   <  4:674999/148‑1 (MQ=255)
      tCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGt                                                                                                                                             <  4:588759/148‑1 (MQ=255)
      acAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtt                                                                                                                                            >  2:139545/2‑149 (MQ=255)
       cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGttt                                                                                                                                           <  1:139545/149‑1 (MQ=255)
         gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttttt                                                                                                                                         <  1:975574/149‑2 (MQ=255)
           ggtgtCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGTAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTAt                                                                                                                                        <  4:661422/148‑1 (MQ=255)
           ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTAt                                                                                                                                        <  3:205484/148‑1 (MQ=255)
           ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATc                                                                                                                                       <  1:385938/149‑1 (MQ=255)
              ggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCAt                                                                                                                                     <  4:838570/148‑1 (MQ=255)
              ggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCAt                                                                                                                                     <  3:283926/148‑1 (MQ=255)
                 gcgttATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                 <  2:308603/146‑1 (MQ=255)
                  cATTATAGGGAGTTCTTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                 >  4:1012909/1‑148 (MQ=255)
                  cATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                                <  1:845922/149‑1 (MQ=255)
                          ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttt                                                                                                                                           >  2:1066021/1‑130 (MQ=255)
                          ggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCg                                                                                                                        >  2:886271/1‑149 (MQ=255)
                           gAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCg                                                                                                                        >  3:970156/1‑148 (MQ=255)
                              ttATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACaa                                                                                                                    >  1:26799/1‑149 (MQ=255)
                              ttATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACaa                                                                                                                    >  1:1012034/1‑149 (MQ=255)
                                   ccGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                                <  3:819118/148‑1 (MQ=255)
                                       ccTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCAATTTTCAGGCGAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACCGCGACAATCCGGTCt                                                                                                            <  3:658156/148‑1 (MQ=255)
                                        cTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAGAACATCTTACATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCtt                                                                                                           <  4:727729/148‑1 (MQ=255)
                                         tGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTc                                                                                                          <  3:628606/148‑1 (MQ=255)
                                             aaGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                     >  2:267972/1‑149 (MQ=255)
                                             aaGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                     >  1:594281/1‑149 (MQ=255)
                                               gCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACt                                                                                                   >  1:1063356/1‑149 (MQ=255)
                                                cGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACt                                                                                                   >  4:853634/1‑148 (MQ=255)
                                                       tataAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttttt                                                                                                                                         >  1:1074702/1‑102 (MQ=39)
                                                       tataAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttttt                                                                                                                                         <  2:1074702/103‑2 (MQ=38)
                                                        ataAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGa                                                                                          <  1:512087/149‑1 (MQ=255)
                                                          aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                        >  1:267589/1‑149 (MQ=255)
                                                            aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAt                                                                                                        <  4:577723/131‑1 (MQ=255)
                                                            aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAt                                                                                                        >  3:577723/1‑131 (MQ=255)
                                                            aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                       >  4:219424/1‑148 (MQ=255)
                                                            aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCa                                                                                      <  2:1063356/149‑1 (MQ=255)
                                                            aaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCa                                                                                      <  1:155118/149‑1 (MQ=255)
                                                                 ttATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACttt                                                                                  >  3:487035/1‑148 (MQ=255)
                                                                    tCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTACATTAACTTCCAGGCGAATCACTTTAccc                                                                              <  2:406457/149‑1 (MQ=255)
                                                                      gTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCgg                                                                            >  2:209891/1‑149 (MQ=255)
                                                                           cTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAcc                                                                        >  4:307733/1‑148 (MQ=255)
                                                                            tCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCaa                                                                      >  2:249029/1‑149 (MQ=255)
                                                                              aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCa                                                                                                         <  2:200531/112‑1 (MQ=255)
                                                                              aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCa                                                                                                         >  1:200531/1‑112 (MQ=255)
                                                                                   tCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAcc                                                                <  4:509392/148‑1 (MQ=255)
                                                                                   tCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAcc                                                                <  3:178353/148‑1 (MQ=255)
                                                                                     aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAg                                                              <  3:465868/148‑1 (MQ=255)
                                                                                      ggCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGAc                                                            >  1:129529/1‑149 (MQ=255)
                                                                                      ggCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGAc                                                            >  2:602358/1‑149 (MQ=255)
                                                                                        cAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCAt                                                                                                                                     <  4:117216/74‑1 (MQ=255)
                                                                                        cAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCAt                                                                                                                                     >  3:117216/1‑74 (MQ=255)
                                                                                         aaaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAg                                                          >  3:1045653/1‑148 (MQ=255)
                                                                                           aaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTa                                                        >  2:862770/1‑148 (MQ=255)
                                                                                           aaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTa                                                        <  1:862770/148‑1 (MQ=255)
                                                                                              aTCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTg                                                    <  1:267972/149‑1 (MQ=255)
                                                                                                 ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                 >  2:328904/1‑149 (MQ=255)
                                                                                                  tAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                                <  2:594281/149‑1 (MQ=255)
                                                                                                  tAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                 >  3:1035347/1‑148 (MQ=255)
                                                                                                  tAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                 >  4:65196/1‑148 (MQ=255)
                                                                                                  tAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                 >  3:338004/1‑148 (MQ=255)
                                                                                                      tataGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                     >  4:7359/1‑92 (MQ=255)
                                                                                                      tataGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                     <  3:7359/92‑1 (MQ=255)
                                                                                                        taGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAg                                           <  3:219424/148‑1 (MQ=255)
                                                                                                             ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGt                                      <  4:118797/148‑1 (MQ=255)
                                                                                                              tttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtt                                     >  3:84924/1‑148 (MQ=255)
                                                                                                               ttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                    >  4:236648/1‑148 (MQ=255)
                                                                                                                  cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCg                                 >  3:214069/1‑148 (MQ=255)
                                                                                                                   gTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGa                                <  3:236648/148‑1 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCgg                                        >  1:265176/1‑136 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCgg                                        <  2:265176/136‑1 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                           >  1:279825/1‑149 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                           >  1:145293/1‑149 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                           <  1:249029/149‑1 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                           <  1:849453/149‑1 (MQ=255)
                                                                                                                       aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctg                            <  3:307733/148‑1 (MQ=255)
                                                                                                                           tataGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                        >  4:958765/1‑148 (MQ=255)
                                                                                                                           tataGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                        >  4:271055/1‑148 (MQ=255)
                                                                                                                             taGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGa                     >  2:367321/1‑149 (MQ=255)
                                                                                                                              aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGa                     <  3:958765/148‑1 (MQ=255)
                                                                                                                               gACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAt                   >  1:51240/1‑149 (MQ=255)
                                                                                                                                gataaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAtt                  >  2:920107/4‑149 (MQ=255)
                                                                                                                                  aaaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc                <  2:1012034/149‑1 (MQ=255)
                                                                                                                                   aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc                >  4:944940/1‑148 (MQ=255)
                                                                                                                                    aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                    <  1:174421/95‑1 (MQ=255)
                                                                                                                                    aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                    >  2:174421/1‑95 (MQ=255)
                                                                                                                                    aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCAc              <  1:602358/149‑1 (MQ=255)
                                                                                                                                        gAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGttt          <  1:209891/149‑1 (MQ=255)
                                                                                                                                           ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTc                                  <  2:348601/122‑1 (MQ=255)
                                                                                                                                           ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTc                                  >  1:348601/1‑122 (MQ=255)
                                                                                                                                           ccccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCag       >  2:603366/1‑149 (MQ=255)
                                                                                                                                            cccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga      >  1:556304/1‑149 (MQ=255)
                                                                                                                                            cccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCag       <  3:944940/148‑1 (MQ=255)
                                                                                                                                             ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCagag     <  2:556304/149‑1 (MQ=255)
                                                                                                                                              cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCagag     >  3:867640/1‑148 (MQ=255)
                                                                                                                                              cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAgg    >  1:1046929/1‑149 (MQ=255)
                                                                                                                                               gAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAgg    >  3:321842/1‑148 (MQ=255)
                                                                                                                                                aaGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa   <  4:321842/148‑1 (MQ=255)
                                                                                                                                                aaGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAc  >  1:965956/1‑149 (MQ=255)
                                                                                                                                                 aGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAc  >  3:414151/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGAC  >  NC_000913/2731438‑2731730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑