Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,093,922 A→T 5.8% L82Q (CTG→CAG)  ptsN ← phosphotransferase system subunit IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,093,9220AT 5.8% 196.1 / 2.7 69L82Q (CTG→CAG) ptsNphosphotransferase system subunit IIA
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base T (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GGACCTCCTTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTGGTCGGCGA  >  NC_002947/1093776‑1094069
                                                                                                                                                  |                                                                                                                                                    
ggACCTCCTTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTgcagc                                                                                                                                                    <  7:253699/149‑1 (MQ=255)
 gACCTCCTTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTgcagca                                                                                                                                                   <  2:33840/149‑1 (MQ=255)
 gACCTCCTTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTgcagca                                                                                                                                                   <  2:27700/149‑1 (MQ=255)
       ctTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGgcgc                                                                                                                                             >  2:327217/1‑149 (MQ=255)
          gCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGa                                                                                                                                          >  1:94414/1‑149 (MQ=255)
                    aTGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGcgcg                                                                                                                                 >  5:32637/1‑147 (MQ=255)
                    aTGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGa                                                                                                                               <  2:94414/149‑1 (MQ=255)
                    aTGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGa                                                                                                                               <  4:258083/149‑1 (MQ=255)
                       cgtgCGATCTGCCGCAGAAGTTCCCGGTGGGCATCGGTGGCCGCTTCGCGGACCAGCAGCACAAACACCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAAGCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGc                                                                                                                             <  8:189967/145‑1 (MQ=255)
                        tGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCa                                                                                                                            >  4:175368/1‑148 (MQ=255)
                          gCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATcc                                                                                                                         >  8:257832/1‑149 (MQ=255)
                            gATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCt                                                                                                                        >  3:338835/1‑148 (MQ=255)
                                tGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAg                                                                                                                   >  6:308536/1‑149 (MQ=255)
                                tGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAg                                                                                                                   >  8:217459/1‑149 (MQ=255)
                                 gACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGc                                                                                                                  >  7:332507/1‑149 (MQ=255)
                                  aCGAAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCAACCGGCGCGCCATCCAGGGCGTCGTAATCGATAGGCGCTTCCACGTGCAGCAGGGCGCTGACGGGCGCGGAGCATCAGTCAAgac                                                                                                                 >  2:94823/1‑147 (MQ=255)
                                      aGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCa                                                                                                             <  5:312019/149‑1 (MQ=255)
                                            tCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTgg                                                                                                       >  7:136359/1‑149 (MQ=255)
                                              aaGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTggg                                                                                                      <  6:32637/148‑1 (MQ=255)
                                               aGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGAt                                                                                                    >  2:406347/1‑149 (MQ=255)
                                               aGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGAt                                                                                                    <  1:327217/149‑1 (MQ=255)
                                                      taTCGGTGTCCGCTTCGGGTAGCAGCAGGACAAACAGCAGGTCCACCGTCGCGCCCTCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCAt                                                                                                                           <  1:326124/118‑1 (MQ=255)
                                                      taTCGCTGGCCGATTCGGGTACCAGCACCACAAAAACCAGGTCCACCGGCGCGCCATCCATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGAAGGTCGCTGAATGGCGCGGAGCAt                                                                                                                           >  2:326124/2‑119 (MQ=255)
                                                      cATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATg                                                                                             <  3:16453/149‑1 (MQ=255)
                                                       aTCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCTGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGc                                                                                            >  3:307536/1‑149 (MQ=255)
                                                       aTCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCTGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGc                                                                                            >  3:312997/1‑149 (MQ=255)
                                                          ggTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCAt                                                                                         >  5:415096/1‑149 (MQ=255)
                                                          ggTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCAt                                                                                         >  6:296921/1‑149 (MQ=255)
                                                           gTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCATGCCGGCAGTGTGGGATGGCGATGCCAtt                                                                                        >  7:186906/1‑149 (MQ=255)
                                                           gTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCAtt                                                                                        >  3:45351/1‑149 (MQ=255)
                                                             ggCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGTGATGCCATTGc                                                                                      <  1:406347/149‑1 (MQ=255)
                                                              gCCGGTTGGGGTACCAGCAGCAAAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGGAGGGCGCTGACTGCGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGc                                                                                      <  1:94823/149‑1 (MQ=255)
                                                              gCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGcc                                                                                     >  7:312414/1‑149 (MQ=255)
                                                              gCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGcc                                                                                     <  5:296921/149‑1 (MQ=255)
                                                                cGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCgg                                                                                                                                <  6:323794/104‑1 (MQ=255)
                                                                cGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCgg                                                                                                                                >  5:323794/1‑104 (MQ=255)
                                                                cGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGc                                                                                      >  3:339202/1‑146 (MQ=255)
                                                                cGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGc                                                                                      <  4:339202/146‑1 (MQ=255)
                                                                 gCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGa                                                                                   >  3:322076/1‑148 (MQ=255)
                                                                      gggTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGCTAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCg                                                                              <  6:415096/148‑1 (MQ=255)
                                                                       ggTACCAGCAGCACAGACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCgg                                                                             >  4:174839/1‑148 (MQ=255)
                                                                       ggTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGt                                                                            <  8:312414/149‑1 (MQ=255)
                                                                         tACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAg                                                                          <  8:332507/149‑1 (MQ=255)
                                                                           ccAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGaa                                                                        >  8:48147/1‑149 (MQ=255)
                                                                            cagcagCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAAc                                                                       >  3:135476/1‑149 (MQ=255)
                                                                              gcagcaCAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAAcc                                                                      >  4:79774/1‑148 (MQ=255)
                                                                                       aCAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTcc                                                            <  5:308536/149‑1 (MQ=255)
                                                                                        cagcagGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTccc                                                           <  3:175368/149‑1 (MQ=255)
                                                                                           cagGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCggg                                                        <  2:194351/149‑1 (MQ=255)
                                                                                            agGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGc                                                       <  8:136359/149‑1 (MQ=255)
                                                                                                 cACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCAGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTAc                                                    >  2:26540/1‑147 (MQ=255)
                                                                                                   ccGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGt                                                >  2:52467/1‑149 (MQ=255)
                                                                                                    cGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGtt                                               <  1:385885/149‑1 (MQ=255)
                                                                                                        gcgcCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTc                                           <  4:135476/149‑1 (MQ=255)
                                                                                                         cgcCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCGGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCg                                          <  1:26540/149‑1 (MQ=255)
                                                                                                           ccATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGTGCTACCAGTTTTTCGaa                                        <  6:344931/149‑1 (MQ=255)
                                                                                                           ccATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGaa                                        <  1:52467/149‑1 (MQ=255)
                                                                                                            cATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAg                                       <  3:79774/149‑1 (MQ=255)
                                                                                                                 aTGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTc                                  <  6:144193/149‑1 (MQ=255)
                                                                                                                     cgtcgtAATCGATAGGCGCTTCCAGGTGCTGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGc                              <  4:312997/149‑1 (MQ=255)
                                                                                                                     cgtcgtAATCGATAGGCGCTTCCAGGTGCTGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGc                              <  4:307536/149‑1 (MQ=255)
                                                                                                                     cgtcgtAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGc                              <  4:45351/149‑1 (MQ=255)
                                                                                                                      gtcgtAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCt                                <  5:273907/146‑1 (MQ=255)
                                                                                                                      gtcgtAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCt                                >  6:273907/1‑146 (MQ=255)
                                                                                                                       tcgtAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCAt                            >  6:85333/1‑149 (MQ=255)
                                                                                                                            aTCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCCTCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCAtctc                         >  3:210505/1‑147 (MQ=255)
                                                                                                                                 tAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCt                  <  4:338835/149‑1 (MQ=255)
                                                                                                                                    gcgcTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTgg                >  3:45838/1‑148 (MQ=255)
                                                                                                                                     cgcTTCCAGGTGAAGCAGGGCGCTGACTGGCGCGGATCATCCTTCAAGCCGGCAGTGTGGTATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACTTCCTGCATCTCCAGTTCTCGCa              >  1:69320/1‑149 (MQ=255)
                                                                                                                                        ttCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACtt           >  4:275227/1‑149 (MQ=255)
                                                                                                                                          ccAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTgg         >  2:98052/1‑149 (MQ=255)
                                                                                                                                           cAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTGGt        >  2:210876/1‑149 (MQ=255)
                                                                                                                                                 cagcagGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTGGTCGGCGa  <  5:85333/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGACCTCCTTGCGATCGAGCATGCTGGCGATCTGACGCAGAAGTTCAAGGTGGGCATCGGTGGCCGCTTCGGGTACCAGCAGCACAAACAGCAGGTCCACCGGCGCGCCATCGATGGCGTCGTAATCGATAGGCGCTTCCAGGTGCAGCAGGGCGCTGACTGGCGCGGAGCATCCTTCAAGCCGGCAGTGTGGGATGGCGATGCCATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTGGTCGGCGA  >  NC_002947/1093776‑1094069

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: