Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,264,935 A→T 5.2% N208Y (AAC→TAC)  dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,264,9350AT 5.2% 309.6 / 3.4 94N208Y (AAC→TAC) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (40/49);  new base T (2/3);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACT  >  NC_000913/4264786‑4265083
                                                                                                                                                     |                                                                                                                                                    
tgGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGAAGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACAGTGGCCCG‑ATTTAGCAATTGTTTAAGCAGCCACACGAATGCGTTAACGGGGTtaa                                                                                                                                                     >  2:1363525/1‑150 (MQ=255)
 ggATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTaaa                                                                                                                                                     >  2:1383784/1‑150 (MQ=255)
   aTCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAaca                                                                                                                                                   >  1:968022/1‑150 (MQ=255)
     cTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACAcc                                                                                                                                                 >  2:633545/1‑150 (MQ=255)
      tGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCg                                                                                                                                                >  2:349028/1‑150 (MQ=255)
       ggCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCgg                                                                                                                                               <  1:1453510/150‑1 (MQ=255)
         cTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGtt                                                                                                                                             >  2:26664/1‑150 (MQ=255)
          tGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGTGTTTACCCCGGtaa                                                                                                                                            <  1:1033233/150‑3 (MQ=255)
          tGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTa                                                                                                                                            <  2:1293503/150‑1 (MQ=255)
          tGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTa                                                                                                                                            <  1:1526598/150‑1 (MQ=255)
           gAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTAt                                                                                                                                           <  2:1598144/150‑1 (MQ=255)
            aaTCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATg                                                                                                                                          <  1:163752/150‑1 (MQ=255)
             aTCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGa                                                                                                                                         >  1:776252/1‑150 (MQ=255)
               cccGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGTGTTTACCCCGGtaacgcca                                                                                                                                       <  1:378503/150‑8 (MQ=255)
               cccGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACg                                                                                                                                       <  2:429231/150‑1 (MQ=255)
                  gcgTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGAtc                                                                                                                                    <  2:1152535/150‑1 (MQ=255)
                   cgTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGAtct                                                                                                                                   <  2:43351/150‑1 (MQ=255)
                    gTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGAtctc                                                                                                                                  <  2:335730/150‑1 (MQ=255)
                      cTTTAAATTTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAACAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACGGGTTATGACGATCTcaa                                                                                                                                <  2:516648/150‑1 (MQ=255)
                       tttAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTcaac                                                                                                                               <  1:626569/150‑1 (MQ=255)
                            aaTTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACaaaaa                                                                                                                          <  1:259694/150‑1 (MQ=255)
                                gCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAAcc                                                                                                                      <  1:1374449/150‑1 (MQ=255)
                                gCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAAcc                                                                                                                      >  2:287396/1‑150 (MQ=255)
                                 ccGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCg                                                                                                                     >  2:438268/1‑150 (MQ=255)
                                   gAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCCCACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCt                                                                                                                   >  2:1257546/1‑150 (MQ=255)
                                      aGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGc                                                                                                                <  1:633545/150‑1 (MQ=255)
                                        tCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTTGCGCGTATTGAGCATTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCtt                                                                                                              <  1:624908/150‑1 (MQ=255)
                                            gCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAg                                                                                                          <  1:1385171/150‑1 (MQ=255)
                                                 cAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAAGCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTc                                                                                                     <  1:1456688/150‑1 (MQ=255)
                                                 cAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTc                                                                                                     >  1:707788/1‑150 (MQ=255)
                                                   aaGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCgg                                                                                                   >  2:712291/1‑150 (MQ=255)
                                                     gACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGACGTCGGAt                                                                                                 >  2:95553/1‑150 (MQ=255)
                                                      aCGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGAtt                                                                                                <  2:297033/150‑1 (MQ=255)
                                                      aCGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACAAGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTTGCTTGCAGCCGTCGGAtt                                                                                                <  2:1503886/150‑1 (MQ=255)
                                                         aaGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGa                                                                                             <  1:26664/150‑1 (MQ=255)
                                                         aaGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGa                                                                                             >  2:168922/1‑150 (MQ=255)
                                                          aGGGCCGAAGATCATCGGCGATGTGCTCGACGCAACCGTGGCTCGTATTGAGCAGTTGTTTCAGCATCCACACGATTGCGTTACCGGGGTAAACAGCGGTTATGACGATCTCAACAAAAAAAACGCTGGCTTGCTGCCGTCGGATTTGat                                                                                            <  2:707788/150‑1 (MQ=255)
                                                            ggCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGatca                                                                                          <  2:968022/150‑1 (MQ=255)
                                                            ggCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGatca                                                                                          >  1:203025/1‑150 (MQ=255)
                                                             gCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGatcat                                                                                         >  1:1052774/1‑150 (MQ=255)
                                                             gCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGatcat                                                                                         <  2:522239/150‑1 (MQ=255)
                                                               cGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCAtcg                                                                                       <  1:1318720/150‑1 (MQ=255)
                                                                     aCATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTcgccg                                                                                 <  1:158770/150‑1 (MQ=255)
                                                                        tCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCcgcgc                                                                              >  2:1201369/1‑150 (MQ=255)
                                                                           ccGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGcgtc                                                                           <  1:287396/150‑1 (MQ=255)
                                                                             gATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGcgtccg                                                                         <  2:231283/150‑1 (MQ=255)
                                                                               tgtgCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGcgtccgtc                                                                       >  1:20255/1‑150 (MQ=255)
                                                                                gtgCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCg                                                                      <  2:846011/150‑1 (MQ=255)
                                                                                gtgCGCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCg                                                                      <  2:851590/150‑1 (MQ=255)
                                                                                  gCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGAt                                                                    >  1:1524220/1‑150 (MQ=255)
                                                                                   cTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATg                                                                   <  1:438268/150‑1 (MQ=255)
                                                                                    tCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATgg                                                                  >  1:296188/1‑150 (MQ=255)
                                                                                     cGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATggg                                                                 <  1:1471930/150‑1 (MQ=255)
                                                                                     cGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATggg                                                                 >  1:1320145/1‑150 (MQ=255)
                                                                                       aCGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTa                                                               <  2:864826/150‑1 (MQ=255)
                                                                                        cGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTaa                                                              >  1:789982/1‑150 (MQ=255)
                                                                                         gCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTaaa                                                             <  1:1279274/150‑1 (MQ=255)
                                                                                          cAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTaaaa                                                            >  1:1523898/1‑150 (MQ=255)
                                                                                           aaCCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAaac                                                           <  2:776252/150‑1 (MQ=255)
                                                                                                   cgcgTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTg                                                   >  2:530386/1‑150 (MQ=255)
                                                                                                     cgTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCCTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCg                                                 <  1:897055/150‑1 (MQ=255)
                                                                                                       tATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGAt                                               <  1:712291/150‑1 (MQ=255)
                                                                                                        aTTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATg                                              >  2:405219/1‑150 (MQ=255)
                                                                                                         ttGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGa                                             >  1:1286364/1‑150 (MQ=255)
                                                                                                            aGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGt                                                                                                                                              <  2:930765/50‑1 (MQ=255)
                                                                                                            aGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGt                                                                                                                                              >  1:930765/1‑50 (MQ=255)
                                                                                                            aGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAAcc                                          <  1:168922/150‑1 (MQ=255)
                                                                                                            aGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAAcc                                          <  1:349028/150‑1 (MQ=255)
                                                                                                               aGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCtcg                                       >  1:496638/1‑150 (MQ=255)
                                                                                                                gttgttTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCtcgt                                      >  1:729769/1‑150 (MQ=255)
                                                                                                                gttgttTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCtcgt                                      >  2:1592305/1‑150 (MQ=255)
                                                                                                                   gttTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGa                                   >  1:160620/1‑150 (MQ=255)
                                                                                                                      tCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAacaaca                                                       <  1:1564826/127‑1 (MQ=255)
                                                                                                                      tCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAacaaca                                                       >  2:1564826/1‑127 (MQ=255)
                                                                                                                       cagcagCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAAc                               >  2:649150/1‑150 (MQ=255)
                                                                                                                       caccagcCACACGATGGCGTTACCGGGGTAAACAACGCTTATGACGATCTCAACAAAAAAAAAAAAAGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAAc                               >  2:903531/4‑150 (MQ=255)
                                                                                                                        agcagcCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAAcg                              >  2:1622847/1‑150 (MQ=255)
                                                                                                                        agcagcCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAAcg                              >  1:261553/1‑150 (MQ=255)
                                                                                                                            gcCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACgcggc                          >  1:1562794/1‑150 (MQ=255)
                                                                                                                            gcCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACgcggc                          >  2:687376/1‑150 (MQ=255)
                                                                                                                              cacaCGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGa                        >  2:1422999/1‑150 (MQ=255)
                                                                                                                               acacGATGGCGTTACCGGTGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGAt                       <  2:1181954/150‑1 (MQ=255)
                                                                                                                                 acGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGt                     <  1:514209/150‑1 (MQ=255)
                                                                                                                                    aTGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGc                  <  1:1422999/150‑1 (MQ=255)
                                                                                                                                     cgtcataaCCGGTGTTTACCCCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCa                 >  2:1271476/8‑150 (MQ=255)
                                                                                                                                      ggCGGTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCGTGCAGCCGTCGGATTTGATCATCGGCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATTTTGCAg                <  2:1621767/150‑1 (MQ=255)
                                                                                                                                        cataaCCGGTGTTTACCCCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGa              >  2:326205/5‑150 (MQ=255)
                                                                                                                                        cGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGa              >  2:1590064/1‑150 (MQ=255)
                                                                                                                                         gTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGAt             <  1:1201369/150‑1 (MQ=255)
                                                                                                                                         ataaCCGGTGTTTACCCCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGAt             <  1:326205/147‑1 (MQ=255)
                                                                                                                                                 ggTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGt     <  2:160620/150‑1 (MQ=255)
                                                                                                                                                  gTAAACACCGGTTATGACGATCTCAACAAAAAACCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACCACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGCTa    >  2:290968/1‑150 (MQ=255)
                                                                                                                                                  gTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTa    <  1:482798/150‑1 (MQ=255)
                                                                                                                                                  gTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTa    <  1:1592305/150‑1 (MQ=255)
                                                                                                                                                  gTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGACGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTa    >  1:323743/1‑150 (MQ=255)
                                                                                                                                                   tAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTAc   <  2:203025/150‑1 (MQ=255)
                                                                                                                                                    aaaCACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACt  <  2:1562794/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACT  >  NC_000913/4264786‑4265083

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: ‑