Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,882,514 G→T 5.7% P406T (CCG→ACG)  dnaA ← chromosomal replication initiator protein DnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,882,5140GT 5.7% 412.9 / 9.9 122P406T (CCG→ACG) dnaAchromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base G (61/54);  new base T (5/2);  total (66/56)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

Only 100 of 124 total aligned reads displayed.

TTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAGT  >  NC_000913/3882365‑3882660
                                                                                                                                                     |                                                                                                                                                  
tttGATATCGTGGCTCTCTTCACGCAACTGCTCTATCTTACGGCAGTCATGAAGCACCGGCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCgg                                                                                                                                                    <  1:35349/150‑1 (MQ=255)
tttGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCgg                                                                                                                                                    >  1:353441/1‑150 (MQ=255)
       tCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGcc                                                                                                                                             >  1:1358267/1‑150 (MQ=255)
        cGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCa                                                                                                                                            >  2:936744/1‑150 (MQ=255)
            gCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGa                                                                                                                                        <  1:524205/150‑1 (MQ=255)
             ctctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagatgg                                                                                                                                       <  1:407349/150‑6 (MQ=255)
                 ctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCggg                                                                                                                                   >  2:616039/1‑150 (MQ=255)
                  ttCCCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGa                                                                                                                                  >  2:1529654/1‑150 (MQ=255)
                    cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGCCTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCCCGGACGCGCCACCGCGCGGGATc                                                                                                                                >  2:462643/1‑150 (MQ=255)
                    cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                                <  1:1007155/150‑1 (MQ=255)
                    cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                                >  2:969633/1‑150 (MQ=255)
                    cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                                >  1:1124672/1‑150 (MQ=255)
                     aCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCg                                                                                                                               >  2:410844/1‑150 (MQ=255)
                         aaCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGc                                                                                                                           <  2:1124672/150‑1 (MQ=255)
                          aCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGAcg                                                                                                                                                 <  1:176548/127‑1 (MQ=255)
                          aCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCt                                                                                                                          <  2:1357257/150‑1 (MQ=255)
                           cTGCTCGATCTTACGGCAGGCATGAACCACCGTCGTGTGGTCACGGCCACGAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCTGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCtt                                                                                                                         <  1:711194/150‑1 (MQ=255)
                            tGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTg                                                                                                                        <  2:741880/150‑1 (MQ=255)
                             gCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTgg                                                                                                                       >  2:1132764/1‑150 (MQ=255)
                              cTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                      <  1:128341/150‑1 (MQ=255)
                               tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATTGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaa                                                                                                                     >  1:1172220/1‑150 (MQ=255)
                                cGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                    >  1:886874/1‑150 (MQ=255)
                                 gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagatgg                                                                                                                                       <  1:1283977/130‑6 (MQ=255)
                                 gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagatgg                                                                                                                                       >  2:1283977/1‑125 (MQ=255)
                                 gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTg                                                                                                                        <  2:812968/145‑1 (MQ=255)
                                 gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTg                                                                                                                        >  1:812968/1‑145 (MQ=255)
                                  aTCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCTACGCTTGGAAAgg                                                                                                                  <  2:343056/150‑1 (MQ=255)
                                   tCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGga                                                                                                                 <  1:991016/150‑1 (MQ=255)
                                    cTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGgag                                                                                                                <  1:470832/150‑1 (MQ=255)
                                         ggcaggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCGTGGAAAGGAGATCcg                                                                                                           >  1:920805/1‑150 (MQ=255)
                                             ggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGAc                                                                                                       <  2:1127365/150‑1 (MQ=255)
                                                  gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGa                                                                                                  <  1:1371819/150‑1 (MQ=255)
                                                  gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGa                                                                                                  >  2:769424/1‑150 (MQ=255)
                                                   aaGCACCGTCGTGTGGTCACCGCCACCAAACGCATCGCCAAACTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCAACGCCATCTGGCCCGGACGCGCCACCCAGCGGGATCCACGCTTGGAAAGCAGATCCGCGACTTTGaa                                                                                                 >  2:1248991/1‑149 (MQ=255)
                                                     gCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCt                                                                                               >  2:354711/1‑150 (MQ=255)
                                                      cACCGTCGTGTGGCCACGGCCACCAAACGCATCGCCAATAACCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCCCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCtt                                                                                              >  2:912520/1‑150 (MQ=255)
                                                        ccGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                                                            <  2:816749/150‑1 (MQ=255)
                                                         cgtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                           <  2:849366/150‑1 (MQ=255)
                                                         cgtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                           >  1:704011/1‑150 (MQ=255)
                                                             gtgtGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTAc                                                                                       >  2:1058777/1‑150 (MQ=255)
                                                              tgtgGTCACGGCCACAAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTGGGAAAGGAGATCCGCGACTTTGATCTTgtagtaat                                                                                      >  2:193370/1‑148 (MQ=255)
                                                                tgGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTcc                                                                                    >  1:1529761/1‑150 (MQ=255)
                                                                 ggTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCg                                                                                   <  2:183331/150‑1 (MQ=255)
                                                                 ggTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCg                                                                                   <  2:704011/150‑1 (MQ=255)
                                                                  gTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCGTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                                  >  1:1344907/1‑150 (MQ=255)
                                                                   tCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGcc                                                                                 >  1:1031093/1‑150 (MQ=255)
                                                                         ccaccaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtc                                                                           >  1:459291/1‑150 (MQ=255)
                                                                             caAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                                          <  1:1044970/147‑1 (MQ=255)
                                                                             caAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtct                                                                          >  2:1044970/1‑147 (MQ=255)
                                                                              aaaCGCATCGCCTATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCGGCCACCGtcttct                                                                       <  2:637919/149‑1 (MQ=255)
                                                                              aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttct                                                                       >  1:637919/1‑149 (MQ=255)
                                                                              aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTg                                                                      >  2:1222199/1‑150 (MQ=255)
                                                                              aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTg                                                                      <  1:969633/150‑1 (MQ=255)
                                                                                 cGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAat                                                                   <  2:44853/150‑1 (MQ=255)
                                                                                       gCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                             >  1:101269/1‑150 (MQ=255)
                                                                                       gCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                             >  1:1395333/1‑150 (MQ=255)
                                                                                        ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGgag                                                                                                                >  2:153669/1‑98 (MQ=255)
                                                                                        ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCg                                                            >  2:174765/1‑150 (MQ=255)
                                                                                        ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCg                                                            >  2:1616240/1‑150 (MQ=255)
                                                                                            tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                                        <  1:376833/150‑1 (MQ=255)
                                                                                            tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                                        <  1:504223/150‑1 (MQ=255)
                                                                                               ccGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa                                                     >  2:1298614/1‑150 (MQ=255)
                                                                                                 ggCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTg                                                      <  1:547011/147‑1 (MQ=255)
                                                                                                  gCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCa                                                  >  1:1370621/1‑150 (MQ=255)
                                                                                                         gtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAGGGTGGCCAGTTTTTc                                           >  2:438023/1‑150 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTGGCCAGCGCCATCGCCTTCTGGCGCGGACGCGCCACCGAGCGGGATGGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          <  1:177526/150‑1 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          <  2:1358267/150‑1 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          <  1:1298614/150‑1 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          <  1:1132764/150‑1 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          <  2:459291/150‑1 (MQ=255)
                                                                                                          tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTcc                                          >  2:962217/1‑150 (MQ=255)
                                                                                                           ggTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                                         >  1:318611/1‑150 (MQ=255)
                                                                                                             ttAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGc                                       <  1:354257/150‑1 (MQ=255)
                                                                                                             gtAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGc                                       <  2:652678/149‑1 (MQ=255)
                                                                                                                gTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAAt                                    >  2:987831/1‑150 (MQ=255)
                                                                                                                 tCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                                   <  2:353441/150‑1 (MQ=255)
                                                                                                                 tCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATg                                   <  1:912520/150‑1 (MQ=255)
                                                                                                                  cAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGc                                  >  2:180958/1‑150 (MQ=255)
                                                                                                                  aaGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGc                                  >  2:1238917/2‑150 (MQ=255)
                                                                                                                      tctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGc                              >  2:581159/1‑150 (MQ=255)
                                                                                                                       ctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCa                             >  2:686067/1‑150 (MQ=255)
                                                                                                                         ttCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAg                           <  1:1457926/150‑1 (MQ=255)
                                                                                                                                   ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg                 >  2:1197486/1‑150 (MQ=255)
                                                                                                                                   ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg                 >  1:934575/1‑150 (MQ=255)
                                                                                                                                   ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg                 <  1:174765/150‑1 (MQ=255)
                                                                                                                                   ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcg                 <  2:590973/150‑1 (MQ=255)
                                                                                                                                    catcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgcgc                >  2:712230/1‑150 (MQ=255)
                                                                                                                                       cgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTc             <  1:340328/150‑1 (MQ=255)
                                                                                                                                        gctcggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTca            >  1:129828/7‑150 (MQ=255)
                                                                                                                                        gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTca            >  1:1549893/1‑150 (MQ=255)
                                                                                                                                        gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTca            >  1:1030301/1‑150 (MQ=255)
                                                                                                                                         ctcggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcac           >  1:1408859/6‑150 (MQ=255)
                                                                                                                                         ctcggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcac           >  1:85496/6‑150 (MQ=255)
                                                                                                                                         ctcggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcac           >  2:283731/6‑150 (MQ=255)
                                                                                                                                              tGGCGCGGAGGCGCCACCGAGCGGGAGCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcacgcacg      <  1:1248991/150‑1 (MQ=255)
                                                                                                                                               ggCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGa     <  2:471464/150‑1 (MQ=255)
                                                                                                                                                 cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCAGCGTCTTCTCAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCTGCGCCTCACGCACGAAg   >  2:751480/1‑150 (MQ=255)
                                                                                                                                                 cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAg   >  2:769223/1‑150 (MQ=255)
                                                                                                                                                 cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAg   <  1:624611/150‑1 (MQ=255)
                                                                                                                                                  gcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAGt  >  1:532770/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAGT  >  NC_000913/3882365‑3882660

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: