Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,923,142 G→A intergenic (‑173/+223) pphA ← / ← yebY serine/threonine‑specific protein phosphatase 1/DUF2511 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,923,1420GA97.9% 365.1 / ‑1.6 95intergenic (‑173/+223)pphA/yebYserine/threonine‑specific protein phosphatase 1/DUF2511 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (50/43);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TCACATAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGGAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGTCATCAGC  >  NC_000913/1922995‑1923284
                                                                                                                                                   |                                                                                                                                              
tCACATAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGaaa                                                                                                                                              >  1:689348/1‑150 (MQ=255)
 cacaTAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAAt                                                                                                                                            >  1:671629/1‑151 (MQ=255)
  acaTAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATg                                                                                                                                           >  2:907238/1‑151 (MQ=255)
    aTAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGc                                                                                                                                         >  1:405702/1‑151 (MQ=255)
         cTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTg                                                                                                                                    >  2:154504/1‑151 (MQ=255)
           tAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTggg                                                                                                                                  >  2:691922/1‑151 (MQ=255)
           tAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTgg                                                                                                                                   <  1:1123608/150‑1 (MQ=255)
              cAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGgaga                                                                                                                               <  1:1009687/151‑1 (MQ=255)
               aaGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGgagag                                                                                                                              <  2:1407665/151‑1 (MQ=255)
                aGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGgagaga                                                                                                                             <  1:1628189/151‑1 (MQ=255)
                 gAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGagagag                                                                                                                            <  2:1616671/151‑1 (MQ=255)
                     gTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGt                                                                                                                        >  1:1029554/1‑151 (MQ=255)
                      ttAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGt                                                                                                                        <  2:1019024/150‑1 (MQ=255)
                       tAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTgc                                                                                                                      >  1:1535217/1‑151 (MQ=255)
                        aaaaaaCGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGGAAATGGCTCTTGGGAGAGAGCCGTgcg                                                                                                                     <  2:1714217/151‑1 (MQ=255)
                           aaaCGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCt                                                                                                                   >  1:1391566/1‑150 (MQ=255)
                            aaCGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTa                                                                                                                  >  1:981375/1‑150 (MQ=255)
                             aCGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTaaa                                                                                                                <  2:1196435/151‑1 (MQ=255)
                               gCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTaaa                                                                                                                >  1:1238015/1‑149 (MQ=255)
                               gCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAg                                                                                                              <  2:162345/151‑1 (MQ=255)
                                 tGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGtt                                                                                                            >  1:836791/1‑151 (MQ=255)
                                  ggACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTg                                                                                                           <  2:1103607/151‑1 (MQ=255)
                                    aCTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGc                                                                                                         <  2:1252244/151‑1 (MQ=255)
                                      tCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCAt                                                                                                       >  2:30966/1‑151 (MQ=255)
                                       cAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGATAGAGCCGTGCGCTAAAAGTTGGCa                                                                                                        <  1:1516218/149‑1 (MQ=255)
                                       cAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCa                                                                                                        <  1:303220/149‑1 (MQ=255)
                                        aGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTa                                                                                                     <  1:924015/151‑1 (MQ=255)
                                             gTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCa                                                                                                <  1:980512/151‑1 (MQ=255)
                                                gagTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGc                                                                                             <  1:1100914/151‑1 (MQ=255)
                                                 agTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCt                                                                                            <  2:951604/151‑1 (MQ=255)
                                                  gtgtgtgtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCtt                                                                                           <  2:671687/151‑1 (MQ=255)
                                                   tgtgtgtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCtt                                                                                           >  2:1691861/1‑150 (MQ=255)
                                                   tgtgtgtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCtt                                                                                           >  1:1268264/1‑150 (MQ=255)
                                                    gtgtgtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTa                                                                                          <  1:1379819/150‑1 (MQ=255)
                                                       tgtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttg                                                                                      >  1:875484/1‑151 (MQ=255)
                                                        gtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttgc                                                                                     <  2:1552571/151‑1 (MQ=255)
                                                        gtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttg                                                                                      <  1:483449/150‑1 (MQ=255)
                                                        gtTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttg                                                                                      <  1:1270669/150‑1 (MQ=255)
                                                          tATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttgcc                                                                                    <  2:1349335/150‑1 (MQ=255)
                                                           aTGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttgcctt                                                                                  <  2:1453671/151‑1 (MQ=255)
                                                             ggTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttgccttgc                                                                                >  2:1142772/1‑151 (MQ=255)
                                                             ggTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGttgccttgc                                                                                >  2:944479/1‑151 (MQ=255)
                                                                tGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCt                                                                             >  2:685425/1‑151 (MQ=255)
                                                                   cTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCttt                                                                           >  1:1387665/1‑150 (MQ=255)
                                                                    tataAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCttt                                                                           >  2:534181/1‑149 (MQ=255)
                                                                     ataAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTa                                                                          >  1:1023660/1‑149 (MQ=255)
                                                                      taAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTa                                                                          >  1:86908/1‑148 (MQ=255)
                                                                       aaaGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGaa                                                                      <  1:613507/151‑1 (MQ=255)
                                                                             aGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATg                                                                >  2:471833/1‑151 (MQ=255)
                                                                              gCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATga                                                               <  1:1163581/151‑1 (MQ=255)
                                                                                gAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATgacg                                                             >  2:571553/1‑151 (MQ=255)
                                                                                gAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATgacg                                                             >  1:275705/1‑151 (MQ=255)
                                                                                gAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATgacg                                                             >  1:257335/1‑151 (MQ=255)
                                                                                 aaGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATgacga                                                            <  2:510697/151‑1 (MQ=255)
                                                                                 aaGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATgacg                                                             >  2:1434020/1‑150 (MQ=255)
                                                                                   ggAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGacgacg                                                          >  2:687283/1‑151 (MQ=255)
                                                                                    gAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGc                                                         <  1:1227194/151‑1 (MQ=255)
                                                                                     aaaTGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGcc                                                        <  2:1563439/151‑1 (MQ=255)
                                                                                     aaaTGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGc                                                         >  1:1578568/1‑150 (MQ=255)
                                                                                       aTGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGcc                                                        >  2:1048002/1‑149 (MQ=255)
                                                                                       aTGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGcc                                                        >  2:145295/1‑149 (MQ=255)
                                                                                           ttCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGttt                                                  >  2:580255/1‑151 (MQ=255)
                                                                                           ttCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGttt                                                  >  2:242434/1‑151 (MQ=255)
                                                                                           ttCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGttt                                                  >  1:771541/1‑151 (MQ=255)
                                                                                            tCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGtttt                                                 >  1:995171/1‑151 (MQ=255)
                                                                                            tCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGttt                                                  <  2:1189741/150‑1 (MQ=255)
                                                                                            tCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGttt                                                  <  2:486252/150‑1 (MQ=255)
                                                                                               ggCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTcc                                               <  2:1698392/150‑1 (MQ=255)
                                                                                                 cTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGt                                            >  2:498153/1‑151 (MQ=255)
                                                                                                    ttAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGttt                                          <  2:1482278/150‑1 (MQ=255)
                                                                                                       aCAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGt                                      <  1:1218257/151‑1 (MQ=255)
                                                                                                        cAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTg                                     <  1:1048038/151‑1 (MQ=255)
                                                                                                         aGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTg                                     >  1:978599/1‑150 (MQ=255)
                                                                                                         aGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGc                                    >  1:1072858/1‑151 (MQ=255)
                                                                                                         aGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGc                                    >  1:1372588/1‑151 (MQ=255)
                                                                                                             aaaaaGAGACCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCaaa                                 >  2:648147/1‑150 (MQ=255)
                                                                                                                  gagaCCGAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTc                           <  1:312976/151‑1 (MQ=255)
                                                                                                                   agaCCGAACACGTTTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTc                           >  1:1044678/1‑150 (MQ=255)
                                                                                                                        gAACACGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATaa                      >  1:71995/1‑150 (MQ=255)
                                                                                                                           cacGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATaaaa                    >  2:485148/1‑149 (MQ=255)
                                                                                                                           cacGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAg                  <  2:689407/151‑1 (MQ=255)
                                                                                                                            acGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGc                 <  1:1549771/151‑1 (MQ=255)
                                                                                                                             cGATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCg                >  2:412922/1‑151 (MQ=255)
                                                                                                                              gATTCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCgt               >  1:1320787/1‑151 (MQ=255)
                                                                                                                                ttCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCgtgg             >  1:953639/1‑151 (MQ=255)
                                                                                                                                 tCCTGTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCgtggt            >  2:806085/1‑151 (MQ=255)
                                                                                                                                     gTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAACTGGTCAATAAAAAGCGTGGTGGtc        >  2:1375634/1‑151 (MQ=255)
                                                                                                                                     gTATTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGtc        <  1:1141887/151‑1 (MQ=255)
                                                                                                                                       aTTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGtcat      <  2:1535315/151‑1 (MQ=255)
                                                                                                                                       aTTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGtcat      <  2:316774/151‑1 (MQ=255)
                                                                                                                                       aTTCGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGtcat      >  1:1150998/1‑151 (MQ=255)
                                                                                                                                          cGGTCCAGGGAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGTCATCAg   <  1:343111/151‑1 (MQ=255)
                                                                                                                                          cGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGtcatca    <  2:613639/150‑1 (MQ=255)
                                                                                                                                          cGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGTCATCAg   <  1:359311/151‑1 (MQ=255)
                                                                                                                                          cGGTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGACAATAAAAAGCGTGGTGGTCATCAg   <  2:980487/151‑1 (MQ=255)
                                                                                                                                           ggTCCAGGAAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGTCATCAGc  >  1:1532132/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCACATAATCTTAACAAGAATGTTAAAAAACGCTGGACTCAGACAGTAGAGTGTGTGTTATGGTTGACTATAAAGTCAGCGAAGGAAATGCTTCTGGCTTTTAACAGATAAAAAGAGACCGAACACGATTCCTGTATTCGGTCCAGGGAAATGGCTCTTGGGAGAGAGCCGTGCGCTAAAAGTTGGCATTAATGCAGGCTTAGTTGCCTTGCCCTTTAAGAATAGATGACGACGCCAGGTTTTCCAGTTTGCGTGCAAAATGGTCAATAAAAAGCGTGGTGGTCATCAGC  >  NC_000913/1922995‑1923284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: