Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,112,5990TG43.4% 142.1 / 21.1 85G95G (GGT→GGGyiiRDUF805 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (11/36);  new base G (36/0);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 6.84e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAG  >  NC_000913/4112450‑4112748
                                                                                                                                                            |                                                                                                                                                            
aCCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGgtg                                                                                                                                                             <  2:1231603/151‑1 (MQ=255)
  cgcgGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTGCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGgt                                                                                                                                                              <  2:347805/148‑1 (MQ=255)
   gcggcgTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCt                                                                                                                                                         <  2:290893/151‑1 (MQ=255)
          tttGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAAt                                                                                                                                                <  1:1103194/150‑1 (MQ=255)
            tGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATgg                                                                                                                                             <  2:230733/150‑1 (MQ=255)
                 tgtgtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGTTTTACCGGGGGTCTGGCAATGGGCgggg                                                                                                                                       >  1:643834/1‑149 (MQ=255)
                    gtgCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTGATCTGGGGGTTTTTTCGGGGGTTTGGCAATTgggggggggg                                                                                                                                    >  1:680473/1‑151 (MQ=255)
                     tgcttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTTGGGGATTTTTCCCGGGGTCTGGCAATGgggggggggg                                                                                                                                    >  1:631296/1‑150 (MQ=255)
                     tgcttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGGGGTTTTTCCGGGGGGTTGGCAATGgggggggggg                                                                                                                                    >  1:324115/1‑150 (MQ=255)
                       cttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACCGGGGGGTTTTTCCCGGGGTCTGGGAAATGGGGGGGGGGCgg                                                                                                                                 >  2:820641/1‑150 (MQ=255)
                        ttgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAAAGGGGGGGTTTTACCGGGGTCTTGCAAAGGGCGGGGGGGCg                                                                                                                                  >  2:684798/1‑149 (MQ=255)
                         tgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGtt                                                                                                                                <  2:820872/150‑1 (MQ=255)
                         tgctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGt                                                                                                                                 <  2:213769/149‑1 (MQ=255)
                          gctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCGGTTTTTCCCGGGGTCTGGCAATGGGGGGGGGGGGGgttt                                                                                                                              >  1:956104/1‑151 (MQ=255)
                          gctCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGGGATTTTTCCCGGGGTCTGGCAATGGGCGGGGGGGGGgttt                                                                                                                              >  1:594411/1‑151 (MQ=255)
                             ctGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTg                                                                                                                             <  2:246930/149‑1 (MQ=255)
                                    aCAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCAGTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGGTTTTACCGGGTGTCTGGCAATGGGCGTTAAGGCGTTTTGTCCCGACg                                                                                                                    <  2:701877/150‑1 (MQ=255)
                                     cAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGtt                                                                                                                  <  2:682124/151‑1 (MQ=255)
                                       gCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACCGGGGGGTTTTTCCGGGGGGTTGGCAATGGGGGGGGGGGGGTTTTGTTCCggc                                                                                                                      >  2:848095/1‑144 (MQ=255)
                                          gCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTc                                                                                                             <  1:1053165/151‑1 (MQ=255)
                                            aGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGGGATTTTTCCGGGTGTCTGGCGATGGGGGGTGGGGGGTTTTGTTCCGGCGGTGtttat                                                                                                            >  1:448339/1‑148 (MQ=255)
                                            aGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAACGGGGGGTTTTTCCCGGGGGCTGGGAAATGGCGGGGGGGGGTTTTTTccccgcctttgttt                                                                                                              >  2:3951/1‑138 (MQ=255)
                                              tAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACAGGGGGGATTTACCGGGGGGGTGGGAAATGGGGGGGGGGGGGTTTGGTCCGACGTTGGTTCTgg                                                                                                          >  2:304979/1‑150 (MQ=255)
                                               aaCTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGt                                                                                                         <  1:341949/150‑1 (MQ=255)
                                                aCTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGGTGCTGCCGGCGGGGAAGCTGGGGCTTTTAACCGGGGGTTGGCCATGGGGGGTGGGGGGGTTTTTTCCCACggctctttctgggg                                                                                                       >  2:182278/1‑138 (MQ=255)
                                                  tGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCGATTTTTCCGGGGGTCTGGGAATGGGCGGGGGGGGGGTTTTTTCCGACCGTGAtttttgggat                                                                                                      >  1:1048268/1‑142 (MQ=255)
                                                  tGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGGGGTTTTTCCCGGGGTCTTGGCAAGGGGGGTGGGGGGGTTTTTCCCCACGTTGGTTTTtgtggt                                                                                                      >  1:73932/1‑148 (MQ=255)
                                                     tAAGCGCTTGCATGATCGCTGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAGTGGGCGGTGGGGCGGTTTGTCCCGACGTTGATTCTGGTgatgat                                                                                                   <  2:825807/150‑1 (MQ=255)
                                                     tAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCGTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTgatgat                                                                                                   <  2:974386/150‑1 (MQ=255)
                                                      aaGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGCGGTTTTTACGGGGGGTTGGCAATGGGGGGGGGGGGGGTTTTGCCCGACGTTTGTTCTGGGGGTGAAGc                                                                                                 >  2:535843/1‑151 (MQ=255)
                                                      aaGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATCGTGGCGTGGATGCTGCTGGCGGGGAACCGGGGGAATTTTCCGGGGGGCCGGCGAAGGGCGGGGGGGGGGTTTGTTCCGGCGGTGGTTCTgggg                                                                                                        >  1:1156953/1‑142 (MQ=255)
                                                           cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGTTTTTCCGGGGGTCTGGCAATGGGGGGGGGGGGGGTTTTTCCCCACGTTGTTTCTGGGGATGATGCTTAt                                                                                             >  1:188725/1‑150 (MQ=255)
                                                           cTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATTGGGGGTTTTAACGGGGGTCTGGGAATGGGCGGGGGGGGGGTTTGTCCCGACGGTGATTTTtgtggttgtggtt                                                                                               >  1:1037330/1‑139 (MQ=255)
                                                             tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATc                                                                                            <  2:689676/149‑1 (MQ=255)
                                                             tGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGAGTTTAACCGGGGGCTGGGCAAGGGCGGGGGGGGGTTTTTTCCCGACTTTGGTTCTTGGGGTGATGCTTTTCgg                                                                                          >  2:426462/1‑150 (MQ=255)
                                                                   aTCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTgg                                                                                    <  1:973680/151‑1 (MQ=255)
                                                                    tCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGATTTTTCGGGGGTTTGGGAATTGGGGGGGGGGGGGTTTTTTCCCCCGTTGATTTTTGGGGTGATTGTTATCGaaatt                                                                                      >  1:25629/1‑144 (MQ=255)
                                                                       cgGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCg                                                                                <  1:854127/151‑1 (MQ=255)
                                                                        gggCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGTTTTTCCGGGGGTCTGGGCATGGGGGGGGGGGGGGTTTTTTCCCGCGGTTATTTTGGGGATGGTGCTTATCGATCtttg                                                                                    >  1:253726/1‑144 (MQ=255)
                                                                           cGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACCTGGGGGTTTTTACCGGGGGCTGGGAAATGGGGGGGGGGGGGTTTTTTCCCGCGGTTGTTTTTGGGGTGATGGTTATTGGTATTGGGGCGttt                                                                             >  1:164943/1‑150 (MQ=255)
                                                                           cGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTttg                                                                            <  1:1143700/151‑1 (MQ=255)
                                                                            gTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCGATTTTTCCCGGGGTCTGGCAATGGGCGGGGGGGGGGTTTGTCCCGACGTTTATTTTGGGGGTGATGCTTTTCGGTTTTTGTGGGGttgt                                                                           >  1:888021/1‑151 (MQ=255)
                                                                            gTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACCGGGCGGTTTTATCCGGGGTTCGGCAAAGGGCGGGGGGGGGGTTTTTTCCCTCCTTGGTTTTGGGGATGATGGTTTTTCATCTTTGGGGGttt                                                                             >  1:525316/1‑149 (MQ=255)
                                                                              tCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTttgtt                                                                          <  2:119684/150‑1 (MQ=255)
                                                                               ccGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTgg                                                                        <  2:119053/151‑1 (MQ=255)
                                                                               ccGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGTTTTTACGGGGGTCTGGCCATGGGGGGGGGGGGGGTTTTTTCCGACGGTGGTTCTGGGGATGATGCTTTTCGATCTTTGTGCGGttttttg                                                                        >  1:687476/1‑149 (MQ=255)
                                                                                cGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGGGGTTTTACCCGGGGGCTGGCAATGGGGGGGGGGGCGGTTTTTCCCGACGGTGATTTTGGGGGTGATGCTTATCGATCTTTGTGGGTttgtttg                                                                        >  1:1125424/1‑148 (MQ=255)
                                                                                cGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAACTGGGCGGTTTTTCCCGGTGGTTGGCAATGGGGGGGGGGGGGGTTTTTTCCCCCGGTGTTTTTGGTGGTGGTGGTTTTTGATCTTGGGGCGgttgtt                                                                          >  2:6157/1‑148 (MQ=255)
                                                                                 ggCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGca                                                                      <  2:376276/151‑1 (MQ=255)
                                                                                   cgcATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAACGGGGGGTTTTTCCGGGGGGCTGGGAAAGGGGGGGGGGGGGGTTTTTTCCCGCGGTTATTTTGGTGGTGGTGGTTTTCGGTCTTTGGGCGGTTGTTGGca                                                                      >  2:279822/1‑149 (MQ=255)
                                                                                      aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACaa                                                                <  1:823777/151‑1 (MQ=255)
                                                                                      aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACaa                                                                <  1:585588/151‑1 (MQ=255)
                                                                                      aTGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGGACTGGGGGGATTTTCCGGGGGGCTGGCAATGGGGGGGGGGGCGGTTTGTTCCCACGGTTGTTTTTGGGGTGATGGTTATTGTTTTTGGTGCGTTTTTTTGgacac                                                                   >  1:1068144/1‑149 (MQ=255)
                                                                                        gggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAg                                                               <  2:1125404/150‑1 (MQ=255)
                                                                                         ggCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGc                                                             <  1:1219990/151‑1 (MQ=255)
                                                                                          gCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTTAACTTGGGAGTTTTCCGGGGGTCTGGCCAAGGGGGGGGGGGGGGTTTTTCCCCAACGTTTTTTTTGTGGTGGTGGTTTTCGATCTTTGGGGGGTTTTTGGCATACCAgggg                                                            >  1:804579/1‑149 (MQ=255)
                                                                                            aTTTCTGATGATTGTGGCGTGGATGCTGCTGGGGGGTAACATGTGGAGTTTTCCCGGGGGCCTGCAACTGGGGGGGGGGGGGTTTGTTCCGAGCTTGGTTTTGGTGGAGGGGCTCGTCGGTCTTTGTGGGGTTGTTGGGAATCAAgg                                                              >  1:353043/1‑147 (MQ=255)
                                                                                              ttCTGATGATTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGCGCTTTTTCCCGGGGGCTGGGCATTGGCGGGGGGGGGGTTTTTCCCCGCGTTTATTTTGGGGGTGGTGGTTTTTGGTCTTTGGGGGGTTGTTGGgaccccaagcggaa                                                          >  2:933373/1‑139 (MQ=255)
                                                                                                  gatgatTGTGGCGTGGATGCTGCTGGCGGGGAACCGGGGGGATTTTACCGGGGTCTGGGGATGGGGGGGGGGGGGGTTTTTGCCGACGGTGGTTCTTGGGGTGGTGGTTATCGGTCTTTGTGGGGTTGTTGGGACACAAAGGGGaaaaaa                                                      >  2:358962/1‑148 (MQ=255)
                                                                                                   atgatTGTGGCGTGGATGCTGCTGGCGGGGAACCGGGGGATTTTTCCCGGGGTCTGGCAATGGGGGGGGGGGGGGTTTTTCCCGGCGGTGATTCTGGGGGTGATGCTTTTCGATTTTTGTGCGTTTTTTTGGACACAAGGCGGaaataaa                                                    >  1:1185195/1‑150 (MQ=255)
                                                                                                    tgatTGTGGCGTGGATGCTGCTGTCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAaaataaat                                                   <  1:1219480/150‑1 (MQ=255)
                                                                                                    tgatTGTGGCGTGGATGCTGCTGGCGGGTAAATTGGCGGTTTTTACCGGGGGGTTGGAATGGGGGGGGGGGGGGTTTTTCCCCGCGTTGGTTTTTGGGGTGGTGGTTTTCGGTCTTTGGGGGGTTGTTGGCACACCAAGCGGaaataaa                                                    >  2:99917/1‑149 (MQ=255)
                                                                                                     gatTGTGGCGTGGATGCTGCTGGCGGGGAACTAGGGCGTTTTTCCGGGGGTGTGGGAAAGGGGGGGGGGGGGGTTTTTCCCGGCGGTGAATTTTGGGGTGATGGTTTTTGATTTTGGTGCGGTTGGTGGCACACAAGGGGAaaataaaaa                                                  >  2:179226/1‑148 (MQ=255)
                                                                                                     gatTGTGGCGTGGATGCTGCTGGCGGGGAAATGGGCGGTTTTTCCGGGGGTCTGGCAATGGGGGGGGGGGGGGTTTTTTCCCGCGTTGATTTTGGGGGTGGTGCTTTTCGGTCTTGGGGCGGTTTTTGGGACACAAGGGGGaaataaa                                                    >  1:875558/1‑148 (MQ=255)
                                                                                                        tgtgGCGTGGATGCTGCTGGCGGGGAAATGGGGGGTTTTTACGGGGGTCTTGGAATGGGGGGGGGGGGGGTTTTTTCCCACGGTTATTTTTGGGGTGGTGCTTTTCCGTCTTTGTGGGGTTGTTTGGACAAAAAGCGGAAATAAATATgg                                               >  1:1145439/1‑150 (MQ=255)
                                                                                                           ggCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTaa                                            >  1:1060584/1‑150 (MQ=255)
                                                                                                              gTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGAt                                        <  1:589614/151‑1 (MQ=255)
                                                                                                                 gATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACt                                     <  1:1080362/151‑1 (MQ=255)
                                                                                                                  aTGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTATTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACt                                     <  2:39682/150‑1 (MQ=255)
                                                                                                                    gctgctGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCa                                   <  2:179834/150‑1 (MQ=255)
                                                                                                                    gctgctGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCa                                   <  1:601264/150‑1 (MQ=255)
                                                                                                                     ctgctgGCGGGTAAATGGGCGCTTTTAACGGGGGTCTGGCAATGGGCGGGGGGGCGGTTTGTGCCGACGTTTGTTCTGGTGGTGATGCTTATCGGTTTTGGTGCGTTTGTTGGCACAGAAGGCGAAAATAAATATGGTAAAGATACTCCgg                                 >  2:104237/1‑151 (MQ=255)
                                                                                                                       gctgGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGa                                <  2:770984/150‑1 (MQ=255)
                                                                                                                         tgGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgt                             >  2:630976/1‑151 (MQ=255)
                                                                                                                          ggCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGAtgtg                            >  1:90403/1‑151 (MQ=255)
                                                                                                                              gggTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGa                           >  2:1179654/1‑149 (MQ=255)
                                                                                                                              gggTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAg                         <  1:985622/151‑1 (MQ=255)
                                                                                                                                  tAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGtat                      >  1:1233887/1‑151 (MQ=255)
                                                                                                                                    aaCTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGtata                     <  2:685306/151‑1 (MQ=255)
                                                                                                                                       cTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATaa                    >  1:1105015/1‑150 (MQ=255)
                                                                                                                                         tGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATaa                    >  2:125832/1‑149 (MQ=255)
                                                                                                                                         tGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATaa                    >  1:24697/1‑149 (MQ=255)
                                                                                                                                           gggCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAg                  <  2:487459/150‑1 (MQ=255)
                                                                                                                                            ggCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGcc                >  1:318941/1‑151 (MQ=255)
                                                                                                                                                   tttACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACaa           <  2:1231049/150‑1 (MQ=255)
                                                                                                                                                    ttACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAAt          <  1:3951/150‑1 (MQ=255)
                                                                                                                                                        cGGGTGTCTGGCAGTGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAAt      <  1:590173/150‑1 (MQ=255)
                                                                                                                                                           gtgtCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAg  <  2:750457/151‑1 (MQ=255)
                                                                                                                                                           gtgtCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAg  <  1:397437/151‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                            
ACCGCGGCGTTTTGCCTTGTGTGCTTGCTCTGGCCAACAGCGGCAGTAACTGTTAAGCGCTTGCATGATCGCGGGCGTTCCGGCGCATGGGCATTTCTGATGATTGTGGCGTGGATGCTGCTGGCGGGTAACTGGGCGATTTTACCGGGTGTCTGGCAATGGGCGGTGGGGCGTTTTGTCCCGACGTTGATTCTGGTGATGATGCTTATCGATCTTGGTGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAG  >  NC_000913/4112450‑4112748

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