Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,491,904 A→T 9.0% intergenic (+54/+18) ydcA → / ← hokB protein YdcA/toxin HokB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,491,9040AT 9.0% 230.3 / 9.4 77intergenic (+54/+18)ydcA/hokBprotein YdcA/toxin HokB
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base T (3/4);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACC  >  NC_000913/1491755‑1492051
                                                                                                                                                      |                                                                                                                                                   
ggTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaa                                                                                                                                                     >  2:375485/1‑150 (MQ=255)
  ttCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATTCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                   <  1:1477904/150‑1 (MQ=255)
  ttCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                   <  2:748572/150‑1 (MQ=255)
   tCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAAt                                                                                                                                                  >  1:442075/1‑150 (MQ=255)
     gAAAGGGGGGATCTCGCACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcc                                                                                                                                                >  1:283283/1‑150 (MQ=255)
     gAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcc                                                                                                                                                >  1:296456/1‑150 (MQ=255)
            gggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgca                                                                                                                                         <  2:1363745/150‑2 (MQ=255)
              gATCTCACACTGTACGGCCGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAACACATCCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgcatt                                                                                                                                       <  2:558250/150‑4 (MQ=255)
              gATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgcatt                                                                                                                                       <  2:761316/150‑4 (MQ=255)
               aTCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgcattt                                                                                                                                      <  1:1436876/150‑5 (MQ=255)
                 ctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                    <  2:340782/150‑1 (MQ=255)
                   cacaCTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCtt                                                                                                                                  >  2:1548139/1‑150 (MQ=255)
                       cTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTAcc                                                                                                                               <  1:1177885/150‑1 (MQ=255)
                       cTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCt                                                                                                                              >  1:1577127/1‑150 (MQ=255)
                         gTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGGAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                                                                            >  2:1335500/1‑150 (MQ=255)
                         gTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                                                                            >  2:1010240/1‑150 (MQ=255)
                          tACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGCGGAAAATCCCCGCGCTTGCCCTTACCTGGa                                                                                                                           <  1:217546/150‑1 (MQ=255)
                           aCGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGAc                                                                                                                          <  2:718779/150‑1 (MQ=255)
                            cGGCAGGCTGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACg                                                                                                                         <  1:1924926/150‑1 (MQ=255)
                                  gcggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAgg                                                                                                                   >  2:1825639/1‑150 (MQ=255)
                                   cggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                  >  1:1382874/1‑150 (MQ=255)
                                   cggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                  <  1:195061/150‑1 (MQ=255)
                                      cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                               <  2:1382299/150‑1 (MQ=255)
                                      cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                               >  1:1380197/1‑150 (MQ=255)
                                       aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                              <  2:455872/150‑1 (MQ=255)
                                       aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                              <  2:1035430/150‑1 (MQ=255)
                                       aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATCCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                              <  2:855640/150‑1 (MQ=255)
                                                tGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGGAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                                                                     <  1:1335500/150‑1 (MQ=255)
                                                tGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                  >  1:1168485/1‑137 (MQ=255)
                                                tGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGc                                                                                                                  <  2:1168485/137‑1 (MQ=255)
                                                tGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATCAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCCCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                                                                     >  2:420661/1‑150 (MQ=255)
                                                       atgGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctcct                                                                                              <  1:380117/150‑1 (MQ=255)
                                                         ggTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTa                                                                                            <  2:1436580/150‑1 (MQ=255)
                                                             cTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCAc                                                                                        >  2:195761/1‑150 (MQ=255)
                                                               aTTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCg                                                                                      <  1:1548139/150‑1 (MQ=255)
                                                                    tGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCg                                                                                 <  1:195761/150‑1 (MQ=255)
                                                                    tGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCg                                                                                 <  1:554158/150‑1 (MQ=255)
                                                                     gCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCgg                                                                                >  2:1097422/1‑150 (MQ=255)
                                                                      cATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGa                                                                               >  1:1706120/1‑150 (MQ=255)
                                                                       aTCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCAAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGaa                                                                              >  1:843041/1‑150 (MQ=255)
                                                                            aaaaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGc                                                                          <  2:442075/150‑1 (MQ=255)
                                                                              aaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAg                                                                        >  1:1778817/1‑150 (MQ=255)
                                                                                  aCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCg                                                                    <  2:405079/150‑1 (MQ=255)
                                                                                  aCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCg                                                                    <  1:1097422/150‑1 (MQ=255)
                                                                                       cACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTagag                                                               <  2:1380197/150‑1 (MQ=255)
                                                                                       cACTAACTGAAGTGTAAAAGGGGAGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGCGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTACCGGAACCGCAGTTCCAagag                                                               >  1:1798153/1‑150 (MQ=255)
                                                                                         cTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCg                                                             <  2:1357356/150‑1 (MQ=255)
                                                                                           aaCTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGtt                                                           >  1:475771/1‑150 (MQ=255)
                                                                                            aCTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGttt                                                          >  1:277139/1‑150 (MQ=255)
                                                                                              tGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt                                                        >  1:108729/1‑150 (MQ=255)
                                                                                               gAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTc                                                       >  1:1359057/1‑150 (MQ=255)
                                                                                                 aGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCgg                                                     >  2:1926819/1‑150 (MQ=255)
                                                                                                 aGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCgg                                                     >  1:507030/1‑150 (MQ=255)
                                                                                                    gtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGta                                                  >  2:846622/1‑149 (MQ=255)
                                                                                                    gtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGGTCGTAGAGCGTTTGTCGGGTc                                                  >  2:1304693/1‑150 (MQ=255)
                                                                                                      aaaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGt                                                                   >  1:1102328/1‑131 (MQ=255)
                                                                                                      aaaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGt                                                                   <  2:1102328/131‑1 (MQ=255)
                                                                                                       aaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGCCGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAgg                                               >  2:1948876/1‑150 (MQ=255)
                                                                                                         aGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAg                                             <  1:128661/150‑1 (MQ=255)
                                                                                                            ggTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAgtgt                                          <  1:375485/150‑1 (MQ=255)
                                                                                                                 cATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATg                                     >  2:370660/1‑150 (MQ=255)
                                                                                                                   tGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTc                                   <  2:1577127/150‑1 (MQ=255)
                                                                                                                    gagaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCa                                  >  1:1310495/1‑150 (MQ=255)
                                                                                                                    gagaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCa                                  >  2:1377518/1‑150 (MQ=255)
                                                                                                                     agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAg                                 <  1:1971782/150‑1 (MQ=255)
                                                                                                                     agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAg                                 >  1:905857/1‑150 (MQ=255)
                                                                                                                        aaaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAAt                              >  1:1676381/1‑150 (MQ=255)
                                                                                                                               tCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATg                       <  2:108729/150‑1 (MQ=255)
                                                                                                                                   cTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGCGAATGTCAGAATCGTAATGCAGa                   >  2:1444334/1‑150 (MQ=255)
                                                                                                                                       tCAAATGCGCGCGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAAt               <  1:846622/150‑1 (MQ=255)
                                                                                                                                        cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg              >  1:1215693/1‑150 (MQ=255)
                                                                                                                                         aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa             >  2:185867/1‑150 (MQ=255)
                                                                                                                                          caagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAg            >  1:1201927/4‑150 (MQ=255)
                                                                                                                                          caagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAg            >  1:1568707/4‑150 (MQ=255)
                                                                                                                                            tGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCa          >  1:1876914/1‑150 (MQ=255)
                                                                                                                                             gcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg         >  1:891568/1‑150 (MQ=255)
                                                                                                                                              cgcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa        >  2:494780/1‑150 (MQ=255)
                                                                                                                                                  gggAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGTGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacaca    <  2:283283/150‑1 (MQ=255)
                                                                                                                                                    gAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACAcacc  >  2:135427/1‑150 (MQ=255)
                                                                                                                                                    gAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACAcacc  >  2:1048567/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACC  >  NC_000913/1491755‑1492051

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