Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,986,9950AC42.2% 186.4 / 23.2 103Y24S (TAC→TCC) ftnBferritin B, putative ferrous iron reservoir
Reads supporting (aligned to +/- strand):  ref base A (45/14);  new base C (0/43);  total (45/58)
Fisher's exact test for biased strand distribution p-value = 9.54e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCC  >  NC_000913/1986846‑1987140
                                                                                                                                                       |                                                                                                                                                 
   attattCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCtt                                                                                                                                               >  4:184063/1‑151 (MQ=255)
   attattCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCtt                                                                                                                                               <  1:1229955/151‑1 (MQ=255)
      attCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCAc                                                                                                                                            >  3:167342/1‑151 (MQ=255)
      attCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCAc                                                                                                                                            >  1:523694/1‑151 (MQ=255)
          tGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTg                                                                                                                                         >  1:99367/1‑150 (MQ=255)
            cAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACTGATGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCGACCTTCACCTGa                                                                                                                                        >  4:569042/1‑149 (MQ=255)
             aaaCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTa                                                                                                                                     <  3:62706/151‑1 (MQ=255)
              aaCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTaa                                                                                                                                    <  4:602779/151‑1 (MQ=255)
               aCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAAc                                                                                                                                   <  4:309729/151‑1 (MQ=255)
                cAGTAATTCTGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACt                                                                                                                                  >  1:898845/1‑151 (MQ=255)
                cAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACt                                                                                                                                  <  4:1342097/151‑1 (MQ=255)
                 aGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTg                                                                                                                                 <  4:63806/151‑1 (MQ=255)
                     aTTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGgtgt                                                                                                                             >  4:1473243/1‑151 (MQ=255)
                           gtgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGaa                                                                                                                       <  4:468441/151‑1 (MQ=255)
                           gtgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGa                                                                                                                        >  4:593827/1‑150 (MQ=255)
                                    tttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTg                                                                                                              >  1:482714/1‑151 (MQ=255)
                                      tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGa                                                                                                             >  4:369625/1‑150 (MQ=255)
                                      tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGa                                                                                                             >  4:1205721/1‑150 (MQ=255)
                                      tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGa                                                                                                             >  2:75176/1‑150 (MQ=255)
                                      tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGa                                                                                                             >  1:302535/1‑150 (MQ=255)
                                       aTCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAc                                                                                                           >  1:747532/1‑151 (MQ=255)
                                        tCTTGCACCTCTCCACTTCTGGATATAAGGCTATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACg                                                                                                          >  4:1280657/1‑151 (MQ=255)
                                           tGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAc                                                                                                           <  2:10564/147‑1 (MQ=255)
                                           tGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGc                                                                                                        <  4:344851/150‑1 (MQ=255)
                                               ccTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGc                                                                                                   >  1:1074962/1‑151 (MQ=255)
                                                  ctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCAGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGcc                                                                                                  <  1:1012631/149‑1 (MQ=255)
                                                   tcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACt                                                                                               >  4:591583/1‑151 (MQ=255)
                                                    ccACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                              <  1:621691/151‑1 (MQ=255)
                                                    ccACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCa                                                                                                 <  4:1238404/148‑1 (MQ=255)
                                                         tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCtt                                                                                         >  1:543692/1‑151 (MQ=255)
                                                         tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCtt                                                                                         >  1:987681/1‑151 (MQ=255)
                                                           tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTcg                                                                                       >  1:967330/1‑151 (MQ=255)
                                                              ataaAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTAAACCTCAACTCAAAAATGAC‑CCGCGATTTTTA‑CGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTcgcg                                                                                     <  1:781980/150‑1 (MQ=255)
                                                                ataAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGccc                                                                                  >  4:940544/1‑151 (MQ=255)
                                                                 taAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGccc                                                                                  >  4:928471/1‑150 (MQ=255)
                                                                  aaGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAg                                                                                <  1:1054615/151‑1 (MQ=255)
                                                                   aGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAgg                                                                               >  4:148084/1‑151 (MQ=255)
                                                                      atatTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGcac                                                                            >  1:263084/1‑151 (MQ=255)
                                                                         ttAGATATGGCAACCGCTGAAATGCTTCTCAACCTCAACTCAAAAAGGAC‑CCGCGATTTTTA‑CGCTCCCAACCTCCCCCTTCCCTTGAGTACCGGGTGTCCTAAACAGATTCTGACGGGCACCGCCACTTTCCTTCGCGCCCAGGcaca                                                                           <  4:785192/149‑1 (MQ=255)
                                                                           aGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       >  1:810381/1‑151 (MQ=255)
                                                                            gttATGGCAACCGCTGGAATGCTTCTCAAACTCACCCCTCAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTCCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       <  1:198742/148‑1 (MQ=255)
                                                                            ggTATGGCACCCGCTGGAATGCTTCCCAACCTCACCTCTCAAATGAC‑CCGCGAGTTTTA‑CGCACCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTTTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       <  1:920518/150‑1 (MQ=255)
                                                                            ggTAGGGAACCCGCTGGAATGCTTCCCAAACTCACCTCTCAAATGAA‑CCGCGAGTTTTA‑CGCTTCCAATCTCTCCCTCCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAGGGGCACCGCCACTTTCCTTCGCGCCCAGGCACagag                                                                        <  4:1472472/149‑1 (MQ=255)
                                                                              tATGGCAACCGCTGGAATGCTTCTCAAACTCACCTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTa                                                                      >  1:764713/1‑149 (MQ=255)
                                                                              tATGGCAACCGCTGGAATGCTTCTAAAACTCACCCCTAAAATGAC‑CCGCGATTTTTA‑CGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTAAACAGAGTCTAACCGGCGCCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       <  4:766836/148‑1 (MQ=255)
                                                                                   cacccGCTGGAATGTTTCTAAAACCCAACTCCCAAAGAAC‑CCGCGATTTTTA‑CGCATCAAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGAc                                                               <  1:366026/148‑1 (MQ=255)
                                                                                    aaCGGCGGGATTGCTTCCCAACCTCACCCTGAAAATGAG‑CGGCG‑GTTTTATGGCACCCAATCTCTCCCTTCCCCTAGGTAACTGGTTTTCCGACAAGAGTCTAACCGGGACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGa                                                                <  2:140540/149‑1 (MQ=255)
                                                                                    aaCCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGAcc                                                              >  1:1466503/1‑151 (MQ=255)
                                                                                        gCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCaa                                                           >  4:1307553/1‑150 (MQ=255)
                                                                                         cTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAat                                                         >  1:439486/1‑151 (MQ=255)
                                                                                           ggAATGCTTCTCAAACTCACCTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatg                                                        <  1:998226/150‑1 (MQ=255)
                                                                                           ggAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGATTTTTA‑CGCATCCAATCTCTCCCTTCACCTGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCaaa                                                          <  1:329250/148‑1 (MQ=255)
                                                                                           ggAATGCTTCTAAAACTCAACTCTCAAATGAC‑CCGCGAGTTTTA‑CGCATCAAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatga                                                       <  1:575271/151‑1 (MQ=255)
                                                                                             aaTGCTCCTCAAACTCAACTCTCAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatgat                                                      <  1:308550/150‑1 (MQ=255)
                                                                                               ttttccccAACCTCACCTCT‑AAACGAACCCGCGTGTTTTA‑CGCACCCATTCGCCCCCTCCCCCGGATTACTGGGTGTTCGGAACAGATTCGGAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                    <  4:1460563/143‑1 (MQ=255)
                                                                                               tGCTTCTCAACCTCAACTCTCAAATGAC‑CCGCGATTTTTA‑CGCATCCAATCTCTCCCTCCACCGGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                    <  4:682973/150‑1 (MQ=255)
                                                                                                gCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcg                                                   <  1:177075/150‑1 (MQ=255)
                                                                                                gCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcgc                                                  >  4:1358218/1‑151 (MQ=255)
                                                                                                  ttCTCAACCTCAACTCTCAAATGAC‑CCGGGATTTTTA‑CGCACCCATTCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcg                                                   <  4:135329/148‑1 (MQ=255)
                                                                                                  ttCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCAt                                                >  4:798166/1‑151 (MQ=255)
                                                                                                   tctcAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                               >  4:54195/1‑151 (MQ=255)
                                                                                                   tctcAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                               >  4:1406306/1‑151 (MQ=255)
                                                                                                    ctcAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGt                                              >  4:1126765/1‑151 (MQ=255)
                                                                                                       aaccccaactcTAAAAGGAG‑CCGCGAGTTTTA‑CCCACCCAATCTCTCCCTCCCCCTAAGTAACTGGTTTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTa                                           <  4:737514/146‑1 (MQ=255)
                                                                                                       aaaCTCACCTCTCAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGtt                                             <  1:298569/149‑1 (MQ=255)
                                                                                                       aaaCTCAACTCTCAAATGAC‑CCGCGAGTTTTA‑CGCACCCAATCTCTCCCTTCCCCGGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCa                                                 <  1:696871/145‑1 (MQ=255)
                                                                                                        aactcaactcTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTa                                           >  1:1096953/1‑150 (MQ=255)
                                                                                                        aacccacctcTCAAATCAA‑CCGCGATTTTCG‑CGCTTCCAATCTCTCCCTTCCCCTGAGTAGCGGGTGTTTTGAACGGAGTCTGAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTaa                                          <  4:1075334/151‑1 (MQ=255)
                                                                                                          ctcaactcTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACt                                        >  4:1264176/1‑151 (MQ=255)
                                                                                                            cacctctAAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTCCCCCTGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACtt                                       <  4:125025/147‑1 (MQ=255)
                                                                                                            caactcTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTTCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACt                                        <  1:884517/149‑1 (MQ=255)
                                                                                                             acctctAAAATAAC‑CCGCGATTTTTA‑CCCTTCCATTCTCTCCCTCCCCCTGAGTAACGGGTTTTCTAAAAAGAGTCTGACCGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACtt                                       <  4:926623/147‑1 (MQ=255)
                                                                                                             aactcTCAAATGAC‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTa                                     <  1:539093/151‑1 (MQ=255)
                                                                                                               ctctCAAATGAC‑CCGCGAGTTTTA‑CGCACCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTa                                     <  1:738854/149‑1 (MQ=255)
                                                                                                               ccctcAAATGAA‑CCGCGAGTTTAA‑CGCACCCAATCCCGCCCTCCCCCGGAGTAACGGTTGTTCTGAACAGAGTATGACCGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACt                                        <  4:446394/144‑1 (MQ=255)
                                                                                                                tctcAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATg                                   >  1:463455/1‑150 (MQ=255)
                                                                                                                    aaaTGAC‑CCGCGATTTTTA‑CGCACCCATTCTCTCCCTCCCCCTGAGTAACGGGTTTTCTGAACAGAGTCTGAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAag                                <  4:651041/149‑1 (MQ=255)
                                                                                                                    aaaTGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAagag                              <  1:517033/151‑1 (MQ=255)
                                                                                                                     aaTGAA‑CCGCGATTTTTA‑CGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                             <  1:346179/151‑1 (MQ=255)
                                                                                                                      aaaaacccGCGAGTTTTA‑CGCACCCATTCTCTCCCTTCCCCGGAGTACCTGGTGTTCTGAACAGAGTCTGAGCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                             <  4:1454250/145‑1 (MQ=255)
                                                                                                                            ccGCGATTTTTA‑CGCACCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgc                       <  1:628257/151‑1 (MQ=255)
                                                                                                                             cgggatttttA‑CGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTAAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgcg                      <  4:209957/145‑1 (MQ=255)
                                                                                                                             cgcgAGTTTTA‑CGCATCCAATCTCTACCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgcg                      <  1:1121695/151‑1 (MQ=255)
                                                                                                                              gggggTTTTA‑CGCCTCCAATCTCTCCCTTCCCCGGGGTAACTGGTGTTCAGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgcg                      <  4:1144233/146‑1 (MQ=255)
                                                                                                                              gcgATTTTTA‑CGCATCCAATCTCTCCCTCCCCCTGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgc                       <  1:589847/149‑1 (MQ=255)
                                                                                                                               cgAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCt                    >  1:515987/1‑151 (MQ=255)
                                                                                                                                gAGTTTTA‑CCCATCCAATCTCTTCCTCCCCCGGAGTAACTGGTGTTCTGAACAGAGTCTGAGCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCgg                        <  2:35663/146‑1 (MQ=255)
                                                                                                                                  gTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAcc                 >  1:32156/1‑151 (MQ=255)
                                                                                                                                  gTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAcc                 >  3:184802/1‑151 (MQ=255)
                                                                                                                                   ttttC‑CGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAc                  <  1:1136348/149‑1 (MQ=255)
                                                                                                                                    ttttccccctCCAATCTCTCCCTCCCCCTGATTACCTGGTGTTCTAACCAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTa                   <  1:721817/139‑1 (MQ=255)
                                                                                                                                    tttA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAcccc               >  4:1463584/1‑151 (MQ=255)
                                                                                                                                     ttA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAccccc              >  4:1185692/1‑151 (MQ=255)
                                                                                                                                         cGCACCCAATCTCTCCCTTCACCGGAGTAACGGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATc           <  4:945691/151‑1 (MQ=255)
                                                                                                                                              ccAATCTCTCCCTTCCCCTGATTAACTGGTGTTCTGAACAGAGTCGGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTa       <  4:72509/150‑1 (MQ=255)
                                                                                                                                              ccAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGtt        <  4:619581/149‑1 (MQ=255)
                                                                                                                                                aaTCTCTCCCTTCACCTGAGTACCTGGTTTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTaaa     <  1:369657/150‑1 (MQ=255)
                                                                                                                                                  tctcTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTGTGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGcc  >  1:783279/1‑151 (MQ=255)
                                                                                                                                                    tctcCCTTCCCCTGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAg    <  4:851273/147‑1 (MQ=255)
                                                                                                                                                     ctACCTTCCCCTGATAAACTGTTTTTCTGAACAGAGTCTGAACGTCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGc   <  4:754325/147‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAA‑CCGCGAGTTTTA‑CGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCC  >  NC_000913/1986846‑1987140

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