Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,419,395 C→G 7.6% V132L (GTG→CTG)  yfcE ← phosphodiesterase activity on bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,419,3950CG 7.6% 318.8 / 14.3 92V132L (GTG→CTG) yfcEphosphodiesterase activity on bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base G (5/2);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419247‑2419542
                                                                                                                                                    |                                                                                                                                                   
gcgACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAc                                                                                                                                                     <  4:722114/149‑1 (MQ=255)
 cgACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGCATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCAc                                                                                                                                                     <  1:588140/148‑1 (MQ=255)
  gACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACt                                                                                                                                                    <  2:10086/148‑1 (MQ=255)
  gACCTGATCAATGCTACTATCGTCATTGAGTGCGATTACACTGGGTACGTCATTATTCAGCATGCCATACCTCGCCTGACTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTTGAAAAGATCCCCCCGCTGTTCTGCCCCt                                                                                                                                                    <  2:743937/148‑1 (MQ=255)
      tGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCa                                                                                                                                                >  1:246079/1‑148 (MQ=255)
        cgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGAt                                                                                                                                             <  4:580401/149‑1 (MQ=255)
         gcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGAt                                                                                                                                             <  2:355516/148‑1 (MQ=255)
          cAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  3:416556/149‑1 (MQ=255)
          cAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATg                                                                                                                                            <  2:631405/148‑1 (MQ=255)
            aTGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGt                                                                                                                                         <  4:798300/149‑1 (MQ=255)
             tGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGt                                                                                                                                         <  1:40272/148‑1 (MQ=255)
               aTACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTAt                                                                                                                                       <  1:572805/148‑1 (MQ=255)
                   tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCg                                                                                                                                  >  3:78540/1‑149 (MQ=255)
                    ttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCg                                                                                                                                  >  1:406409/1‑148 (MQ=255)
                       aTCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTac                                                                                                                               >  2:25123/1‑148 (MQ=255)
                        tCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTacac                                                                                                                             >  4:745628/1‑149 (MQ=255)
                         cATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTg                                                                                                                                                   >  1:104677/1‑126 (MQ=255)
                         cATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTg                                                                                                                                                   <  2:104677/126‑1 (MQ=255)
                            tGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAg                                                                                                                          >  2:198153/1‑148 (MQ=255)
                             gAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGc                                                                                                                         <  2:261608/148‑1 (MQ=255)
                                tGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACAcc                                                                                                                            <  4:503586/142‑1 (MQ=255)
                                tGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACAcc                                                                                                                            >  3:503586/1‑142 (MQ=255)
                                tGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACAt                                                                                                                     >  4:426324/1‑149 (MQ=255)
                                tGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACAt                                                                                                                     >  3:57926/1‑149 (MQ=255)
                                     ttACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGtt                                                                                                                 >  1:244171/1‑148 (MQ=255)
                                      tACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTc                                                                                                                <  2:406409/148‑1 (MQ=255)
                                       acacTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                                                               <  1:275511/148‑1 (MQ=255)
                                            aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGttt                                                                                                          >  2:727433/1‑148 (MQ=255)
                                            aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGt                                                                                                            >  1:98309/1‑146 (MQ=255)
                                            aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGt                                                                                                            <  2:98309/146‑1 (MQ=255)
                                            aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTa                                                                                                         <  3:745628/149‑1 (MQ=255)
                                              gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTaca                                                                                                                              >  3:450736/1‑126 (MQ=255)
                                              gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTaca                                                                                                                              <  4:450736/126‑1 (MQ=255)
                                               tACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAg                                                                                                      >  3:402852/1‑149 (MQ=255)
                                                 cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgca                                                                                                    >  3:878672/1‑149 (MQ=255)
                                                 cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgc                                                                                                     <  2:648589/148‑1 (MQ=255)
                                                   tCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  2:577315/1‑148 (MQ=255)
                                                   tCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  1:348517/1‑148 (MQ=255)
                                                        aTCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAt                                                                                             >  3:789557/1‑149 (MQ=255)
                                                         tCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAt                                                                                             <  1:697738/148‑1 (MQ=255)
                                                          ccAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAttt                                                                                           <  3:100604/149‑1 (MQ=255)
                                                            aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAtttt                                                                                          <  1:198153/148‑1 (MQ=255)
                                                                   cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaa                                                                                                        <  2:220411/127‑1 (MQ=255)
                                                                   cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaa                                                                                                        >  1:220411/1‑127 (MQ=255)
                                                                   cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCa                                                                                  >  4:156908/1‑149 (MQ=255)
                                                                    aTAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCa                                                                                  >  2:803057/1‑148 (MQ=255)
                                                                    aTAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCa                                                                                  >  2:559368/1‑148 (MQ=255)
                                                                       aCTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaa                                                                               <  1:11390/148‑1 (MQ=255)
                                                                       aCTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaa                                                                               <  2:379540/148‑1 (MQ=255)
                                                                            ccGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATg                                                                          >  1:800157/1‑148 (MQ=255)
                                                                            ccGGATTAACGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACAACAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGa                                                                         <  4:422640/149‑1 (MQ=255)
                                                                                      gCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCa                                                                      >  1:229486/1‑142 (MQ=255)
                                                                                      gCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCa                                                                      <  2:229486/142‑1 (MQ=255)
                                                                                          ttCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCaaaa                                                            >  2:429209/1‑148 (MQ=255)
                                                                                            cGGAATACTCACCGAGCCGGGGTTGAAGTGGAAATTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACAGCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAAc                                                          <  1:559368/148‑1 (MQ=255)
                                                                                            cGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACa                                                         <  4:900876/149‑1 (MQ=255)
                                                                                            cGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACa                                                         <  3:426324/149‑1 (MQ=255)
                                                                                             ggAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACa                                                         >  1:37236/1‑148 (MQ=255)
                                                                                                  aCTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgt                                                    <  2:246079/148‑1 (MQ=255)
                                                                                                    tCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgt                                                 <  3:156908/149‑1 (MQ=255)
                                                                                                       ccGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGtttt                                              >  3:883940/1‑149 (MQ=255)
                                                                                                         gAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttttt                                             >  1:319364/1‑148 (MQ=255)
                                                                                                           gCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTcc                                           >  1:39596/1‑148 (MQ=255)
                                                                                                            ccGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa                                          <  1:577315/148‑1 (MQ=255)
                                                                                                                 gTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAAt                                     <  1:803057/148‑1 (MQ=255)
                                                                                                                   tGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAcc                                  <  4:883940/149‑1 (MQ=255)
                                                                                                                      aGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTg                                <  1:518628/148‑1 (MQ=255)
                                                                                                                      aGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTATTCCAGTAATACCTg                                >  2:264299/1‑148 (MQ=255)
                                                                                                                         ggAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTg                             >  1:383223/1‑148 (MQ=255)
                                                                                                                         ggAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTg                             >  2:46199/1‑148 (MQ=255)
                                                                                                                          gAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGcc                           >  3:566659/1‑149 (MQ=255)
                                                                                                                           aaaaTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            >  1:798617/1‑131 (MQ=255)
                                                                                                                           aaaaTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTc                                            <  2:798617/131‑1 (MQ=255)
                                                                                                                               tctcCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACgg                       <  2:37236/148‑1 (MQ=255)
                                                                                                                                  ccccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATa                                    >  1:746780/1‑132 (MQ=255)
                                                                                                                                  ccccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATa                                    <  2:746780/132‑1 (MQ=255)
                                                                                                                                  ccccccGCTGCTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGgcgc                    >  1:288014/1‑148 (MQ=255)
                                                                                                                                    ccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAg                                                        >  1:701495/1‑110 (MQ=255)
                                                                                                                                    ccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAg                                                        <  2:701495/110‑1 (MQ=255)
                                                                                                                                     cccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGt                               <  1:155387/134‑1 (MQ=255)
                                                                                                                                     cccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGt                               >  2:155387/1‑134 (MQ=255)
                                                                                                                                      ccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTg                                <  1:694871/132‑1 (MQ=255)
                                                                                                                                      ccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTg                                >  2:694871/1‑132 (MQ=255)
                                                                                                                                        gCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTAt              <  1:25123/148‑1 (MQ=255)
                                                                                                                                         cccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATc             >  2:854621/5‑148 (MQ=255)
                                                                                                                                         cTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATc             >  1:259897/1‑148 (MQ=255)
                                                                                                                                          ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGtt                              <  4:178081/127‑1 (MQ=255)
                                                                                                                                          ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGtt                              >  3:178081/4‑130 (MQ=255)
                                                                                                                                          ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCgg           >  3:772046/4‑149 (MQ=255)
                                                                                                                                          ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCgg           >  4:811320/4‑149 (MQ=255)
                                                                                                                                           catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCggg          <  3:811320/147‑1 (MQ=255)
                                                                                                                                           catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCggg          >  4:856816/3‑149 (MQ=255)
                                                                                                                                                  cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgc    >  2:105348/1‑148 (MQ=255)
                                                                                                                                                   aCTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgcag  <  4:566659/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419247‑2419542

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