Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,491,905 A→T 8.5% intergenic (+55/+17) 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,9050AT 8.5% 227.2 / 7.0 77intergenic (+55/+17)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 = 1.00e+00

TTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACC  >  NC_000913/1491757‑1492051
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                               >  1:1380197/1‑150 (MQ=255)
                                    cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                               <  2:1382299/150‑1 (MQ=255)
                                     aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                              <  2:1035430/150‑1 (MQ=255)
                                     aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                              <  2:455872/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                                                                   <  2:1102328/131‑1 (MQ=255)
                                                                                                    aaaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGt                                                                   >  1:1102328/1‑131 (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:905857/1‑150 (MQ=255)
                                                                                                                   agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAg                                 <  1:1971782/150‑1 (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:1568707/4‑150 (MQ=255)
                                                                                                                                        caagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAg            >  1:1201927/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)
                                                                                                                                                     |                                                                                                                                                  
TTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCG‑AAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACC  >  NC_000913/1491757‑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: ‑