Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,418 C→A 6.9% intergenic (+18/+226) intA → / ← yfjH CP4‑57 prophage; integrase/CP4‑57 prophage; protein YfjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,4180CA 6.9% 356.1 / 9.5 108intergenic (+18/+226)intA/yfjHCP4‑57 prophage; integrase/CP4‑57 prophage; protein YfjH
Reads supporting (aligned to +/- strand):  ref base C (48/51);  new base A (5/3);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757269‑2757567
                                                                                                                                                     |                                                                                                                                                     
cgcgCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGc                                                                                                                                                       >  2:278838/1‑150 (MQ=255)
   gcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctct                                                                                                                                                    >  2:1738877/1‑150 (MQ=255)
   gcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctct                                                                                                                                                    >  2:665394/1‑150 (MQ=255)
     ttATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctctct                                                                                                                                                  <  1:1527647/150‑1 (MQ=255)
      tataACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTg                                                                                                                                                 >  2:1719374/1‑150 (MQ=255)
      tataACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTg                                                                                                                                                 >  1:1550854/1‑150 (MQ=255)
        taACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTggg                                                                                                                                               >  1:1579769/1‑150 (MQ=255)
        taACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTggg                                                                                                                                               >  1:1126410/1‑150 (MQ=255)
         aaCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGc                                                                                                                                              <  2:175185/150‑1 (MQ=255)
          aCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCt                                                                                                                                             <  1:1071717/150‑1 (MQ=255)
          aCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCt                                                                                                                                             >  1:490714/1‑150 (MQ=255)
           ccGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTg                                                                                                                                            >  2:1977697/1‑150 (MQ=255)
             gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAAAGCTCGGCTGATCTCGTAAAAGCAGAAAATACTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGat                                                                                                                                          >  2:1268004/1‑150 (MQ=255)
              cAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGata                                                                                                                                         <  1:499161/150‑1 (MQ=255)
                gTGATTATCTTGAGCAACGTCGTCCGATGATGAAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATTTTTATCAGCCCAGCTCTCTGGGCTGATAtt                                                                                                                                       >  2:836269/1‑150 (MQ=255)
                gTGATTATCTTGAGCAACGTCGGCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGataaa                                                                                                                                       <  1:1038676/150‑3 (MQ=255)
                 tGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGataaat                                                                                                                                      <  2:302830/150‑4 (MQ=255)
                  gATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTc                                                                                                                                     >  1:1893109/1‑150 (MQ=255)
                   aTTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCt                                                                                                                                    <  2:1569201/150‑1 (MQ=255)
                    ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   >  2:415173/1‑150 (MQ=255)
                       tCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaaaa                                                                                                                                <  2:1550854/150‑1 (MQ=255)
                        cTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAg                                                                                                                               <  1:154561/150‑1 (MQ=255)
                         ttGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGa                                                                                                                              <  2:392682/150‑1 (MQ=255)
                           gAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATc                                                                                                                            >  1:1729640/1‑150 (MQ=255)
                             gCAACGTCGTCCGATGATGCAATGGTGGTCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          <  2:1188821/150‑1 (MQ=255)
                               aaCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        <  1:965111/150‑1 (MQ=255)
                                 cgtcgtCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTa                                                                                                                      >  2:810126/1‑150 (MQ=255)
                                   tcgtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCTTCTTTTCGGATAAAATATTTATGAGCCCAGCTCTCTGGGATGATATTTCTAAAAGATCGGAATaat                                                                                                                    <  2:19381/150‑2 (MQ=255)
                                     gtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATa                                                                                                                  <  1:1526035/150‑1 (MQ=255)
                                       ccGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGc                                                                                                                >  2:389967/1‑150 (MQ=255)
                                             atgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcgcg                                                                                                          >  1:876579/1‑150 (MQ=255)
                                             atgCAAGGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTACATAGCTAcgcg                                                                                                          <  2:1638151/150‑1 (MQ=255)
                                              tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                         >  2:384305/1‑150 (MQ=255)
                                               tcAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGa                                                                                                        <  1:1738877/149‑1 (MQ=255)
                                               gCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGa                                                                                                        >  1:973089/1‑150 (MQ=255)
                                                 aaTGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAAt                                                                                                      <  2:911208/150‑1 (MQ=255)
                                                     gtggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTg                                                                                                  <  1:1719374/150‑1 (MQ=255)
                                                       gggCTGATCTCGTAAAAGCAGCAGATAGCGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGct                                                                                                >  2:704089/1‑150 (MQ=255)
                                                        ggCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCTCACTAAATAGCTACGCGGAATAGTGctc                                                                                               <  2:1893109/150‑1 (MQ=255)
                                                          cTGATCTCGTAAAAGCAGCAGATAGTTGTAGAATAGTTTTAACTCATTTCAGCAAAATTCTACTTATAAGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTa                                                                                             >  2:654556/1‑150 (MQ=255)
                                                          cTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTa                                                                                             >  2:449157/1‑150 (MQ=255)
                                                               ctcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTc                                                                                        >  2:1359324/1‑150 (MQ=255)
                                                                tcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCg                                                                                       <  1:433644/150‑1 (MQ=255)
                                                                 cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                      >  1:586720/1‑150 (MQ=255)
                                                                    aaaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTa                                                                                   <  1:679712/150‑1 (MQ=255)
                                                                     aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  <  2:1741430/150‑1 (MQ=255)
                                                                      aaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGc                                                                                 <  1:389967/150‑1 (MQ=255)
                                                                            cagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTc                                                                           <  2:291374/150‑1 (MQ=255)
                                                                               aTAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCaga                                                                        <  1:384305/150‑1 (MQ=255)
                                                                                  gTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGt                                                                     <  1:1977697/150‑1 (MQ=255)
                                                                                  gTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGt                                                                     <  2:1917673/150‑1 (MQ=255)
                                                                                  gTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGt                                                                     <  2:398880/150‑1 (MQ=255)
                                                                                   tGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTa                                                                    >  1:945085/1‑150 (MQ=255)
                                                                                   tGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTa                                                                    <  2:3896/150‑1 (MQ=255)
                                                                                    ggTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTaa                                                                   <  1:665394/150‑1 (MQ=255)
                                                                                    ggTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTaa                                                                   >  2:1494954/1‑150 (MQ=255)
                                                                                      tAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGa                                                                 >  2:1406266/1‑150 (MQ=255)
                                                                                       ggCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        >  2:288437/2‑94 (MQ=255)
                                                                                       cgCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        <  1:288437/93‑1 (MQ=255)
                                                                                         cATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACAt                                                              >  1:834239/1‑150 (MQ=255)
                                                                                             gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          <  2:171999/150‑1 (MQ=255)
                                                                                             gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          >  1:1050107/1‑150 (MQ=255)
                                                                                              ttttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAggg                                                         <  1:460721/150‑1 (MQ=255)
                                                                                              ttttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAggg                                                         <  2:501023/150‑1 (MQ=255)
                                                                                                ttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCa                                                       >  2:560877/1‑150 (MQ=255)
                                                                                                  aaCTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAAt                                                     >  2:487640/1‑150 (MQ=255)
                                                                                                       aTTTGAGCAAAATTCGTCTTGTCGTATAAAATATTTATCAGCCAAGATCTCTGGGCTGATATTTATAAAAGATAGCACTAAATAGCTACGCGGAATAGTGCTCTAATTTCCTTTAGCCATTTCAGAGGTAAGACATTAGGGCAATTTAAc                                                >  2:1156156/1‑150 (MQ=255)
                                                                                                       aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                >  2:1238463/1‑150 (MQ=255)
                                                                                                       aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                >  2:1542700/1‑150 (MQ=255)
                                                                                                        tttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACt                                               <  1:449157/150‑1 (MQ=255)
                                                                                                        tttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACt                                               <  1:704089/150‑1 (MQ=255)
                                                                                                         ttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                              <  1:1936620/150‑1 (MQ=255)
                                                                                                           gAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTa                                            >  2:620575/1‑150 (MQ=255)
                                                                                                             ccAAAATTCGTCTTGTAGTATAAAATATTTATCAGCCCAGGTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:1503206/149‑1 (MQ=255)
                                                                                                                 aaTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTg                                      <  2:1963245/150‑1 (MQ=255)
                                                                                                                      gTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAg                                 <  2:1579769/150‑1 (MQ=255)
                                                                                                                      gTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAg                                 <  1:1414526/150‑1 (MQ=255)
                                                                                                                           gTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAta                            >  1:477985/1‑149 (MQ=255)
                                                                                                                           gTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATc                            >  2:1354432/1‑150 (MQ=255)
                                                                                                                            tCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCg                           <  2:1470243/150‑1 (MQ=255)
                                                                                                                             cGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCCCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGa                          >  2:1171758/1‑150 (MQ=255)
                                                                                                                              ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAg                         >  2:702738/1‑150 (MQ=255)
                                                                                                                              ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAg                         >  1:225584/1‑150 (MQ=255)
                                                                                                                                  aaaaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtg                     >  2:1598710/1‑150 (MQ=255)
                                                                                                                                  aaaaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtg                     >  2:752579/1‑150 (MQ=255)
                                                                                                                                    tagaaatATCAGCCCAGAGAGCTGGGCTGATATTTCTAATAGATCGCACTAAATAGCGACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAGGCCAGTAGGGCATTTTCACTGTATGACTGCGCAGTCATCGAGTCtgtg                   >  1:1652231/7‑150 (MQ=255)
                                                                                                                                    aaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtgtg                   <  1:797081/150‑1 (MQ=255)
                                                                                                                                    aaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtgtg                   <  1:1494954/150‑1 (MQ=255)
                                                                                                                                     atatTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgg                  >  2:648727/1‑150 (MQ=255)
                                                                                                                                      gaaatATCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTggg                 >  1:1497532/5‑150 (MQ=255)
                                                                                                                                       atTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgggg                <  1:1268004/150‑1 (MQ=255)
                                                                                                                                          tATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGgtgt             <  1:1354432/150‑1 (MQ=255)
                                                                                                                                           aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            <  2:876579/150‑1 (MQ=255)
                                                                                                                                           aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            >  1:1427748/1‑150 (MQ=255)
                                                                                                                                            tCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTg           >  1:41116/1‑150 (MQ=255)
                                                                                                                                               gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTc        >  2:962946/1‑150 (MQ=255)
                                                                                                                                                cccAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTcc       >  1:83614/1‑150 (MQ=255)
                                                                                                                                                 ccAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCt      >  2:825905/1‑150 (MQ=255)
                                                                                                                                                  cAGCTCTCTGGGCTTATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCtt     <  2:1050107/150‑1 (MQ=255)
                                                                                                                                                     ctctctGGGCTGATATTTCTAAAAGATCGCACTAAATCGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg  <  1:752579/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757269‑2757567

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