Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,183 A→G 6.5% R101R (CGT→CGC)  yfcF ← glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,1830AG 6.5% 449.5 / 13.8 138R101R (CGT→CGC) yfcFglutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (59/70);  new base G (5/4);  total (64/74)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 141 total aligned reads displayed.

GTTCTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAC  >  NC_000913/2420035‑2420331
                                                                                                                                                    |                                                                                                                                                    
gTTCTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCa                                                                                                                                                      >  3:956084/1‑149 (MQ=255)
   cTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACCACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCAcgc                                                                                                                                                   <  4:937962/149‑1 (MQ=255)
             gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                           <  4:617407/147‑3 (MQ=255)
             gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                           >  3:617407/1‑145 (MQ=255)
             gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtttt                                                                                                                                         <  2:957765/149‑1 (MQ=255)
             gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtttt                                                                                                                                         <  2:137504/149‑1 (MQ=255)
              aaCAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTc                                                                                                                                        <  1:592377/149‑1 (MQ=255)
                cAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTa                                                                                                                                      <  4:839725/149‑1 (MQ=255)
                 aGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGacaga                                                                                                                                         <  2:898370/145‑5 (MQ=255)
                  gTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTaaa                                                                                                                                    <  4:252690/149‑1 (MQ=255)
                   tttCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAAt                                                                                                                                   >  3:1147453/1‑149 (MQ=255)
                      ctctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCaa                                                                                                                                >  3:190015/1‑149 (MQ=255)
                       tctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAg                                                                                                                               >  3:785365/1‑149 (MQ=255)
                        ctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                              <  4:956084/149‑1 (MQ=255)
                        ctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                              <  1:385180/149‑1 (MQ=255)
                             cTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGGTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAATTTCTAAATCAAGc                                                                                                                              >  2:579281/1‑144 (MQ=255)
                             cTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                              <  1:579281/144‑1 (MQ=255)
                              tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGTCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCa                                                                                                                                                >  4:1112179/1‑125 (MQ=255)
                              tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGTCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCa                                                                                                                                                <  3:1112179/125‑1 (MQ=255)
                              tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTaaa                                                                                                                                    >  4:673404/1‑137 (MQ=255)
                              tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaa                                                                                                                        >  3:566371/1‑149 (MQ=255)
                                gCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                      >  1:888434/1‑149 (MQ=255)
                                gCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                      <  3:683369/149‑1 (MQ=255)
                                 cTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCACCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATa                                                                                                                     <  2:1106494/149‑1 (MQ=255)
                                  ttttCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATAc                                                                                                                    <  2:297711/149‑1 (MQ=255)
                                     tCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGtt                                                                                                                 <  2:450966/149‑1 (MQ=255)
                                       ccTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTcc                                                                                                               >  2:384655/1‑149 (MQ=255)
                                        cTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTccc                                                                                                              <  4:284583/149‑1 (MQ=255)
                                         tCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCa                                                                                                             >  2:492665/1‑149 (MQ=255)
                                                 tAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTg                                                                                                     <  4:1004338/149‑1 (MQ=255)
                                                  aGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTgg                                                                                                    <  4:1135700/149‑1 (MQ=255)
                                                    tGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                           >  3:1141420/1‑106 (MQ=255)
                                                    tGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                           <  4:1141420/108‑3 (MQ=255)
                                                      gcgcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCa                                                                                                >  2:429351/1‑149 (MQ=255)
                                                        gcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCaaa                                                                                              >  4:815338/1‑149 (MQ=255)
                                                        gcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCaaa                                                                                              >  1:120201/1‑149 (MQ=255)
                                                        gcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCaaa                                                                                              >  2:443586/1‑149 (MQ=255)
                                                         cTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAAt                                                                                             <  4:9577/149‑1 (MQ=255)
                                                            ttttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGa                                                                                          <  4:670782/149‑1 (MQ=255)
                                                             tttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGTACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGAt                                                                                         <  1:247162/149‑1 (MQ=255)
                                                             tttCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGAt                                                                                         <  4:246620/149‑1 (MQ=255)
                                                                 gCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTc                                                                                     >  4:857419/1‑149 (MQ=255)
                                                                   ccccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGac                                                                                                                                            <  3:418046/92‑2 (MQ=255)
                                                                     ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGAt                                                                                         >  3:94787/1‑141 (MQ=255)
                                                                     ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGa                                                                                 >  2:647999/1‑149 (MQ=255)
                                                                          aaGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTc                                                                            <  3:51621/149‑1 (MQ=255)
                                                                          aaGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTc                                                                            <  4:566371/149‑1 (MQ=255)
                                                                            gACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtt                                                                                                                                           >  1:923482/1‑84 (MQ=255)
                                                                                acaTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGTACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCGGCAAt                                                                      <  4:726972/149‑1 (MQ=255)
                                                                                acaTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAAt                                                                      <  1:338013/149‑1 (MQ=255)
                                                                                   tCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGc                                                                   <  1:432688/149‑1 (MQ=255)
                                                                                      gTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCaga                                                                <  1:647999/149‑1 (MQ=255)
                                                                                         ggACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                             <  3:248309/149‑1 (MQ=255)
                                                                                          gACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATgg                                                                    >  4:443076/1‑141 (MQ=255)
                                                                                          gACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATgg                                                                    <  3:443076/141‑1 (MQ=255)
                                                                                           aCGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGAtt                                                           <  3:528329/149‑1 (MQ=255)
                                                                                              ctctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCac                                                        <  4:510272/149‑1 (MQ=255)
                                                                                               tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                       >  2:219336/1‑149 (MQ=255)
                                                                                               tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                       <  1:492665/149‑1 (MQ=255)
                                                                                               tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACATTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                       <  3:561652/149‑1 (MQ=255)
                                                                                                 ttCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCTAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactca                                                     >  2:717220/1‑149 (MQ=255)
                                                                                                  tCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactcaa                                                    <  4:424280/149‑1 (MQ=255)
                                                                                                  tCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactcaa                                                    <  1:441228/149‑1 (MQ=255)
                                                                                                   cgcgGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGt                                                                                                      <  4:706821/98‑1 (MQ=255)
                                                                                                   cgcgGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGt                                                                                                      >  3:706821/1‑98 (MQ=255)
                                                                                                   cgcgGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactcaac                                                   <  2:548396/149‑1 (MQ=255)
                                                                                                        aTGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCaaa                                              <  2:653027/149‑1 (MQ=255)
                                                                                                          gggCATCAGATCGCTGCGCAGCCCGGCCTGGATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAat                                            >  4:624008/1‑149 (MQ=255)
                                                                                                             cATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcat                                         >  1:2134/1‑149 (MQ=255)
                                                                                                             cATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcat                                         >  3:302594/1‑149 (MQ=255)
                                                                                                              aTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcatc                                        <  2:888434/149‑1 (MQ=255)
                                                                                                              aTCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcatc                                        <  1:384655/149‑1 (MQ=255)
                                                                                                                   aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGTGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                   <  2:775168/149‑1 (MQ=255)
                                                                                                                   aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATAc                                                                              <  4:555449/106‑1 (MQ=255)
                                                                                                                   aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                   <  3:815338/149‑1 (MQ=255)
                                                                                                                   aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                   >  1:785844/1‑149 (MQ=255)
                                                                                                                   aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                   <  3:1120013/149‑1 (MQ=255)
                                                                                                                        tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                             >  3:625813/1‑118 (MQ=255)
                                                                                                                        tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaat                              >  1:856540/1‑149 (MQ=255)
                                                                                                                           gcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTcc                                                                                    <  3:77721/92‑1 (MQ=255)
                                                                                                                           gcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTcc                                                                                    >  4:77721/1‑92 (MQ=255)
                                                                                                                                 aGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                           <  2:118623/143‑1 (MQ=255)
                                                                                                                                 aGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                           >  1:118623/1‑143 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                             <  1:961595/105‑1 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                             >  2:961595/1‑105 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                 <  4:190015/149‑1 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                 >  2:971583/1‑149 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                 >  2:1023037/1‑149 (MQ=255)
                                                                                                                                     cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                 >  3:351941/1‑149 (MQ=255)
                                                                                                                                       gAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgt               <  3:624008/149‑1 (MQ=255)
                                                                                                                                       gAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgt               <  1:692753/149‑1 (MQ=255)
                                                                                                                                        aaTCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgtg              >  4:943145/1‑149 (MQ=255)
                                                                                                                                            aatCGTGCGCGTGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTg          >  4:1026917/3‑149 (MQ=255)
                                                                                                                                               cGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGAcc       <  1:219336/149‑1 (MQ=255)
                                                                                                                                                gTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCg      >  4:136357/1‑149 (MQ=255)
                                                                                                                                                gTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCg      >  1:1015093/1‑149 (MQ=255)
                                                                                                                                                  gCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                           >  1:263184/1‑126 (MQ=255)
                                                                                                                                                  gCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                           <  2:263184/126‑1 (MQ=255)
                                                                                                                                                   cACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTaa   <  2:785844/149‑1 (MQ=255)
                                                                                                                                                    aCGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAc  >  2:295040/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTCTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAC  >  NC_000913/2420035‑2420331

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