Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,419,395 C→G 7.0% 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.0% 455.5 / 21.1 129V132L (GTG→CTG) yfcEphosphodiesterase activity on bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (56/64);  new base G (6/3);  total (62/67)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-02

Only 100 of 129 total aligned reads displayed.

GACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419249‑2419542
                                                                                                                                                  |                                                                                                                                                   
gACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTg                                                                                                                                                   <  1:813477/149‑1 (MQ=255)
     gcgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGa                                                                                                                                              <  2:362055/149‑1 (MQ=255)
     gcgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGa                                                                                                                                              >  4:906717/1‑149 (MQ=255)
      cgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGAt                                                                                                                                             >  2:388977/1‑149 (MQ=255)
        cAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  3:711684/149‑1 (MQ=255)
        cAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATgg                                                                                                                                           <  3:358295/149‑1 (MQ=255)
         aaTGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATggg                                                                                                                                          >  1:96416/1‑149 (MQ=255)
          aTGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGt                                                                                                                                         <  1:429300/149‑1 (MQ=255)
             aTACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATggg                                                                                                                                          >  4:911252/1‑145 (MQ=255)
             aTACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATg                                                                                                                                      <  1:642507/149‑1 (MQ=255)
              tACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGa                                                                                                                                     <  3:273668/149‑1 (MQ=255)
                cTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGTGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGAcc                                                                                                                                   <  1:602309/149‑1 (MQ=255)
                 tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCg                                                                                                                                  >  4:583353/1‑149 (MQ=255)
                      tCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTacac                                                                                                                             >  1:641112/1‑149 (MQ=255)
                      tCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTacac                                                                                                                             >  4:65184/1‑149 (MQ=255)
                         ttGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAg                                                                                                                          >  2:335736/1‑149 (MQ=255)
                          tGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGc                                                                                                                         <  4:837655/149‑1 (MQ=255)
                            aGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGcac                                                                                                                       >  3:365868/1‑149 (MQ=255)
                               gCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATc                                                                                                                    <  3:421714/149‑1 (MQ=255)
                                  aTTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGtt                                                                                                                 <  1:660481/149‑1 (MQ=255)
                                   ttACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTc                                                                                                                <  3:235194/149‑1 (MQ=255)
                                    tACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                                                               <  3:4287/149‑1 (MQ=255)
                                      cacTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTgg                                                                                                             <  3:583353/149‑1 (MQ=255)
                                       acTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGt                                                                                                            >  2:332728/1‑149 (MQ=255)
                                       acTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGt                                                                                                            >  2:55920/1‑149 (MQ=255)
                                          aaGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTgg                                                                                                             <  3:904829/145‑1 (MQ=255)
                                           aGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaa                                                                                                        <  1:645764/149‑1 (MQ=255)
                                            gTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaaa                                                                                                       >  3:355714/1‑149 (MQ=255)
                                             tACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAg                                                                                                      >  1:52163/1‑149 (MQ=255)
                                               cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgca                                                                                                    <  4:896149/149‑1 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  2:661084/1‑149 (MQ=255)
                                                gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  4:558033/1‑149 (MQ=255)
                                                gTCATAATCCAGCATGCCATAACTCGCCGGATTGCCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                                                   >  1:266695/1‑149 (MQ=255)
                                                 tCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcagg                                                                                                  <  1:335736/149‑1 (MQ=255)
                                                   aTTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcaggca                                                                                                >  3:111739/1‑149 (MQ=255)
                                                    ttATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcaggcag                                                                                               >  2:431862/1‑149 (MQ=255)
                                                     tATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGa                                                                                              <  1:55920/149‑1 (MQ=255)
                                                      aTCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGCTCTGGTTTAAAGCAGGCAGAt                                                                                             >  3:444242/1‑149 (MQ=255)
                                                       tCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAtt                                                                                            >  1:50100/1‑149 (MQ=255)
                                                         cAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGAtttt                                                                                          >  1:703230/1‑149 (MQ=255)
                                                          aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         >  1:401626/1‑149 (MQ=255)
                                                          aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         <  4:227033/149‑1 (MQ=255)
                                                          aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         <  3:246578/149‑1 (MQ=255)
                                                          aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTc                                                                                         <  1:661221/149‑1 (MQ=255)
                                                            cATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCg                                                                                       >  1:568144/1‑149 (MQ=255)
                                                             aTGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCgg                                                                                      <  4:397071/149‑1 (MQ=255)
                                                             aTGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCgg                                                                                      <  4:656055/149‑1 (MQ=255)
                                                               gCCCTAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGc                                                                                    <  3:380076/149‑1 (MQ=255)
                                                                ccATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGcc                                                                                   >  1:168365/1‑149 (MQ=255)
                                                                     aCTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaaaa                                                                              <  4:16043/149‑1 (MQ=255)
                                                                       tCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                                                            <  4:641537/149‑1 (MQ=255)
                                                                        cGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGAt                                                                           <  4:381027/149‑1 (MQ=255)
                                                                         gCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATg                                                                          >  1:470549/1‑149 (MQ=255)
                                                                          ccGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGa                                                                         <  2:96416/149‑1 (MQ=255)
                                                                            ggATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAcc                                                                       >  1:231427/1‑149 (MQ=255)
                                                                             gATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCa                                                                      >  3:138030/1‑149 (MQ=255)
                                                                             gATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCa                                                                      >  1:935037/1‑149 (MQ=255)
                                                                               ttACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATg                                                                    >  4:158063/1‑149 (MQ=255)
                                                                                tACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGc                                                                   <  2:26433/149‑1 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                                                            >  2:764999/1‑138 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                                                            <  1:764999/138‑1 (MQ=255)
                                                                                  ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGt                                                                 >  4:28794/1‑149 (MQ=255)
                                                                                   cGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTc                                                                >  3:438993/1‑149 (MQ=255)
                                                                                   cGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTc                                                                >  1:596455/1‑149 (MQ=255)
                                                                                            gAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGa                                                       <  2:856564/149‑1 (MQ=255)
                                                                                                 cTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttg                                                  >  1:929568/1‑149 (MQ=255)
                                                                                                  tCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgt                                                 >  3:192626/1‑149 (MQ=255)
                                                                                                  tCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgt                                                 <  3:152275/149‑1 (MQ=255)
                                                                                                  tCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttgt                                                 >  2:450313/1‑149 (MQ=255)
                                                                                                    aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGTTGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAATGATGACCATGCGTCAAAAACAGACGTTGttt                                               <  4:566178/149‑1 (MQ=255)
                                                                                                     ccGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGtttt                                              >  3:550520/1‑149 (MQ=255)
                                                                                                        aGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTcc                                           <  2:641600/149‑1 (MQ=255)
                                                                                                         gCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa                                          <  1:711270/149‑1 (MQ=255)
                                                                                                         gCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa                                          <  2:568144/149‑1 (MQ=255)
                                                                                                             gggTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTaa                                      <  2:875779/149‑1 (MQ=255)
                                                                                                                     gTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGtt                              <  4:355714/149‑1 (MQ=255)
                                                                                                                      tGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGAc                                                      <  4:917060/124‑1 (MQ=255)
                                                                                                                      tGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGAc                                                      >  3:917060/1‑124 (MQ=255)
                                                                                                                           aaTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACg                        >  3:675999/1‑149 (MQ=255)
                                                                                                                            aTCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACgg                       <  2:703230/149‑1 (MQ=255)
                                                                                                                              ctcCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGgcg                     >  1:336359/1‑149 (MQ=255)
                                                                                                                               tcCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGgcgc                    <  3:158063/149‑1 (MQ=255)
                                                                                                                                ccccccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGcgcg                   <  2:336359/149‑1 (MQ=255)
                                                                                                                                   cccGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGt                               >  3:99036/1‑134 (MQ=255)
                                                                                                                                     cGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTAt              <  1:388977/149‑1 (MQ=255)
                                                                                                                                       cTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCg            <  4:365868/149‑1 (MQ=255)
                                                                                                                                       cTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCg            <  2:747313/149‑1 (MQ=255)
                                                                                                                                        ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCgg           >  4:613809/4‑149 (MQ=255)
                                                                                                                                        ccatCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCgg           >  1:441632/4‑149 (MQ=255)
                                                                                                                                         gTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCggg          <  1:151276/149‑1 (MQ=255)
                                                                                                                                         catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGc                            >  1:848083/3‑131 (MQ=255)
                                                                                                                                         catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGc                            <  2:848083/129‑1 (MQ=255)
                                                                                                                                         catCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCggg          >  4:790656/3‑149 (MQ=255)
                                                                                                                                          atCTGCCAGTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCg            <  2:362403/145‑1 (MQ=255)
                                                                                                                                             tGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGaaaa      >  3:80223/1‑148 (MQ=255)
                                                                                                                                              gCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATg     <  2:641112/149‑1 (MQ=255)
                                                                                                                                               ccACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgc    <  4:438993/149‑1 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgca   >  2:932318/1‑149 (MQ=255)
                                                                                                                                                cACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgca   <  1:332728/149‑1 (MQ=255)
                                                                                                                                                 aCTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATgcag  >  2:75223/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATACCTGTTGCCACGGCGCGGTTATCGGGAAATGCAG  >  NC_000913/2419249‑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: