Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,147,304 G→C 6.3% R56P (CGG→CCG)  frwD → putative PTS enzyme IIB component FrwD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,3040GC 6.3% 518.5 / 15.5 154R56P (CGG→CCG) frwDputative PTS enzyme IIB component FrwD
Reads supporting (aligned to +/- strand):  ref base G (76/68);  new base C (5/5);  total (81/73)
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.93e-01

Only 100 of 163 total aligned reads displayed.

GTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCC  >  NC_000913/4147156‑4147451
                                                                                                                                                    |                                                                                                                                                   
gTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCgg                                                                                                                                                    <  1:1154287/150‑1 (MQ=255)
gTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCgg                                                                                                                                                    >  1:106791/1‑150 (MQ=255)
 tAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGAGGATGAGGATATCCGCCggg                                                                                                                                                   >  2:181213/1‑150 (MQ=255)
    ccGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCgg                                                                                                                                                >  1:202332/1‑150 (MQ=255)
     cGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                               <  1:1165630/150‑1 (MQ=255)
     cGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGAGGAACGGCTGGAAAAGTTGTGCATGTTAGAGAAGTGGGGAGTAAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCAGATGAGGATATACGCCGGGCGGa                                                                                                                                               >  1:766176/1‑150 (MQ=255)
           cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcct                                                                                                                                         <  2:580061/150‑6 (MQ=255)
            gTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGTGGAGTCAGCAGTGAATCTGCGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATCTCCGCCCGGCGGAtatccgc                                                                                                                                        <  2:828311/150‑7 (MQ=255)
                  ggcgTGGCGCATACTTATATGGCGGTGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTg                                                                                                                                  <  2:340985/150‑1 (MQ=255)
                    cgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAt                                                                                                                                <  1:681346/150‑1 (MQ=255)
                      tGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTa                                                                                                                              <  1:1303957/150‑1 (MQ=255)
                         cgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACgg                                                                                                                           <  1:72062/150‑1 (MQ=255)
                          gcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                          <  2:1137936/150‑1 (MQ=255)
                            aTACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCNGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGata                                                                                                                        <  2:59939/150‑1 (MQ=255)
                             tACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                   <  2:1256572/138‑1 (MQ=255)
                               cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCg                                                                                                                     >  1:1145236/1‑150 (MQ=255)
                                   tatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCt                                                                                                                 <  2:106791/150‑1 (MQ=255)
                                   tatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCt                                                                                                                 <  1:484131/150‑1 (MQ=255)
                                     tGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTg                                                                                                               >  2:1179843/1‑150 (MQ=255)
                                      ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                                                              >  2:968612/1‑150 (MQ=255)
                                      ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                                                              >  2:833384/1‑150 (MQ=255)
                                        cggcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCg                                                                                                            <  1:1352071/150‑1 (MQ=255)
                                          gcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGt                                                                                                          >  1:335541/1‑150 (MQ=255)
                                                 ggCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGc                                                                                                   <  1:1032618/150‑1 (MQ=255)
                                                   cTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGa                                                                                                 <  1:1295317/150‑1 (MQ=255)
                                                       aaaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTg                                                                                             >  1:402899/1‑150 (MQ=255)
                                                            ttGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACAt                                                                                        <  1:712295/150‑1 (MQ=255)
                                                               tgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGc                                                                                     >  1:651177/1‑150 (MQ=255)
                                                               tgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGc                                                                                     <  2:951940/150‑1 (MQ=255)
                                                                gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGcc                                                                                    >  1:1111253/1‑150 (MQ=255)
                                                                 ccTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCg                                                                                   <  1:272883/150‑1 (MQ=255)
                                                                   tGTTAGAGAAGTTGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCt                                                                                 <  1:151856/150‑1 (MQ=255)
                                                                   tGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCt                                                                                 >  2:765023/1‑150 (MQ=255)
                                                                    gTTAGAGAAGTGGGTAGTAAGCATTGAAAATCAGGGAGCGATGGGAACGGAGAATCGTTTAGCGGATGAGGATATAAGCAGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTa                                                                                >  2:34825/1‑150 (MQ=255)
                                                                      tAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                          <  1:617002/106‑1 (MQ=255)
                                                                       agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                             >  1:619207/1‑150 (MQ=255)
                                                                        gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgtg                                                                            >  1:308161/1‑150 (MQ=255)
                                                                        gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgtg                                                                            >  1:589411/1‑150 (MQ=255)
                                                                        gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgtg                                                                            >  2:563472/1‑150 (MQ=255)
                                                                         agaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACGTTGCCGCTATGTGc                                                                           <  2:332250/150‑1 (MQ=255)
                                                                         agaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGc                                                                           <  2:651177/150‑1 (MQ=255)
                                                                          gaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCa                                                                          <  2:1026811/150‑1 (MQ=255)
                                                                           aaGTGGGGAGTCATCATTGAAACTCATGGCGATCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                                                                         <  1:124297/150‑1 (MQ=255)
                                                                            aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAAt                                                                        >  1:372521/1‑150 (MQ=255)
                                                                            aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAAt                                                                        >  1:1216974/1‑150 (MQ=255)
                                                                            aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAAt                                                                        <  2:536692/150‑1 (MQ=255)
                                                                               ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                     <  1:181213/150‑1 (MQ=255)
                                                                               ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                     <  1:211479/150‑1 (MQ=255)
                                                                                 ggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagc                                                                   >  1:1413693/1‑150 (MQ=255)
                                                                                 ggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagc                                                                   <  2:301670/150‑1 (MQ=255)
                                                                                    gTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGAGGATGAGGATATCCGCCGGGCCGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATc                                                                >  2:420897/1‑150 (MQ=255)
                                                                                      cAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTa                                                              <  2:619207/150‑1 (MQ=255)
                                                                                        gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                            <  1:515094/150‑1 (MQ=255)
                                                                                         cATTGAAACTCATGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCTAGCTTGCCGGGGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGc                                                           <  1:518159/150‑1 (MQ=255)
                                                                                         cATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGc                                                           >  1:1019392/1‑150 (MQ=255)
                                                                                         cATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGc                                                           >  1:1428960/1‑150 (MQ=255)
                                                                                           ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAt                                                         >  1:961621/1‑150 (MQ=255)
                                                                                           ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAt                                                         <  2:1085621/150‑1 (MQ=255)
                                                                                            tGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAtt                                                        >  1:917428/1‑150 (MQ=255)
                                                                                             gAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc                                                       >  1:815867/1‑150 (MQ=255)
                                                                                             gAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc                                                       >  1:1263182/1‑150 (MQ=255)
                                                                                              aaaCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTcc                                                      <  2:1235923/150‑1 (MQ=255)
                                                                                              aaaCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTcc                                                      <  1:1179843/150‑1 (MQ=255)
                                                                                                   cAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGt                                                 >  1:913947/1‑150 (MQ=255)
                                                                                                    aGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                <  2:385706/150‑1 (MQ=255)
                                                                                                    aGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                <  1:563974/150‑1 (MQ=255)
                                                                                                    aGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                <  1:701732/150‑1 (MQ=255)
                                                                                                       gcgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGc                                             <  2:1413693/150‑1 (MQ=255)
                                                                                                        cgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGcc                                            >  2:88184/1‑150 (MQ=255)
                                                                                                          cgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGc                                          >  2:1071619/1‑150 (MQ=255)
                                                                                                             tGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGc                                       <  1:1006789/150‑1 (MQ=255)
                                                                                                              gggAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCg                                      >  1:1092385/1‑150 (MQ=255)
                                                                                                                gAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCggg                                    >  1:1025954/1‑150 (MQ=255)
                                                                                                                gAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCGCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCggg                                    <  1:989389/150‑1 (MQ=255)
                                                                                                                    ggAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAAt                                >  1:29076/1‑150 (MQ=255)
                                                                                                                       gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAg                             >  1:937813/1‑150 (MQ=255)
                                                                                                                            gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                        <  1:80084/150‑1 (MQ=255)
                                                                                                                            gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                        <  2:1428960/150‑1 (MQ=255)
                                                                                                                            gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                        >  1:306034/1‑150 (MQ=255)
                                                                                                                            gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                        <  2:119051/150‑1 (MQ=255)
                                                                                                                              ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTg                      >  2:675149/1‑150 (MQ=255)
                                                                                                                               tAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                     <  2:92926/150‑1 (MQ=255)
                                                                                                                                 gCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgcgc                   >  2:1288967/1‑150 (MQ=255)
                                                                                                                                   ggATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaa                 >  2:172751/1‑150 (MQ=255)
                                                                                                                                   ggATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaa                 >  1:948652/1‑150 (MQ=255)
                                                                                                                                   ggATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaa                 >  2:1368660/1‑150 (MQ=255)
                                                                                                                                      tGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGt              >  2:465556/1‑150 (MQ=255)
                                                                                                                                        aGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc            >  2:1242824/1‑150 (MQ=255)
                                                                                                                                        aGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc            >  1:54412/1‑150 (MQ=255)
                                                                                                                                         caacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgcg                    >  2:1267992/5‑141 (MQ=255)
                                                                                                                                         caacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCt           >  1:1077857/5‑150 (MQ=255)
                                                                                                                                         caacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGAGCt           >  1:243480/5‑150 (MQ=255)
                                                                                                                                          gATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCtt          <  2:815867/150‑1 (MQ=255)
                                                                                                                                             atCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCt       >  2:618329/1‑150 (MQ=255)
                                                                                                                                              tCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTg      >  1:1232330/1‑150 (MQ=255)
                                                                                                                                              tCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTg      >  1:867722/1‑150 (MQ=255)
                                                                                                                                               ccgccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgc     >  2:445375/1‑150 (MQ=255)
                                                                                                                                               ccgccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgc     >  2:384254/1‑150 (MQ=255)
                                                                                                                                                cgccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgcg    >  1:1280956/1‑150 (MQ=255)
                                                                                                                                                  ccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGcgcc  <  2:503220/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCC  >  NC_000913/4147156‑4147451

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