Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,049 A→T 7.6% intergenic (+17/+79) trpT → / ← hdfR tRNA‑Trp/DNA‑binding transcriptional dual regulator HdfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0490AT 7.6% 509.8 / 22.5 166intergenic (+17/+79)trpT/hdfRtRNA‑Trp/DNA‑binding transcriptional dual regulator HdfR
Reads supporting (aligned to +/- strand):  ref base A (79/74);  new base T (7/6);  total (86/80)
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

Only 100 of 170 total aligned reads displayed.

TGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCA  >  NC_000913/3946901‑3947198
                                                                                                                                                      |                                                                                                                                                     
 gcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaaa                                                                                                                                                     >  1:1536337/1‑150 (MQ=255)
  cctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAg                                                                                                                                                    <  1:399845/150‑1 (MQ=255)
   ctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGGTCCGCCACCCTAATTAGG‑‑GGCGTAGTGCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGc                                                                                                                                                   <  2:547976/150‑1 (MQ=255)
       cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                               >  2:1777808/1‑150 (MQ=255)
       cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                               <  2:1119082/150‑1 (MQ=255)
          gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAt                                                                                                                                            <  2:1526500/150‑1 (MQ=255)
           cAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAtg                                                                                                                                           <  1:647340/150‑2 (MQ=255)
            aGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAtga                                                                                                                                          <  2:174886/150‑3 (MQ=255)
               gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCTGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                       <  1:1229207/150‑1 (MQ=255)
                 gTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCAACGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttttt                                                                                                                                     >  1:2153498/1‑150 (MQ=255)
                  tCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                    >  2:496068/1‑150 (MQ=255)
                  gcgcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                    <  1:433172/149‑1 (MQ=255)
                   cgcgGGTTCGAGGCCCGTCCGTTCCGCGACCCGAATTAGG‑‑GGGGGAGTTCAGGTGGTAGAGCACCGGTCTCCAAAGCCGGGTGTTGGGAGTTCGGGTCTCTCCGCCCCGGCCAGAAATCATCCTTAGCGATAGCTAAGGAGTTTTGTTAt                                                                                                                                   >  2:1467027/1‑150 (MQ=255)
                      gggTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTg                                                                                                                                <  2:1073619/150‑1 (MQ=255)
                       ggTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTGCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGa                                                                                                                               <  2:1005822/150‑1 (MQ=255)
                         ttCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaaa                                                                                                                             >  1:757762/1‑150 (MQ=255)
                          tCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAAt                                                                                                                            <  1:66121/150‑1 (MQ=255)
                          tCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAAt                                                                                                                            <  2:2458385/150‑1 (MQ=255)
                            gAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTAGAGTCTCTCCGCCCCTGCCACAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATATGAAATaa                                                                                                                          >  2:586234/1‑150 (MQ=255)
                              gTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAAcc                                                                                                                        <  1:2039824/150‑1 (MQ=255)
                                 ccgtccgtTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCctc                                                                                                                     >  1:1083979/1‑150 (MQ=255)
                                   gtccgtTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCtctc                                                                                                                   <  1:898352/150‑1 (MQ=255)
                                      cgtTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGa                                                                                                                >  2:2390673/1‑150 (MQ=255)
                                       gtTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGaa                                                                                                               <  1:2411685/150‑1 (MQ=255)
                                       gtTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGaa                                                                                                               >  2:1534699/1‑150 (MQ=255)
                                          ccGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                            <  2:688934/150‑1 (MQ=255)
                                           cGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTaa                                                                                                           >  1:531872/1‑150 (MQ=255)
                                                cccTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                      >  1:730591/1‑150 (MQ=255)
                                                cccTAATTAGG‑‑GGCGTAGTTCAATTGGTAAAGCACCGGTCTCCAAACCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                      >  1:1569774/1‑150 (MQ=255)
                                                  cTAATTAGGGTGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                      <  1:518025/150‑1 (MQ=255)
                                                  cTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                                                                                                    >  1:1262392/1‑150 (MQ=255)
                                                  cTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCATCTCCGAAGTAAATCCtta                                                                                                    >  2:1700837/1‑149 (MQ=255)
                                                          g‑‑ggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCa                                                                                            <  2:364155/150‑1 (MQ=255)
                                                          g‑‑ggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCa                                                                                            >  2:2250265/1‑150 (MQ=255)
                                                                  aGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTg                                                                                      >  2:228191/1‑150 (MQ=255)
                                                                   gTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgc                                                                                     >  1:1740953/1‑150 (MQ=255)
                                                                    ttCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcc                                                                                    >  2:2174309/1‑150 (MQ=255)
                                                                       aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTCTCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 <  2:1411632/150‑1 (MQ=255)
                                                                       aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 >  1:2070541/1‑150 (MQ=255)
                                                                             tAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                           <  1:995490/150‑1 (MQ=255)
                                                                             tAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                           <  1:610374/150‑1 (MQ=255)
                                                                             tAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                           <  2:1555001/150‑1 (MQ=255)
                                                                                 gCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATAtt                                                                       <  1:1663126/150‑1 (MQ=255)
                                                                                 gCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGACATTCATAtt                                                                       <  1:1496203/150‑1 (MQ=255)
                                                                                   aCCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTTCCAGCCATTCATATTaa                                                                     <  2:2365290/150‑1 (MQ=255)
                                                                                    ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                    >  1:337987/1‑150 (MQ=255)
                                                                                    ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                    <  1:1056388/150‑1 (MQ=255)
                                                                                    ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                    >  1:1980801/1‑150 (MQ=255)
                                                                                    ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                    <  1:2393291/150‑1 (MQ=255)
                                                                                    ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                    >  2:588435/1‑150 (MQ=255)
                                                                                      ggTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATac                                                                  <  1:228191/150‑1 (MQ=255)
                                                                                      ggTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCATTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATACCCCTCTCCGAAGTAAATCCTTCTACCGGAATCCGTGCCAGCCATTCATATTAATac                                                                  >  2:1198444/1‑150 (MQ=255)
                                                                                         ctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACt                                                               <  1:2060413/150‑1 (MQ=255)
                                                                                         ctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACt                                                               <  1:1146127/150‑1 (MQ=255)
                                                                                         ctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACt                                                               <  1:934706/150‑1 (MQ=255)
                                                                                           ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                             >  2:598886/1‑150 (MQ=255)
                                                                                            cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCa                                                            >  2:1976194/1‑150 (MQ=255)
                                                                                            cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCa                                                            <  1:784087/150‑1 (MQ=255)
                                                                                                  cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGc                                                      >  1:1825661/1‑150 (MQ=255)
                                                                                                  cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGc                                                      <  1:1431439/150‑1 (MQ=255)
                                                                                                  cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGATAAAGATAAGGATTTTTTTTATCTGAAATAACCATCTACGAACTAAATCCTTCTACCGGCATCCTTGCCAGCGATTCATATTAATACACTTCATCCAGc                                                      >  1:538381/1‑150 (MQ=255)
                                                                                                   gggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCa                                                     <  1:1707671/150‑1 (MQ=255)
                                                                                                    gttgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                    <  2:869379/148‑1 (MQ=255)
                                                                                                    ggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                    <  1:2368028/150‑1 (MQ=255)
                                                                                                       gtTGGCAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTATATCTCAAATAACCCTCTCCGAAGTAAATACTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGtt                                                 >  2:1114080/1‑150 (MQ=255)
                                                                                                            gAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAAttt                                            <  2:202991/150‑1 (MQ=255)
                                                                                                            gAGTTCGAGTCTCTACGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAAttt                                            <  1:2321093/150‑1 (MQ=255)
                                                                                                              gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTc                                          >  1:1803159/1‑150 (MQ=255)
                                                                                                                tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaa                                        <  2:754662/150‑1 (MQ=255)
                                                                                                                 cGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaa                                       >  1:1078513/1‑150 (MQ=255)
                                                                                                                  gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                      <  2:2470626/150‑1 (MQ=255)
                                                                                                                   aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAg                                     >  2:2445162/1‑150 (MQ=255)
                                                                                                                     tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                   >  2:274973/1‑150 (MQ=255)
                                                                                                                      ctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATc                                  >  2:325022/1‑150 (MQ=255)
                                                                                                                       tctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcg                                 >  2:1919633/1‑150 (MQ=255)
                                                                                                                          ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                              >  2:1282431/1‑150 (MQ=255)
                                                                                                                           cGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGaa                             >  2:1528272/1‑150 (MQ=255)
                                                                                                                            gCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAAt                            >  2:1903335/1‑150 (MQ=255)
                                                                                                                            gCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAAt                            <  2:1512387/150‑1 (MQ=255)
                                                                                                                              cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCa                          <  1:177612/150‑1 (MQ=255)
                                                                                                                              cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCa                          >  1:461393/1‑150 (MQ=255)
                                                                                                                                cTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAAc                        <  2:932109/150‑1 (MQ=255)
                                                                                                                                  gCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGc                      >  1:252365/1‑150 (MQ=255)
                                                                                                                                  gCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGc                      >  1:2227108/1‑150 (MQ=255)
                                                                                                                                      gAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                  <  1:1111712/150‑1 (MQ=255)
                                                                                                                                          tcatcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGACATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATATTTaa              >  2:1554603/1‑149 (MQ=255)
                                                                                                                                           catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCTAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             <  1:2150660/150‑1 (MQ=255)
                                                                                                                                           catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             >  2:1324290/1‑150 (MQ=255)
                                                                                                                                           catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             <  1:2100906/150‑1 (MQ=255)
                                                                                                                                           catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             <  1:2469318/150‑1 (MQ=255)
                                                                                                                                           aatcCTTAGCTTTCGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAATTCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             >  1:273555/2‑150 (MQ=255)
                                                                                                                                           aatcCTTAGCTTTCGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc             >  2:1717186/2‑150 (MQ=255)
                                                                                                                                            atcCTTAGCTTTCGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCg            >  1:1320801/1‑150 (MQ=255)
                                                                                                                                            atcCTTAGCTTTCGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCg            >  2:1388013/1‑150 (MQ=255)
                                                                                                                                              ccTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCt          >  2:2063739/1‑150 (MQ=255)
                                                                                                                                               cTTAGCGAACTATAAGGATTATATTTATCTGAATTATCCCTCTCCGAACTAAATCCTTCTACCGGCATCCTTGCCAGCCATTAATATTAATAGACTTCATCGAGCACGTTAATTTTCAAAAGATCCCGAATCTACGCATTTTTATCGCTa         <  2:428771/150‑1 (MQ=255)
                                                                                                                                                 tAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAtt       <  2:1575739/150‑1 (MQ=255)
                                                                                                                                                   gCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAtttt     >  2:1253863/1‑150 (MQ=255)
                                                                                                                                                    cGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg    >  1:2256968/1‑150 (MQ=255)
                                                                                                                                                      aaaGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTCATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCa  <  2:502000/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGG‑‑GGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCA  >  NC_000913/3946901‑3947198

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: