Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,233 A→T 5.8% intergenic (+36/+111) yjdP → / ← phnP putative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,314,2330AT 5.8% 410.5 / 9.8 119intergenic (+36/+111)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (55/57);  new base T (4/3);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 121 total aligned reads displayed.

CCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACC  >  NC_000913/4314086‑4314380
                                                                                                                                                   |                                                                                                                                                   
ccAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCaa                                                                                                                                                     <  2:159811/148‑1 (MQ=255)
 cAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATg                                                                                                                                                   <  4:20068/149‑1 (MQ=255)
   gTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGTCATCGCCCCTCACCCCAATggg                                                                                                                                                 <  4:534283/149‑1 (MQ=255)
        aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGGGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagag                                                                                                                                            <  3:607919/149‑1 (MQ=255)
        aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagag                                                                                                                                            <  3:901953/149‑1 (MQ=255)
        aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggaga                                                                                                                                             <  1:694213/148‑1 (MQ=255)
          gATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgaggg                                                                                                                                          <  3:897438/149‑1 (MQ=255)
          gATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagagg                                                                                                                                           <  2:639180/148‑1 (MQ=255)
           aTCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagaggga                                                                                                                                         >  4:673476/1‑149 (MQ=255)
            tCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagaggga                                                                                                                                         >  1:141223/1‑148 (MQ=255)
            tCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCATCACCCCAATGgggagag                                                                                                                                            <  2:568420/145‑1 (MQ=255)
            tCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCATCACCCCAATGgggagag                                                                                                                                            >  1:568420/1‑145 (MQ=255)
                cgcAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGaaa                                                                                                                                     >  2:953406/1‑148 (MQ=255)
                   aggcagCTGGACGATCGCGAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                 <  4:87301/149‑1 (MQ=255)
                   aggcagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                  <  2:619338/148‑1 (MQ=255)
                   aggcagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                 <  4:290301/149‑1 (MQ=255)
                     gcagCTGGACGATCGCCAGCGCCAGCTTGATCGGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                               <  4:174391/149‑1 (MQ=255)
                            gACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAata                                                                                                                         >  2:167982/1‑148 (MQ=255)
                              cGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATAtt                                                                                                                       <  1:335975/148‑1 (MQ=255)
                              cGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTc                                                                                                                      <  3:170018/149‑1 (MQ=255)
                              cGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTc                                                                                                                      <  4:222288/149‑1 (MQ=255)
                               gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAAc                                                                                                                                   <  4:502136/135‑1 (MQ=255)
                               gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAAc                                                                                                                                   >  3:502136/1‑135 (MQ=255)
                               gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                     >  4:337022/1‑149 (MQ=255)
                               gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                     <  3:673476/149‑1 (MQ=255)
                               gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                     >  4:610900/1‑149 (MQ=255)
                                aTCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCaa                                                                                                                    <  3:378864/149‑1 (MQ=255)
                                 tCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAat                                                                                                                   >  3:826548/1‑149 (MQ=255)
                                  cgccagcgccagCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagggg                                                                                                                                         >  2:416809/1‑125 (MQ=255)
                                  cgccagcgccagCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagggg                                                                                                                                         <  1:416809/126‑2 (MQ=255)
                                      agcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                               <  2:31157/148‑1 (MQ=255)
                                      agcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGa                                                                                                              >  3:761298/1‑149 (MQ=255)
                                       gcgccagcTTAATCAGGAGCGTCGGAAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGTTGAGGTTGGGTATGGAGAAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGa                                                                                                              <  1:944609/148‑1 (MQ=255)
                                           cagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGcac                                                                                                          >  2:875799/1‑148 (MQ=255)
                                            agcTTGGTCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCAAGcaca                                                                                                         >  2:70027/1‑148 (MQ=255)
                                             gcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACaa                                                                                                        >  1:448745/1‑148 (MQ=255)
                                                tGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCg                                                                                                     >  2:446900/1‑148 (MQ=255)
                                                  aTCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTc                                                                                                  <  4:530123/149‑1 (MQ=255)
                                                   tCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGTACAATCGGTc                                                                                                  >  2:826780/1‑148 (MQ=255)
                                                     aGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTccc                                                                                                <  2:141223/148‑1 (MQ=255)
                                                       gAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAGTCGGTCCCCt                                                                                              >  1:134273/1‑148 (MQ=255)
                                                          cGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGc                                                                                           >  2:679460/1‑148 (MQ=255)
                                                          cGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGc                                                                                           >  2:442632/1‑148 (MQ=255)
                                                            tCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGCGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGcccc                                                                                        <  3:623857/149‑1 (MQ=255)
                                                               gCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCctct                                                                                     >  3:202591/1‑149 (MQ=255)
                                                                 aaCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                   <  4:431519/149‑1 (MQ=255)
                                                                  aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTggg                                                                                  <  4:922337/149‑1 (MQ=255)
                                                                  aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                   <  1:448747/148‑1 (MQ=255)
                                                                    tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgggg                                                                                 <  2:448745/148‑1 (MQ=255)
                                                                     ggaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                >  2:38832/1‑148 (MQ=255)
                                                                      gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                              <  3:950832/149‑1 (MQ=255)
                                                                       aggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                              >  1:733903/1‑148 (MQ=255)
                                                                       aggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                              <  1:268740/148‑1 (MQ=255)
                                                                          aTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                              <  2:193759/145‑1 (MQ=255)
                                                                          aTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                              >  1:193759/1‑145 (MQ=255)
                                                                            gAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                        <  4:676272/149‑1 (MQ=255)
                                                                             aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                 >  3:135973/1‑75 (MQ=255)
                                                                             aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                 <  4:135973/75‑1 (MQ=255)
                                                                             aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCctc                                                                                      >  1:472884/1‑134 (MQ=255)
                                                                             aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCctc                                                                                      <  2:472884/134‑1 (MQ=255)
                                                                             aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                        <  2:9289/148‑1 (MQ=255)
                                                                                  cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                 <  4:918590/70‑1 (MQ=255)
                                                                                  cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                 >  3:918590/1‑70 (MQ=255)
                                                                                      gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                              <  4:35177/149‑1 (MQ=255)
                                                                                         tGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGcc                                                           >  4:244416/1‑149 (MQ=255)
                                                                                         tGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGc                                                            >  1:930232/1‑148 (MQ=255)
                                                                                         tGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGc                                                            >  2:312064/1‑148 (MQ=255)
                                                                                          ggAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGcc                                                           >  2:634573/1‑148 (MQ=255)
                                                                                          ggAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCa                                                          >  4:469788/1‑149 (MQ=255)
                                                                                             aGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAgcg                                                       >  3:589414/1‑149 (MQ=255)
                                                                                                    tATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                <  3:866674/133‑1 (MQ=255)
                                                                                                    tATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                >  4:866674/1‑133 (MQ=255)
                                                                                                          cGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                          <  3:949283/149‑1 (MQ=255)
                                                                                                           gagcaGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACTATCGCAATAGTCAATATCGAGCACAATCGGACGCCTCGCCCCGCTTGGGAGAGGGTTAGGGTGAGGGTGCCAGCGCGCACATTCCctct                                         >  4:500114/5‑149 (MQ=255)
                                                                                                           gATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                         >  3:394158/1‑149 (MQ=255)
                                                                                                           gATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                         >  3:78947/1‑149 (MQ=255)
                                                                                                            aTGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                         <  2:848438/148‑1 (MQ=255)
                                                                                                             tGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCtctc                                        <  1:211073/148‑1 (MQ=255)
                                                                                                              gAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcat                                      <  3:500114/149‑1 (MQ=255)
                                                                                                               aGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcat                                      <  1:591417/148‑1 (MQ=255)
                                                                                                                  ttGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCAc                                   >  1:934359/1‑148 (MQ=255)
                                                                                                                           gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCa                          >  1:269246/1‑148 (MQ=255)
                                                                                                                              cAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAtct                      <  3:738993/149‑1 (MQ=255)
                                                                                                                                ggCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcc                    <  3:469788/149‑1 (MQ=255)
                                                                                                                                  cATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                  <  4:31827/149‑1 (MQ=255)
                                                                                                                                  cATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcca                   <  1:122837/148‑1 (MQ=255)
                                                                                                                                  cATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcca                   >  2:815318/1‑148 (MQ=255)
                                                                                                                                   aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGc                                 >  3:316929/1‑133 (MQ=255)
                                                                                                                                      gCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca              <  3:610900/149‑1 (MQ=255)
                                                                                                                                       ccccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca              >  1:248286/1‑148 (MQ=255)
                                                                                                                                          ctctCCCCATTGGGGTGAGGGGGAAAACGAGCGCAATATTCAATATc                                                                                                                <  4:346320/47‑1 (MQ=21)
                                                                                                                                          ctctCCCCATTGGGGTGAGGGGGAAAACGAGCGCAATATTCAATATc                                                                                                                >  3:346320/1‑47 (MQ=21)
                                                                                                                                              ccccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcc       >  2:457400/1‑148 (MQ=255)
                                                                                                                                               cccAATGGGGAGAGGGAGAAAACGAGTGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGCGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcccc     >  3:942372/1‑149 (MQ=255)
                                                                                                                                               cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                   <  3:766440/71‑1 (MQ=255)
                                                                                                                                               cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcccc     >  4:475552/1‑149 (MQ=255)
                                                                                                                                               cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcccc     >  3:426888/1‑149 (MQ=255)
                                                                                                                                                ccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcccc     >  2:942939/1‑148 (MQ=255)
                                                                                                                                                 cAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAc   >  3:810324/1‑149 (MQ=255)
                                                                                                                                                  aaTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAcc  >  4:374283/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACC  >  NC_000913/4314086‑4314380

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