Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,326 A→G 6.1% T24A (ACA→GCA)  yffL → CPZ‑55 prophage; uncharacterized protein YffL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3260AG 6.1% 350.3 / 9.5 115T24A (ACA→GCA) yffLCPZ‑55 prophage; uncharacterized protein YffL
Reads supporting (aligned to +/- strand):  ref base A (51/57);  new base G (3/4);  total (54/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 116 total aligned reads displayed.

GGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGAC  >  NC_000913/2560181‑2560474
                                                                                                                                                 |                                                                                                                                                    
ggCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATg                                                                                                                                                  >  2:1472859/1‑150 (MQ=255)
  cttttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                                                                                >  2:768687/1‑150 (MQ=255)
     ttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                                                             <  1:825339/150‑1 (MQ=255)
     ttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                                                             <  1:1463340/150‑1 (MQ=255)
      tcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                                                                            >  1:246725/1‑150 (MQ=255)
      tcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                                                                            <  1:1432353/150‑1 (MQ=255)
       cttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                                                                           >  2:942053/1‑150 (MQ=255)
       cttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                                                                           <  1:1472859/150‑1 (MQ=255)
         ttttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                                                                         >  1:563195/1‑150 (MQ=255)
          tttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                                                                        <  1:863366/150‑1 (MQ=255)
            tCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                                                                      <  1:597521/150‑1 (MQ=255)
             cGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                                                                                     >  2:12297/1‑150 (MQ=255)
              gCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                                                                    <  2:66443/150‑1 (MQ=255)
              gCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                                                                    >  1:529581/1‑150 (MQ=255)
              gCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                                                                    <  2:678250/150‑1 (MQ=255)
                tAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                                                                  <  2:413682/150‑1 (MQ=255)
                 aaTTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                                                                 <  1:261513/150‑1 (MQ=255)
                  attattccattattccAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                                                                <  1:1207884/150‑1 (MQ=255)
                     attccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                                                                             >  2:365079/1‑150 (MQ=255)
                     attccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                                                                             >  1:1095970/1‑150 (MQ=255)
                        ccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                                                                          <  1:1204362/150‑1 (MQ=255)
                        ccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                                                                          <  2:455194/150‑1 (MQ=255)
                        ccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                                                                          >  1:396673/1‑150 (MQ=255)
                        acattattccaATTAAGTGTAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                                                                          <  2:417841/149‑1 (MQ=255)
                         cattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                                                                         >  1:776271/1‑150 (MQ=255)
                            aattccaATTAACTCCAATAAATTATTCAATATCAACAACATAAATCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                                                                      <  1:323395/149‑1 (MQ=255)
                             attccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                                                                     >  1:103947/1‑150 (MQ=255)
                             attccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                                                                     <  1:350887/150‑1 (MQ=255)
                             attccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                                                                     >  1:392911/1‑150 (MQ=255)
                              ttccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAgg                                                                                                                    <  2:49649/150‑1 (MQ=255)
                              ttccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAgg                                                                                                                    <  2:811701/150‑1 (MQ=255)
                                        gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                                                                          <  2:1285/150‑1 (MQ=255)
                                        gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                                                                          <  2:958754/150‑1 (MQ=255)
                                          tgAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCTGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                                                                        <  1:228220/149‑1 (MQ=255)
                                          ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                                                                        <  2:1075362/150‑1 (MQ=255)
                                           gAATAAATTATTCAATATCAAGAACTTAAACCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATg                                                                                                       <  1:679749/150‑1 (MQ=255)
                                            aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATACAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGt                                                                                                      >  1:1629437/1‑150 (MQ=255)
                                               aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                                                                   <  2:187713/150‑1 (MQ=255)
                                               aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                                                                   <  2:193611/150‑1 (MQ=255)
                                                         atatCAAGAACTTAAATCCATGTTTGCAAGGAAATATCCATAGTTCGAAAATAGGATCAAAACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGATAAACAAGCTAAAGGAGAATTTCTCATTTCTGATATGCTAGc                                                                                         <  2:722071/150‑1 (MQ=255)
                                                            tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                                                                      >  2:1372106/1‑150 (MQ=255)
                                                             cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAACCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                                                                     >  1:1385013/1‑150 (MQ=255)
                                                              aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                                                                                    <  2:246725/150‑1 (MQ=255)
                                                               aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                                                                            <  1:1072677/141‑1 (MQ=255)
                                                                gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                                                                                  >  2:334533/1‑150 (MQ=255)
                                                                   cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                                                                               >  1:1083085/1‑150 (MQ=255)
                                                                     tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                                                                             <  1:365079/150‑1 (MQ=255)
                                                                      aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                                                                            <  2:1303214/150‑1 (MQ=255)
                                                                       aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                                                                           <  2:174828/150‑1 (MQ=255)
                                                                       aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                                                                           <  2:1432445/150‑1 (MQ=255)
                                                                       aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                                                                           >  2:647764/1‑150 (MQ=255)
                                                                        aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                                                                          <  2:634801/150‑1 (MQ=255)
                                                                          ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                                                                        <  2:453408/150‑1 (MQ=255)
                                                                              gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAacaaacaaacaa                                                                    >  1:1577984/1‑150 (MQ=255)
                                                                              gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                                                                    >  1:130999/1‑150 (MQ=255)
                                                                               tttGCCACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                                                                   <  2:1029288/150‑1 (MQ=255)
                                                                                 tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                                                                 <  2:1095970/150‑1 (MQ=255)
                                                                                 tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                                                                 <  1:1015559/150‑1 (MQ=255)
                                                                                   cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg                                                               >  2:1551216/1‑150 (MQ=255)
                                                                                    aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt                                                              <  2:1364898/150‑1 (MQ=255)
                                                                                       gAAAGATCCAGAGTTCGAAAATAGCATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATCTTCg                                                           <  2:854466/150‑1 (MQ=255)
                                                                                        aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGt                                                          >  2:681647/1‑150 (MQ=255)
                                                                                         aaGATCAAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTATGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc                                                         >  2:465496/1‑150 (MQ=255)
                                                                                            aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAt                                                      >  2:407952/1‑150 (MQ=255)
                                                                                              ccAGAGTTCGAAAATAGGATAAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTACAACAAGAAACAACCAATGATGTTCGTCAAttt                                                    >  2:5668/1‑150 (MQ=255)
                                                                                                agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttttt                                                  <  2:396673/150‑1 (MQ=255)
                                                                                                agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGATAAAGGAGAATTTATCATGTCTGATATGATAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttttt                                                  >  2:404175/1‑150 (MQ=255)
                                                                                                 gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGGCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa                                                 <  2:1156796/150‑1 (MQ=255)
                                                                                                 gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa                                                 >  1:123071/1‑150 (MQ=255)
                                                                                                 gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa                                                 <  2:357173/150‑1 (MQ=255)
                                                                                                        aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACaaaa                                          <  2:371568/150‑1 (MQ=255)
                                                                                                            tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTaa                                      >  1:211239/1‑150 (MQ=255)
                                                                                                             aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAAc                                     >  2:1045360/1‑150 (MQ=255)
                                                                                                               gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGa                                   >  2:1068434/1‑150 (MQ=255)
                                                                                                               gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAGTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGa                                   >  1:220620/1‑150 (MQ=255)
                                                                                                                aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGTGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAg                                  <  2:3734/150‑1 (MQ=255)
                                                                                                                 tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGt                                 >  1:58973/1‑150 (MQ=255)
                                                                                                                     cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGAGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACg                             >  2:113058/1‑150 (MQ=255)
                                                                                                                      acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACgg                            <  2:529581/150‑1 (MQ=255)
                                                                                                                       caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGa                           <  1:1043831/150‑1 (MQ=255)
                                                                                                                       caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGa                           >  1:836048/1‑150 (MQ=255)
                                                                                                                       caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGACAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGAAGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGa                           >  1:1126297/1‑150 (MQ=255)
                                                                                                                        aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGATTTACGGaa                          <  2:1281303/150‑1 (MQ=255)
                                                                                                                            aaaaGCCCCCGGAATGATGCAACATGCCGGTGGCGGTATGAGAAACAATCTAAATTAGAATTTCTCATGTCTGATATGATAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATaaa                      <  1:647764/150‑1 (MQ=255)
                                                                                                                             aaaGCCCCCGGAATGCTGCAACATGCCGGGGTCGGTATGAGAAACAAGCTAAATGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAg                     <  1:659294/150‑1 (MQ=255)
                                                                                                                             aaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAg                     >  1:1242784/1‑150 (MQ=255)
                                                                                                                             aaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTACTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAg                     >  2:1063348/1‑150 (MQ=255)
                                                                                                                              ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACg                             <  1:1248561/139‑1 (MQ=255)
                                                                                                                               cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCa                   >  2:1333396/2‑150 (MQ=255)
                                                                                                                               aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCa                   <  2:103947/150‑1 (MQ=255)
                                                                                                                                gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAc                  <  1:830361/150‑1 (MQ=255)
                                                                                                                                  ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAccc                <  2:1083085/150‑1 (MQ=255)
                                                                                                                                  ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAccc                >  2:556711/1‑150 (MQ=255)
                                                                                                                                      ggcaTGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGAACCCaaaa            >  2:1199279/4‑150 (MQ=255)
                                                                                                                                           gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGaa       >  1:30030/1‑150 (MQ=255)
                                                                                                                                           gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGaa       >  2:48072/1‑150 (MQ=255)
                                                                                                                                            ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAag      >  1:1256541/1‑150 (MQ=255)
                                                                                                                                             tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAaga     >  1:801862/1‑150 (MQ=255)
                                                                                                                                               cAACATGCCGCGGGCGGTATGAGAAACAAGCTAAAGCAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAgaga   >  2:532060/1‑150 (MQ=255)
                                                                                                                                                aaCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGAc  <  1:942053/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGAC  >  NC_000913/2560181‑2560474

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: