Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913283,3380GA55.5% 141.3 / inf 476M46I (ATG→ATA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (97/115);  new base A (131/133);  total (228/248)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 481 total aligned reads displayed.

CCGGACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCATGCAGATG  >  NC_000913/283188‑283488
                                                                                                                                                      |                                                                                                                                                      
 cGGACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATg                                                                                                                                                        >  5:1153918/1‑150 (MQ=255)
    aCGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGc                                                                                                                                                    >  4:452384/1‑151 (MQ=255)
     cGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCa                                                                                                                                                   <  3:1536561/151‑1 (MQ=255)
        tAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGa                                                                                                                                                 >  2:1216618/1‑150 (MQ=255)
         aaCGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCAGGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAAc                                                                                                                                               <  4:1776405/151‑1 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGAGGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACGCCGCATATGCTGATGGAATCTCCCAGCGTCAACCTGTTCTGCATAACGCAGAACGcc                                                                                                                                            <  1:1708093/151‑1 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGcc                                                                                                                                            >  8:86262/1‑151 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGcc                                                                                                                                            >  2:1343605/1‑151 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGcc                                                                                                                                            <  5:925134/151‑1 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGcc                                                                                                                                            <  2:759809/151‑1 (MQ=255)
            gATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGcc                                                                                                                                            >  5:433375/1‑151 (MQ=255)
             aTGACCATAGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCg                                                                                                                                           <  2:702829/151‑1 (MQ=255)
               gACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCggg                                                                                                                                         >  3:884698/1‑151 (MQ=255)
                aCCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCggg                                                                                                                                         >  5:1877477/1‑150 (MQ=255)
                 ccATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGAt                                                                                                                                       >  6:1746881/1‑151 (MQ=255)
                 ccATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGAt                                                                                                                                       <  4:1483895/151‑1 (MQ=255)
                    ttGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATggg                                                                                                                                    >  4:1127939/1‑151 (MQ=255)
                    ttGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATggg                                                                                                                                    >  5:1592896/1‑151 (MQ=255)
                     tGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGc                                                                                                                                   >  6:307740/1‑151 (MQ=255)
                      gagaAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCt                                                                                                                                  <  3:1648713/151‑1 (MQ=255)
                       agaAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTg                                                                                                                                 <  4:1052634/151‑1 (MQ=255)
                         aaaaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTggg                                                                                                                               <  8:337980/151‑1 (MQ=255)
                          aaaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGa                                                                                                                              <  7:812594/151‑1 (MQ=255)
                           aaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGAc                                                                                                                             >  3:1399558/1‑151 (MQ=255)
                            aTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACg                                                                                                                            <  6:1609078/151‑1 (MQ=255)
                            aTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACg                                                                                                                            <  2:1611288/151‑1 (MQ=255)
                                tCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCca                                                                                                                         >  5:743339/1‑150 (MQ=255)
                                tCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCca                                                                                                                         >  2:1930981/1‑150 (MQ=255)
                                tCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCGCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCca                                                                                                                         <  5:558431/150‑1 (MQ=255)
                                 cACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCcaac                                                                                                                       >  4:139578/1‑151 (MQ=255)
                                 cACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCcaac                                                                                                                       <  4:980943/151‑1 (MQ=255)
                                  aCCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCcaaca                                                                                                                      <  5:999546/151‑1 (MQ=255)
                                   ccccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGACCGCCGGGATGGGCTGGGACGCcaaca                                                                                                                      >  2:494610/1‑150 (MQ=255)
                                   ccccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCcaaca                                                                                                                      >  1:617862/1‑150 (MQ=255)
                                       gCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTc                                                                                                                 >  2:1141594/1‑151 (MQ=255)
                                         aGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCAc                                                                                                               <  1:700362/151‑1 (MQ=255)
                                                cgcgTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCaaa                                                                                                        <  8:1842561/151‑1 (MQ=255)
                                                 gcgTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAag                                                                                                       <  4:1943257/151‑1 (MQ=255)
                                                   gTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAagag                                                                                                     <  7:596295/151‑1 (MQ=255)
                                                    ttttATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAgg                                                                                                    >  1:312463/1‑151 (MQ=255)
                                                       tATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTg                                                                                                  <  8:1695271/150‑1 (MQ=255)
                                                        aTGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCt                                                                                                <  3:329347/151‑1 (MQ=255)
                                                             gttaTCACTCTCGCGGCGGGGCCTCAGGGCGCTGTGGCGCTGACCCGGCAGATGCTGATGGAATCTCGCAGCGGCAACCTGTTCGGCATGACGCAGATCGCCGGGATGGGCTGGGACTCCAACAAGGTCACCGGCAAAGAGGTGCTGTTTa                                                                                           <  8:541648/148‑1 (MQ=255)
                                                             gTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGAtt                                                                                            >  7:1987930/1‑150 (MQ=255)
                                                                aTCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCg                                                                                        >  7:846744/1‑151 (MQ=255)
                                                                aTCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCg                                                                                        >  8:667441/1‑151 (MQ=255)
                                                                aTCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCg                                                                                        >  6:492286/1‑151 (MQ=255)
                                                                 tCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCgg                                                                                       >  6:746351/1‑151 (MQ=255)
                                                                  cACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGc                                                                                      <  5:1828160/151‑1 (MQ=255)
                                                                   aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCa                                                                                     <  5:856137/151‑1 (MQ=255)
                                                                   aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGc                                                                                      >  5:1732141/1‑150 (MQ=255)
                                                                   aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCa                                                                                     <  6:643066/151‑1 (MQ=255)
                                                                   aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGAGGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCa                                                                                     <  5:1464629/151‑1 (MQ=255)
                                                                      cACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTc                                                                                  >  1:82647/1‑151 (MQ=255)
                                                                      cACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTc                                                                                  >  7:697002/1‑151 (MQ=255)
                                                                       ccgcgGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCa                                                                                 <  7:43297/150‑1 (MQ=255)
                                                                       aCGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTc                                                                                  >  1:234918/1‑150 (MQ=255)
                                                                         gcggcgGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAgg                                                                               >  8:1288111/1‑151 (MQ=255)
                                                                          cggcggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCAGTCAggg                                                                              >  1:1002117/1‑151 (MQ=255)
                                                                            gcggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggcg                                                                            <  3:1769229/151‑1 (MQ=255)
                                                                             cggGGCCGCAGGGGGCTCGGCCGCTGACCGCGCAGAGGCTGAGGGCATCTCCCAGCGGCAACCTCTGCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCCGggcgg                                                                           <  7:1261678/151‑1 (MQ=255)
                                                                                ggCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCAt                                                                        <  4:782409/151‑1 (MQ=255)
                                                                                ggCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCAt                                                                        <  1:431609/151‑1 (MQ=255)
                                                                                  ccGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATcc                                                                      <  1:1138832/151‑1 (MQ=255)
                                                                                    gCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCcgc                                                                    >  4:675679/1‑151 (MQ=255)
                                                                                        ggCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCg                                                                <  3:578249/151‑1 (MQ=255)
                                                                                              ctGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacg                                                          >  2:301830/1‑151 (MQ=255)
                                                                                                     tGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGccc                                                   >  2:1399043/1‑151 (MQ=255)
                                                                                                     tGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGccc                                                   >  6:859543/1‑151 (MQ=255)
                                                                                                         cccGTAGATGGTGATGGAATCACCCAGCGGCAACCTGTTCGGCATAACTCAGAACTCCTGGATGGGCTGGGACTCCAACAAGCTCACCGGCAAAGAGGTGATGATTATCTGAACTCAGGGCGGCATCCGCGCCTGAGACGGACGCCCAATc                                               <  1:133322/151‑1 (MQ=255)
                                                                                                         cccGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATc                                               <  4:164209/151‑1 (MQ=255)
                                                                                                          ccGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATcg                                              <  7:295445/151‑1 (MQ=255)
                                                                                                           cGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATcgc                                             >  8:1231465/1‑151 (MQ=255)
                                                                                                             cAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCG‑‑‑CGGAGACGGACGCCCAATCGCGCTgg                                        <  8:1738018/151‑1 (MQ=255)
                                                                                                               gATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTg                                         >  3:1055542/1‑151 (MQ=255)
                                                                                                               gATGCTGATGGAATCTCCCAGCGGCAACCTGTGCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTg                                         <  7:1683276/151‑1 (MQ=255)
                                                                                                                  gCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTggg                                       <  3:236361/150‑1 (MQ=255)
                                                                                                                      aTGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTAcc                                  >  7:182228/1‑151 (MQ=255)
                                                                                                                      aTGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTAcc                                  <  3:1182475/151‑1 (MQ=255)
                                                                                                                        ggAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACcac                                <  6:1342986/151‑1 (MQ=255)
                                                                                                                        ggAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACcac                                <  4:889642/151‑1 (MQ=255)
                                                                                                                         gAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACcaca                               <  5:1979479/151‑1 (MQ=255)
                                                                                                                         gAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACcac                                >  4:1343121/1‑150 (MQ=255)
                                                                                                                               cccAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCg                            <  6:1131500/148‑1 (MQ=255)
                                                                                                                               cccAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGc                         >  6:995453/1‑151 (MQ=255)
                                                                                                                               cccAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGc                         <  2:895199/151‑1 (MQ=255)
                                                                                                                                  aGCGGCAAGCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt                      <  5:67546/151‑1 (MQ=255)
                                                                                                                                    ggtcAACCTGTTCGGCATAACGCAGAACGCCGGGATGTGCTTGGAAGCCAACAAGCTCACGGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCTGCATCCGCGCCGGAGACGTACGCCCAATCGCGCTTGGCTACCACACCGGGCATTgg                    <  8:1407336/148‑1 (MQ=255)
                                                                                                                                    cGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTgg                    >  6:1263623/1‑151 (MQ=255)
                                                                                                                                         aCCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                   >  6:1849020/1‑147 (MQ=255)
                                                                                                                                          ccTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                   >  6:1581311/1‑146 (MQ=255)
                                                                                                                                          ccTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                   >  2:1002605/1‑146 (MQ=255)
                                                                                                                                          ccTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATc              <  8:1639358/151‑1 (MQ=255)
                                                                                                                                          ccTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATc              <  7:1592281/151‑1 (MQ=255)
                                                                                                                                            tGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                   >  1:1619312/1‑144 (MQ=255)
                                                                                                                                             gTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGc           >  7:1409610/1‑151 (MQ=255)
                                                                                                                                             gTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCg             <  1:518958/149‑1 (MQ=255)
                                                                                                                                              ttCGGCATAACGCTGAACGCCGGGATGGGCTGGGACCCCAACAAGCTCACCCGAAAAGACCTCCTGAGAAGCTGCACTCACGGCGGCATCCGCGCCGGAGACGGACGCGCAAGCGCGCGGGCCTCCCCCACCGGGGATTGGGAGATCGGc           >  3:540367/1‑150 (MQ=255)
                                                                                                                                               tCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCa          <  3:657427/150‑1 (MQ=255)
                                                                                                                                                cGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCatg        >  8:1082894/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGGACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCATGCAGATG  >  NC_000913/283188‑283488

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