Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,582 A→C 9.0% intergenic (‑425/+18) 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,5820AC 9.0% 637.8 / 51.4 201intergenic (‑425/+18)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base A (93/89);  new base C (10/8);  total (104/97)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 209 total aligned reads displayed.

TTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731434‑2731729
                                                                                                                                                    |                                                                                                                                                   
 tgttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGaa                                                                                                                                                    <  1:508720/149‑1 (MQ=255)
 tgttCCGTGTCACTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCAAAAAATATAATAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGaa                                                                                                                                                    <  3:1582755/149‑1 (MQ=255)
  gttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAg                                                                                                                                                   <  3:862435/149‑1 (MQ=255)
  gttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAg                                                                                                                                                   <  4:589841/149‑1 (MQ=255)
   ttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAgg                                                                                                                                                  >  3:1000563/1‑149 (MQ=255)
    tCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAggg                                                                                                                                                 >  1:1417882/1‑149 (MQ=255)
         gtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCAGAAGGGCTCGt                                                                                                                                            <  1:1300423/149‑1 (MQ=255)
          tCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtt                                                                                                                                           <  4:800441/149‑1 (MQ=255)
           cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttt                                                                                                                                          <  1:193095/149‑1 (MQ=255)
             gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTa                                                                                                                                        <  2:729602/149‑1 (MQ=255)
               ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttttgt                                                                                                                                      <  1:148331/149‑4 (MQ=39)
               ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATc                                                                                                                                      <  4:420237/149‑1 (MQ=255)
                gtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCa                                                                                                                                     >  4:606859/1‑149 (MQ=255)
                  ggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCAtt                                                                                                                                   <  4:724997/149‑1 (MQ=255)
                     gcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGttttt                                                                                                                                        <  1:1242307/141‑2 (MQ=255)
                      cATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                               <  3:1289843/149‑1 (MQ=255)
                             gggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATGGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGc                                                                                                                        >  2:381545/1‑149 (MQ=255)
                                 gTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtt                                                                                                                                           >  1:980731/1‑126 (MQ=255)
                                 gTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACa                                                                                                                    >  2:395474/1‑149 (MQ=255)
                                  ttATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACaa                                                                                                                   >  1:269637/1‑149 (MQ=255)
                                    aTTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATc                                                                                                                 >  1:118403/1‑149 (MQ=255)
                                     ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACgg                                                                                                                         <  4:566796/140‑1 (MQ=255)
                                     ttCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACgg                                                                                                                         >  3:566796/1‑140 (MQ=255)
                                      tCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCg                                                                                                               <  3:451978/149‑1 (MQ=255)
                                        cGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGt                                                                                                             <  3:179437/149‑1 (MQ=255)
                                            cTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTc                                                                                                         <  1:1221209/149‑1 (MQ=255)
                                                 aaGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                    >  2:632844/1‑149 (MQ=255)
                                                 aaGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTaa                                                                                                    >  2:952634/1‑149 (MQ=255)
                                                     gAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTc                                                                                                >  1:706893/1‑149 (MQ=255)
                                                     gAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTc                                                                                                >  2:1493945/1‑149 (MQ=255)
                                                       aaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                              >  3:257087/1‑149 (MQ=255)
                                                       aaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                              <  2:1417882/149‑1 (MQ=255)
                                                         aaTATAAAAACGTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAgg                                                                                            >  1:916112/1‑149 (MQ=255)
                                                              aaaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                       >  2:877813/1‑149 (MQ=255)
                                                               aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  4:1417199/1‑149 (MQ=255)
                                                               aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATc                                                                                      >  1:1030738/1‑149 (MQ=255)
                                                                 aaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCg                                                                                                                       <  2:255513/114‑1 (MQ=255)
                                                                 aaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCAc                                                                                    >  4:1285150/1‑149 (MQ=255)
                                                                    tttATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACttt                                                                                 >  1:1543316/1‑149 (MQ=255)
                                                                        tCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTAccc                                                                             <  2:1607317/149‑1 (MQ=255)
                                                                          gTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttt                                                                                                                                         <  1:36202/87‑1 (MQ=38)
                                                                          gTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCgg                                                                           >  4:1093141/1‑149 (MQ=255)
                                                                            tCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGaa                                                                         <  2:1543316/149‑1 (MQ=255)
                                                                                tCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCaa                                                                     <  1:441865/149‑1 (MQ=255)
                                                                                  aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGaa                                                                                                                                                    <  3:225301/68‑1 (MQ=255)
                                                                                  aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGaa                                                                                                                                                    >  4:225301/1‑68 (MQ=255)
                                                                                  aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                   <  3:1487766/149‑1 (MQ=255)
                                                                                  aCTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                   <  4:1000563/149‑1 (MQ=255)
                                                                                         aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGa                                                            <  3:122294/149‑1 (MQ=255)
                                                                                         aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGa                                                            <  4:1668520/149‑1 (MQ=255)
                                                                                             aaaaCATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGt                                                        >  1:863656/1‑149 (MQ=255)
                                                                                                     ttAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                >  4:730122/1‑149 (MQ=255)
                                                                                                     ttAAATAGAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtg                                                >  2:1659234/1‑149 (MQ=255)
                                                                                                      gaaaTATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  2:1312782/148‑1 (MQ=255)
                                                                                                        aaTATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                               <  3:79598/67‑1 (MQ=255)
                                                                                                        aaTATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGc                                             >  1:202208/1‑149 (MQ=255)
                                                                                                            taGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAc                                                                        >  3:1645202/1‑118 (MQ=255)
                                                                                                            taGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAc                                                                        <  4:1645202/118‑1 (MQ=255)
                                                                                                             aGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAAc                                                                                                   <  2:1130804/90‑1 (MQ=255)
                                                                                                             aGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAAc                                                                                                   >  1:1130804/1‑90 (MQ=255)
                                                                                                               tCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCggg                                      <  2:202208/149‑1 (MQ=255)
                                                                                                                cTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGt                                     <  1:1597052/149‑1 (MQ=255)
                                                                                                                 ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtt                                    >  3:1155081/1‑149 (MQ=255)
                                                                                                                 ttttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtt                                    >  2:1052211/1‑149 (MQ=255)
                                                                                                                  tttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGttt                                   <  4:827390/149‑1 (MQ=255)
                                                                                                                   ttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTa                                                       >  1:1055950/1‑128 (MQ=255)
                                                                                                                       gTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              <  1:240721/149‑1 (MQ=255)
                                                                                                                       gTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              >  3:42309/1‑149 (MQ=255)
                                                                                                                         cTAACTTATAGACAAAAACGAGCCCCGTAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTATTCATTAACTTCCAGCCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATct                            >  4:1213819/1‑149 (MQ=255)
                                                                                                                         cTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATct                            <  1:395474/149‑1 (MQ=255)
                                                                                                                         cTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATct                            >  4:689488/1‑149 (MQ=255)
                                                                                                                          tAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctg                           >  3:1105510/1‑149 (MQ=255)
                                                                                                                           aaCTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                          <  3:730122/149‑1 (MQ=255)
                                                                                                                              ttATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAtt                                                                                                      <  3:305020/70‑1 (MQ=255)
                                                                                                                              ttATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                       >  2:12248/1‑149 (MQ=255)
                                                                                                                              ttATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                       >  1:467870/1‑149 (MQ=255)
                                                                                                                              ttATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                       >  4:159695/1‑149 (MQ=255)
                                                                                                                               tataGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgct                      <  3:1417199/149‑1 (MQ=255)
                                                                                                                                ataGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgctg                     <  3:693654/149‑1 (MQ=255)
                                                                                                                                  aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGaa                   >  2:1060875/1‑149 (MQ=255)
                                                                                                                                  aGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGaa                   <  4:1497052/149‑1 (MQ=255)
                                                                                                                                   gACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCtgt                                               <  1:1446962/120‑1 (MQ=255)
                                                                                                                                    aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTAttt                                                    >  2:749746/1‑114 (MQ=255)
                                                                                                                                    aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAtt                 >  4:700420/1‑149 (MQ=255)
                                                                                                                                     cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                >  2:472281/1‑149 (MQ=255)
                                                                                                                                     ataaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTg                >  1:1724489/3‑149 (MQ=255)
                                                                                                                                       aaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGTTGCTGAATTGCa              >  3:59632/1‑149 (MQ=255)
                                                                                                                                       aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              <  3:1542705/149‑1 (MQ=255)
                                                                                                                                       aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              <  4:207935/149‑1 (MQ=255)
                                                                                                                                       aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              <  2:334923/149‑1 (MQ=255)
                                                                                                                                        aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCAc             <  3:947474/149‑1 (MQ=255)
                                                                                                                                        aaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCAc             >  1:842983/1‑149 (MQ=255)
                                                                                                                                         aaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACg            >  3:905544/1‑149 (MQ=255)
                                                                                                                                              gCCCCGAAGGGCTCGTTTTATCATTTCCTGGACGGCAACAATCCGTTCTTCATTAACTTCCAGTCGAATCACTTTAACCGGACCCAATTCACCAGACAGTATTTGCTGTGCCAGCCGGTTGTCGATCTGCTGCTGAATTGCACGTTTCa       <  4:1508980/149‑1 (MQ=255)
                                                                                                                                                 ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAtt                                                                                                      <  4:1718693/51‑1 (MQ=255)
                                                                                                                                                 ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAtt                                                                                                      >  3:1718693/1‑51 (MQ=255)
                                                                                                                                                 ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAGCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              <  3:1361333/123‑1 (MQ=255)
                                                                                                                                                 ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              >  4:1361333/1‑123 (MQ=255)
                                                                                                                                                 ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCaga     >  4:873156/1‑148 (MQ=255)
                                                                                                                                                   gAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa  >  1:318221/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGTTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731434‑2731729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: