Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,059,963 A→C 5.1% I37L (ATC→CTC)  iciA → chromosome initiation inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,059,9630AC 5.1% 248.0 / 4.3 78I37L (ATC→CTC) iciAchromosome initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (2/2);  total (40/38)
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.15e-01

AACAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGT  >  NC_002947/1059821‑1060099
                                                                                                                                               |                                                                                                                                          
aaCAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGt                                                                                                                                      >  6:54775/1‑149 (MQ=255)
aaCAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGt                                                                                                                                      <  3:195733/149‑1 (MQ=255)
aaCAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGt                                                                                                                                      >  3:103632/1‑149 (MQ=255)
 aCAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGtt                                                                                                                                     >  6:280347/1‑149 (MQ=255)
  cAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCCTCAAgcaa                                                                                                                                    >  3:6546/1‑146 (MQ=255)
   aaaTAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGc                                                                                                                                   >  2:245960/1‑149 (MQ=255)
     aTAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCtt                                                                                                                                 >  7:228821/1‑149 (MQ=255)
                cACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCgg                                                                                                                        <  4:64885/147‑1 (MQ=255)
                 aCCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCCTCAAGCAACTTGATGCGCtggg                                                                                                                      <  4:6546/148‑2 (MQ=255)
                 aCCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTc                                                                                                                     <  5:54775/149‑1 (MQ=255)
                  cccATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCg                                                                                                                    <  7:196965/149‑1 (MQ=255)
                  cccATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCg                                                                                                                    <  8:228821/149‑1 (MQ=255)
                   ccATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCgg                                                                                                                   >  4:19787/1‑149 (MQ=255)
                         ggAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGGAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCCGCGCCTCAAGCAACTTGATgcgctg                                                                                                                        <  6:26486/138‑3 (MQ=255)
                         ggAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCCGCGCCTCAAGCAACTTGATgcgctg                                                                                                                        >  5:26486/1‑136 (MQ=255)
                         ggAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGcc                                                                                                             <  6:109736/149‑1 (MQ=255)
                          gAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCg                                                                                                            >  2:140801/1‑149 (MQ=255)
                            aaaTCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGt                                                                                                          <  3:121134/149‑1 (MQ=255)
                             aaTCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTg                                                                                                         <  6:307656/149‑1 (MQ=255)
                                 gATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTgg                                                                                                   >  2:306290/1‑149 (MQ=255)
                                 gATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTg                                                                                                    <  4:103632/148‑1 (MQ=255)
                                       cGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGt                                                                                              <  3:19787/148‑1 (MQ=255)
                                        gACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGt                                                                                                  <  7:299609/143‑1 (MQ=255)
                                        gACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGc                                                                                            >  5:84705/1‑149 (MQ=255)
                                        gACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGCGTCGGCCAGACGGTGGTGGTaaga                                                                                              >  8:299609/1‑144 (MQ=255)
                                             cAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTcgc                                                                                       >  8:199882/1‑149 (MQ=255)
                                             cAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTcgc                                                                                       >  4:191540/1‑149 (MQ=255)
                                                 ctgctgGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGc                                                                                            >  5:50661/1‑140 (MQ=255)
                                                 ctgctgGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGc                                                                                            <  6:50661/140‑1 (MQ=255)
                                                 ctgctgGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCgccgcc                                                                                   >  6:71959/1‑149 (MQ=255)
                                                       gctgCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGgc                                                                               <  5:155301/147‑1 (MQ=255)
                                                              tCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCg                                                                        >  5:115235/1‑147 (MQ=255)
                                                              tCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCg                                                                        >  4:240771/1‑147 (MQ=255)
                                                                gcgGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAg                                                                     <  7:199882/148‑1 (MQ=255)
                                                                 cgGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTc                                                                   <  1:306290/149‑1 (MQ=255)
                                                                        gATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGgccgc                                                             >  5:71229/1‑148 (MQ=255)
                                                                         aTCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTTCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCa                                                           <  1:245960/149‑1 (MQ=255)
                                                                         aTCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGgccgc                                                             >  8:324809/1‑147 (MQ=255)
                                                                          tCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGgccgc                                                             >  1:240504/1‑147 (MQ=255)
                                                                           cGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGc                                                         >  5:240891/1‑149 (MQ=255)
                                                                            gAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGc                                                         <  6:240891/148‑1 (MQ=255)
                                                                             aaCAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGgcgc                                                                                                                          >  4:321996/1‑84 (MQ=255)
                                                                             aaCAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGgcgc                                                                                                                          <  3:321996/84‑1 (MQ=255)
                                                                                aaGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACt                                                     >  6:42200/1‑148 (MQ=255)
                                                                                  ggCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGaa                                                                      <  5:7387/129‑1 (MQ=255)
                                                                                  ggCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGaa                                                                      >  6:7387/1‑129 (MQ=255)
                                                                                  ggCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCaa                                                  >  1:233438/1‑149 (MQ=255)
                                                                                     ggTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCa                                               >  1:40123/1‑149 (MQ=255)
                                                                                       tttCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGCGGCGCGGGTCGGGCAGCCGGTGCTTGGTACGTGCCTCGCCGCGCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCAt                                              <  1:267668/149‑1 (MQ=255)
                                                                                       tttCGAGCGTGCGGAGCATGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGACGGTGCTGGGTACGTGCCTCGCCGCCCGGCGCGACCGAATTCGGCAGCCAGCTAATCAACCAt                                              <  1:263458/149‑1 (MQ=255)
                                                                                           gAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCAtgtg                                           >  5:217992/1‑147 (MQ=255)
                                                                                           gAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTgc                                          <  5:42200/148‑1 (MQ=255)
                                                                                            aGCGTGCGGCGCAGGTGATGGGCTTGTAGCAATCGGCGATTTCCCAGCGCATCAAGATGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTgcag                                        >  6:74016/1‑149 (MQ=255)
                                                                                             gCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGTGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTgcag                                        >  2:291765/1‑148 (MQ=255)
                                                                                                 gcggcgCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGt                                   <  8:277574/149‑1 (MQ=255)
                                                                                                 gcggcgCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGt                                   >  8:185695/1‑149 (MQ=255)
                                                                                                      gcAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGATGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGcc                              >  2:242712/1‑149 (MQ=255)
                                                                                                       cAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCt                             <  6:217992/149‑1 (MQ=255)
                                                                                                           tGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGcgc                         >  6:189186/1‑147 (MQ=255)
                                                                                                            gCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCg                        <  5:280347/149‑1 (MQ=255)
                                                                                                                ggCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGGCAGCTACTGAACCATGTGCAGCAGGTGCGCCTGCTCGAAc                     <  6:115235/148‑1 (MQ=255)
                                                                                                                ggCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAAc                     <  6:71229/148‑1 (MQ=255)
                                                                                                                  cTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAAcgc                   >  3:277406/1‑148 (MQ=255)
                                                                                                                       tCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGAcc               <  5:71959/148‑1 (MQ=255)
                                                                                                                       tCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCt              >  7:11247/1‑149 (MQ=255)
                                                                                                                       tCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCt              <  1:242712/149‑1 (MQ=255)
                                                                                                                         gCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGAcc               <  6:152176/146‑1 (MQ=255)
                                                                                                                         gCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGAcc               >  5:152176/1‑146 (MQ=255)
                                                                                                                         gCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTg             <  3:191540/148‑1 (MQ=255)
                                                                                                                           aaTCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACTGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCt              >  8:71433/1‑145 (MQ=255)
                                                                                                                           aaTCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACTGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCt              <  7:71433/145‑1 (MQ=255)
                                                                                                                            aTCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAgc         <  2:233438/149‑1 (MQ=255)
                                                                                                                               ggCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAgcgc       >  8:199161/1‑148 (MQ=255)
                                                                                                                               ggCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAgcgc       >  3:270559/1‑148 (MQ=255)
                                                                                                                               ggCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAgcgc       >  4:237100/1‑148 (MQ=255)
                                                                                                                                gCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCa     <  7:185695/149‑1 (MQ=255)
                                                                                                                                  gATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAgg   <  4:270559/149‑1 (MQ=255)
                                                                                                                                   aTTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGt  <  2:40123/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
AACAAATAAAGCTCTTCACCCATTAGGAAAATCGATGTTCGACTACAAACTGCTGGCTGCCCTCGCGGCAGTGATCGAACAAGGCGGTTTCGAGCGTGCGGCGCAGGTGCTGGGCTTGTCGCAATCGGCGATTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGT  >  NC_002947/1059821‑1060099

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