Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,186 C→T 6.5% A100A (GCG→GCA yfcF ← glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,1860CT 6.5% 451.0 / 14.6 138A100A (GCG→GCAyfcFglutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (59/70);  new base T (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.21e-01

Only 100 of 142 total aligned reads displayed.

CTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAACCTT  >  NC_000913/2420038‑2420334
                                                                                                                                                    |                                                                                                                                                    
cTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACCACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCAcgc                                                                                                                                                      <  4:937962/149‑1 (MQ=255)
 tGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCAcgcg                                                                                                                                                     >  2:361625/1‑149 (MQ=255)
          gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                              >  3:617407/1‑145 (MQ=255)
          gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtttt                                                                                                                                            <  2:957765/149‑1 (MQ=255)
          gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtttt                                                                                                                                            <  2:137504/149‑1 (MQ=255)
          gAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtttt                                                                                                                                            >  3:246620/1‑149 (MQ=255)
           aaCAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTc                                                                                                                                           <  1:592377/149‑1 (MQ=255)
             cAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTa                                                                                                                                         <  4:839725/149‑1 (MQ=255)
                tttCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAAt                                                                                                                                      <  4:124806/149‑1 (MQ=255)
                   ctctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCaa                                                                                                                                   >  3:190015/1‑149 (MQ=255)
                    tctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAg                                                                                                                                  >  1:8497/1‑149 (MQ=255)
                     ctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                                 <  1:385180/149‑1 (MQ=255)
                     ctGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                                 <  4:956084/149‑1 (MQ=255)
                          cTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGc                                                                                                                                 <  1:579281/144‑1 (MQ=255)
                           tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGTCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCa                                                                                                                                                   >  4:1112179/1‑125 (MQ=255)
                           tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGTCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCa                                                                                                                                                   <  3:1112179/125‑1 (MQ=255)
                           tGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTaaa                                                                                                                                       <  3:673404/137‑1 (MQ=255)
                             gCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATaatt                                                                                                                         <  4:785365/149‑3 (MQ=255)
                             gCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                         <  3:683369/149‑1 (MQ=255)
                             gCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAAt                                                                                                                         >  1:888434/1‑149 (MQ=255)
                              cTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCACCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATa                                                                                                                        <  2:1106494/149‑1 (MQ=255)
                                  tCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGtt                                                                                                                    <  2:450966/149‑1 (MQ=255)
                                  tCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAAGCAAGCGGATAAATACGtt                                                                                                                    <  2:496518/149‑1 (MQ=255)
                                    ccTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTcc                                                                                                                  >  2:384655/1‑149 (MQ=255)
                                     cTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTccc                                                                                                                 <  4:284583/149‑1 (MQ=255)
                                      tCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCa                                                                                                                >  3:510272/1‑149 (MQ=255)
                                              tAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTg                                                                                                        <  4:1004338/149‑1 (MQ=255)
                                                 tGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                              >  3:1141420/1‑106 (MQ=255)
                                                 tGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGaca                                                                                                                                              <  4:1141420/108‑3 (MQ=255)
                                                   gcgcTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCa                                                                                                   >  1:548396/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)
                                                             cGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCtt                                                                                         >  3:1051202/1‑149 (MQ=255)
                                                              gCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTc                                                                                        >  4:857419/1‑149 (MQ=255)
                                                               cccccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTcc                                                                                       <  2:8497/149‑1 (MQ=255)
                                                                ccccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCGCGTGCACGac                                                                                                                                               >  4:418046/1‑91 (MQ=255)
                                                                ccccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCa                                                                                      >  3:892727/1‑149 (MQ=255)
                                                                  ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGAt                                                                                            >  3:94787/1‑141 (MQ=255)
                                                                  ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGa                                                                                    >  4:248309/1‑149 (MQ=255)
                                                                  ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGa                                                                                    >  2:647999/1‑149 (MQ=255)
                                                                  ccGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGa                                                                                    <  2:120201/149‑1 (MQ=255)
                                                                    gCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATa                                                                                  >  2:432688/1‑149 (MQ=255)
                                                                       aaGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTc                                                                               <  4:566371/149‑1 (MQ=255)
                                                                       aaGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTc                                                                               <  3:51621/149‑1 (MQ=255)
                                                                         gACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtt                                                                                                                                              >  1:923482/1‑84 (MQ=255)
                                                                         gACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGAtt                                                                                                                                              <  2:923482/84‑1 (MQ=255)
                                                                             acaTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAAt                                                                         <  1:338013/149‑1 (MQ=255)
                                                                               aTCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATgg                                                                       >  3:724656/1‑149 (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)
                                                                                           ctctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCac                                                           <  4:510272/149‑1 (MQ=255)
                                                                                            tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                          >  2:219336/1‑149 (MQ=255)
                                                                                            tctTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACATTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCact                                                          <  3:561652/149‑1 (MQ=255)
                                                                                               tCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactcaa                                                       <  1:441228/149‑1 (MQ=255)
                                                                                                cgcgGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGt                                                                                                         <  4:706821/98‑1 (MQ=255)
                                                                                                cgcgGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCactcaac                                                      <  2:548396/149‑1 (MQ=255)
                                                                                                     aTGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCaaa                                                 <  2:653027/149‑1 (MQ=255)
                                                                                                       gggCATCAGATCGCTGCGCAGCCCGGCCTGGATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAat                                               >  4:624008/1‑149 (MQ=255)
                                                                                                       gggCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAat                                               >  2:583972/1‑149 (MQ=255)
                                                                                                          cATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAatcat                                            >  1:2134/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                                                                                 >  3:555449/1‑106 (MQ=255)
                                                                                                                aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                      >  4:485665/1‑149 (MQ=255)
                                                                                                                aTCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGAttt                                      <  3:815338/149‑1 (MQ=255)
                                                                                                                   gCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGta                                   >  4:1121915/1‑149 (MQ=255)
                                                                                                                     tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                                >  3:625813/1‑118 (MQ=255)
                                                                                                                     tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                                <  4:625813/118‑1 (MQ=255)
                                                                                                                     tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaat                                 >  1:856540/1‑149 (MQ=255)
                                                                                                                     tGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaat                                 >  3:306989/1‑149 (MQ=255)
                                                                                                                       cgcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGtaataa                               >  2:759216/1‑149 (MQ=255)
                                                                                                                        gcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTcc                                                                                       <  3:77721/92‑1 (MQ=255)
                                                                                                                        gcAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTcc                                                                                       >  4:77721/1‑92 (MQ=255)
                                                                                                                              aGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                              >  1:118623/1‑143 (MQ=255)
                                                                                                                              aGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGcaca                        >  3:1136358/1‑149 (MQ=255)
                                                                                                                                  cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                                >  2:961595/1‑105 (MQ=255)
                                                                                                                                  cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCagaaga                                                                <  1:961595/105‑1 (MQ=255)
                                                                                                                                  cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                    >  3:351941/1‑149 (MQ=255)
                                                                                                                                  cTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGc                    >  2:1023037/1‑149 (MQ=255)
                                                                                                                                    gAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgt                  <  1:692753/149‑1 (MQ=255)
                                                                                                                                     aaTCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCgtg                 >  4:943145/1‑149 (MQ=255)
                                                                                                                                         aatCGTGCGCGTGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTg             >  4:1026917/3‑149 (MQ=255)
                                                                                                                                           tCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGc                                                                      >  4:593479/1‑90 (MQ=255)
                                                                                                                                           tCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGc                                                                      <  3:593479/90‑1 (MQ=255)
                                                                                                                                            cGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGAcc          <  1:219336/149‑1 (MQ=255)
                                                                                                                                             gTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCg         >  1:1015093/1‑149 (MQ=255)
                                                                                                                                             gTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCg         >  4:136357/1‑149 (MQ=255)
                                                                                                                                               gcgcgcgCACGATTTTCTGAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCATCTCAAAATCTTCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTa       >  4:88135/4‑149 (MQ=255)
                                                                                                                                               gCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                              >  1:263184/1‑126 (MQ=255)
                                                                                                                                               gCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAAc                              <  2:263184/126‑1 (MQ=255)
                                                                                                                                                 aCGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAc     >  2:295040/1‑149 (MQ=255)
                                                                                                                                                 aCGCGCACGATTTGCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAAc     <  4:724656/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCCATCGCGAACAGTTTCTCTGCACTGGCTTTTCCCTCGGCCGTTAGTGGCGCTTTTTTCGCCCCCGCAAAGACAACATCCGTCGGACGCTCTTCGCGGATGGGCATCAGATCGCTGCGCAGCCAGGCCTGAATCTGTCGTGCACGCGCACGATTTTCTAAATCAAGCGGATAAATACGTTCCCAGGTCGGTGGCGCAAATCGATCTTCCAGATACTCCGCAATGGCAGAAGATTCACTCAACTCAAAATCATCGATTTGTAATAACGGCACACGGCGTGTCTGACCGTAACCTT  >  NC_000913/2420038‑2420334

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: