Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,419,364 A→T 5.2% F142Y (TTC→TAC)  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,3640AT 5.2% 245.4 / 3.4 77F142Y (TTC→TAC) yfcEphosphodiesterase activity on bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base T (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TTGAGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAC  >  NC_000913/2419218‑2419509
                                                                                                                                                  |                                                                                                                                                 
ttGAGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGaa                                                                                                                                                  <  4:65472/148‑1 (MQ=255)
ttGAGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGaa                                                                                                                                                  <  2:72904/148‑1 (MQ=255)
  gAGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGt                                                                                                                                                <  6:68500/148‑1 (MQ=255)
   aGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTgg                                                                                                                                              <  6:63712/149‑1 (MQ=255)
     ttGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGa                                                                                                                                             >  1:20377/1‑148 (MQ=255)
     ttGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGa                                                                                                                                             <  8:244159/148‑1 (MQ=255)
        tgGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAAt                                                                                                                                         <  1:25343/149‑1 (MQ=255)
        tgGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAAt                                                                                                                                         <  6:225715/149‑1 (MQ=255)
            tAAATTACTGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGTTTGAAGTGGAAAATCTcc                                                                                                                                     <  2:330818/149‑1 (MQ=255)
            tAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTcc                                                                                                                                     <  7:18005/149‑1 (MQ=255)
               aTTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTccccc                                                                                                                                  <  8:102361/149‑1 (MQ=255)
                  aCGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGc                                                                                                                               >  8:65239/1‑149 (MQ=255)
                   cGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGc                                                                                                                               <  4:2116/148‑1 (MQ=255)
                          aTCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGc                                                                                                                       <  8:285135/149‑1 (MQ=255)
                                 ccTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTgg                                                                                                                 >  4:266765/1‑148 (MQ=255)
                                    gcgcAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTgg                                                                                                                 >  4:149604/1‑145 (MQ=255)
                                          tGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATggg                                                                                                         <  1:80598/147‑1 (MQ=255)
                                          tGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATggg                                                                                                         >  2:80598/1‑147 (MQ=255)
                                               cTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTAACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCAACTGGCAGATGGGCATGAcc                                                                                                  >  6:17482/1‑149 (MQ=255)
                                                tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGTAGTGGAAAATCTCCCCCCGCTGTTCt                                                                                                                         <  6:43126/125‑1 (MQ=255)
                                                tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGTAGTGGAAAATCTCCCCCCGCTGTTCt                                                                                                                         >  5:43126/1‑125 (MQ=255)
                                                tttGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCg                                                                                                 >  2:392394/1‑149 (MQ=255)
                                                  tGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCAGCGAGCCGGGGTTGAAGTGGAAAATCTCCACCAGCTGTTCTGCCACTGGCAGATGGGTATGACCGTa                                                                                               >  8:334631/1‑149 (MQ=255)
                                                   gATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGTAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTa                                                                                               >  5:133337/1‑148 (MQ=255)
                                                   gATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTa                                                                                               <  2:244458/148‑1 (MQ=255)
                                                   gATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTa                                                                                               >  2:325190/1‑148 (MQ=255)
                                                    aTCATTGAGTGCGATTACCCTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCCCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTaca                                                                                             >  1:286473/1‑149 (MQ=255)
                                                    aTCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTaca                                                                                             <  2:119191/149‑1 (MQ=255)
                                                     tCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTacac                                                                                            >  8:62604/1‑149 (MQ=255)
                                                     tCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTacac                                                                                            >  3:274009/1‑149 (MQ=255)
                                                       aTTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACAcc                                                                                           <  3:266765/148‑1 (MQ=255)
                                                       aTTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCa                                                                                          <  5:39373/149‑1 (MQ=255)
                                                          gtgTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGACTACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGc                                                                                        <  5:39423/146‑1 (MQ=255)
                                                          gAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGca                                                                                       >  7:332120/1‑149 (MQ=255)
                                                                 aTTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGTAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGtt                                                                                <  6:133337/149‑1 (MQ=255)
                                                                 aTTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGtt                                                                                >  3:7349/1‑149 (MQ=255)
                                                                  ttACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTc                                                                               >  7:385494/1‑149 (MQ=255)
                                                                     cacTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCa                                                                                                               <  6:62965/114‑1 (MQ=255)
                                                                     cacTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCa                                                                                                               >  5:62965/1‑114 (MQ=255)
                                                                     cacTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGGATGACCGTACACCAGCACAACGTTCTgg                                                                            >  2:190277/1‑149 (MQ=255)
                                                                      acTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTg                                                                             >  4:146937/1‑147 (MQ=255)
                                                                      acTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTg                                                                             <  3:146937/147‑1 (MQ=255)
                                                                      acTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGt                                                                           >  2:364191/1‑149 (MQ=255)
                                                                        tAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGtt                                                                          >  5:285197/1‑148 (MQ=255)
                                                                          aGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTaa                                                                       <  8:338690/149‑1 (MQ=255)
                                                                              cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgca                                                                   <  8:385494/149‑1 (MQ=255)
                                                                              cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgc                                                                    <  7:65239/148‑1 (MQ=255)
                                                                              cGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgc                                                                    <  1:392394/148‑1 (MQ=255)
                                                                               gTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAgcag                                                                  >  6:305206/1‑149 (MQ=255)
                                                                                tCATTATCCAGCATGCCATAACTCGCCGGATTAGCGCCTTTCGGAATTCTCTCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTAGACCAGCACATCGTTCTGGTTTAAAgcagg                                                                 <  2:406640/149‑1 (MQ=255)
                                                                                         aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                              >  1:89256/1‑127 (MQ=255)
                                                                                         aGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCt                                                                              <  2:89256/127‑1 (MQ=255)
                                                                                               ccATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCTCCCGCTGTTCTGCCACTGGCGGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGc                                                   <  4:70260/148‑1 (MQ=255)
                                                                                                cATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGcc                                                  >  8:118772/1‑148 (MQ=255)
                                                                                                      tCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCaaa                                               >  2:140840/1‑145 (MQ=255)
                                                                                                        gCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGa                                           <  1:140840/147‑1 (MQ=255)
                                                                                                        gCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGAtt                                         <  1:325190/149‑2 (MQ=255)
                                                                                                        gCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGAt                                          <  3:149604/148‑1 (MQ=255)
                                                                                                          cGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGAc                                       >  7:285647/1‑149 (MQ=255)
                                                                                                             aTTACCGCCTTTCGGAATTCTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCAt                                    <  8:285647/149‑1 (MQ=255)
                                                                                                             aTTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCAt                                    <  6:285197/149‑1 (MQ=255)
                                                                                                                aCCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAAGCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGc                                  <  5:295258/148‑1 (MQ=255)
                                                                                                                 ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGt                                >  2:57440/1‑149 (MQ=255)
                                                                                                                 ccGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGt                                >  4:160008/1‑149 (MQ=255)
                                                                                                                   gCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTc                               >  8:342778/1‑148 (MQ=255)
                                                                                                                    ccTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCa                              >  8:315342/1‑148 (MQ=255)
                                                                                                                      tttCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCaaa                            >  8:277077/1‑148 (MQ=255)
                                                                                                                           gAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGa                      <  5:305206/149‑1 (MQ=255)
                                                                                                                                cTCACCGAGCCGTGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgtt                  <  7:316857/148‑1 (MQ=255)
                                                                                                                                cTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgttg                 <  5:237896/149‑1 (MQ=255)
                                                                                                                                cTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgtt                  <  7:315342/148‑1 (MQ=255)
                                                                                                                                cTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACgtt                  <  2:20377/148‑1 (MQ=255)
                                                                                                                                   aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGTGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt              >  3:19358/1‑149 (MQ=255)
                                                                                                                                   aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt              <  7:277077/149‑1 (MQ=255)
                                                                                                                                   aCCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGttt              >  1:377557/1‑149 (MQ=255)
                                                                                                                                       gcccGGGTTGGACGGGGAAAACCCCCCCCCGCTGTTCAGCCACTGGCAGATGGGTATGACCGTACACCAGCACATAGATCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTcc          <  6:314453/147‑1 (MQ=255)
                                                                                                                                        gCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCa         <  8:332120/149‑1 (MQ=255)
                                                                                                                                           ggggTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTGTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGt       >  3:107112/1‑148 (MQ=255)
                                                                                                                                               ttGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAc  <  4:274009/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGAGTTGTGGGTAAATTACGGATTAATCGCGACCTGCGCAATGATACTTTGATCATTGAGTGCGATTACACTAAGTACGTCATTATCCAGCATGCCATAACTCGCCGGATTACCGCCTTTCGGAATACTCACCGAGCCGGGGTTGAAGTGGAAAATCTCCCCCCGCTGTTCTGCCACTGGCAGATGGGTATGACCGTACACCAGCACATCGTTCTGGTTTAAAGCAGGCAGATTTTCCGGGCCAAAAAGATGACCATGCGTCAAAAACAGACGTTGTTTTTCCAGTAATAC  >  NC_000913/2419218‑2419509

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: