Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,264,931 G→T 5.4% G206G (GGG→GGT)  dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,264,9310GT 5.4% 301.0 / 2.7 91G206G (GGG→GGT) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base T (2/4);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTGGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGT  >  NC_000913/4264782‑4265080
                                                                                                                                                     |                                                                                                                                                     
ctgctgGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCggg                                                                                                                                                       >  1:451638/1‑150 (MQ=255)
ctgctgGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCggg                                                                                                                                                       >  2:1471930/1‑150 (MQ=255)
  gctgGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGt                                                                                                                                                     >  1:1185701/1‑150 (MQ=255)
    tgGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGAAGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACAGTGGCCCG‑ATTTAGCAATTGTTTAAGCAGCCACACGAATGCGTTAACGGGGTtaa                                                                                                                                                  >  2:1363525/1‑150 (MQ=255)
     ggATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTaaa                                                                                                                                                  >  2:1383784/1‑150 (MQ=255)
       aTCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAaca                                                                                                                                                >  1:968022/1‑150 (MQ=255)
         cTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGgtgtttacccc                                                                                                                                              <  1:974383/150‑8 (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:1320145/1‑150 (MQ=255)
                                                                                         cGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATggg                                                              <  1:1471930/150‑1 (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                                                                                                                                           >  1:930765/1‑50 (MQ=255)
                                                                                                                aGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGt                                                                                                                                           <  2:930765/50‑1 (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                                                    >  2:1564826/1‑127 (MQ=255)
                                                                                                                          tCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAacaaca                                                    <  1:1564826/127‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
CTGCTGGATCTGGCTGAATCCCGCGTCTTTAAAATTGCCGAAAGTCGTGCGAACAAAGACGAAGGGCCGAAGAACATCGCCGATGTGCTCGACGCAACCGTGGCGCGTATTGAGCAGTTGTTTCAGCAGCCACACGATGGCGTTACCGGGGTAAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGT  >  NC_000913/4264782‑4265080

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