Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,833,667 G→T 5.6% I290I (ATC→ATA yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,833,6670GT 5.6% 199.8 / 2.9 72I290I (ATC→ATAyicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base T (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAT  >  NC_000913/3833519‑3833813
                                                                                                                                                    |                                                                                                                                                   
gTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCCGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCg                                                                                                                                                     <  2:434643/149‑1 (MQ=255)
gTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCg                                                                                                                                                     <  2:394335/149‑1 (MQ=255)
gTCCTTTCGCTTTCAGGCGGCGGATCATACCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGGGGCGGGAGATTGCGTTCCGCCATACCATCg                                                                                                                                                     >  1:498966/1‑149 (MQ=255)
   cTTTCGCTTTCAGGCGGCGGATCATGCCTTCCGGGTCAGGCAAAGTCAGCGGGTCCCACTCAAAATCGCCCCACCGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCCGCGGCAGATTGCGTTCCGCCAAACCATCGAtt                                                                                                                                                  >  1:578949/1‑148 (MQ=255)
     ttCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTTCGTTCCGCCATACCATCGATaaa                                                                                                                                                <  1:755299/149‑1 (MQ=255)
     tccGCTTTCAGGGCGCGGATCATACCTTCCGGGTCAGCGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAGCAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATaaa                                                                                                                                                <  2:722495/147‑1 (MQ=255)
        gCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCGTCCGGAAACAGTCAAAGTGGAAAACATGCAGCTGCAGATTGCGTTCCGCCATACCATCGATAAAGCt                                                                                                                                             <  1:220880/149‑1 (MQ=255)
        gCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCt                                                                                                                                             >  1:939553/1‑149 (MQ=255)
         cTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTg                                                                                                                                            >  1:460634/1‑149 (MQ=255)
          tttCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGt                                                                                                                                           <  1:918718/149‑1 (MQ=255)
          tttCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACACTCGATAACGCTGt                                                                                                                                           >  2:222229/1‑149 (MQ=255)
               ggcggcCGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAACTCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTAcc                                                                                                                                      <  2:349491/149‑1 (MQ=255)
                  ggcggATCATCCCTTCCGGGTCAGGGAAAGTCCACGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTc                                                                                                                                   <  2:518574/149‑1 (MQ=255)
                      gATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGGTTACCGTCGCtt                                                                                                                               >  2:915364/1‑149 (MQ=255)
                       atcatcCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGTTGAGGCTTTTTACCGTCGCTtc                                                                                                                              >  1:6288/1‑149 (MQ=255)
                       atcatcCCTTCCGGGTCAGGGAAAGCCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTtc                                                                                                                              <  1:424356/149‑1 (MQ=255)
                             ccTTCCGGGTCAGGGCAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTa                                                                                                                        <  2:404925/149‑1 (MQ=255)
                              cgtCCGGGGCCGGGAAAGTCAGCGGGTCCCCCCCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAg                                                                                                                       <  2:38344/147‑1 (MQ=255)
                               ttCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAACTGTTTAACGTCGCTTCGCcgtcgt                                                                                                                      >  1:982316/1‑149 (MQ=255)
                                  cGGGTCAGGGAAAGTCAGCGGGTCCCCCTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTtgg                                                                                                                   <  1:613890/149‑1 (MQ=255)
                                    ggTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTtggtg                                                                                                                >  1:183485/1‑149 (MQ=255)
                                    ggTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTACGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTtggtg                                                                                                                <  1:244625/149‑1 (MQ=255)
                                           gAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATg                                                                                                         <  1:776002/149‑1 (MQ=255)
                                                tCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTg                                                                                                    >  1:948439/1‑149 (MQ=255)
                                                tCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTg                                                                                                    <  1:1126215/149‑1 (MQ=255)
                                                     gggTTCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAg                                                                                               <  1:141829/149‑1 (MQ=255)
                                                      ggTCCCACTCAAAGTCGCACCACTGGAAGGCTTTCGTCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGc                                                                                              <  2:564865/149‑1 (MQ=255)
                                                      ggTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGACAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGGGGTTAGc                                                                                              >  2:45964/1‑149 (MQ=255)
                                                        tCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATCAAGCTGGGGACCGTCGCTTCGTCGAAGTTGGTGGTAAATGAAGTGGTTAGCca                                                                                            >  1:264724/1‑149 (MQ=255)
                                                            aCTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAgg                                                                                        >  2:581294/1‑149 (MQ=255)
                                                              tCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGGGGTTAGCCACAGGcc                                                                                      >  2:510055/1‑149 (MQ=255)
                                                               cAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCg                                                                                     <  2:837418/149‑1 (MQ=255)
                                                                aaaaTCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGa                                                                                    >  1:864645/1‑149 (MQ=255)
                                                                aaaaTCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGa                                                                                    >  1:135985/1‑149 (MQ=255)
                                                                        accacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGGCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAc                                                                            >  2:236019/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAAGAGTCCAAGTGGAAAAAATGCACCGGCAGATTGCGTTCCGCCATAAAAAAGATAAAGCTGGTGAACGTAGCTAAGGGGTAGCTGGTGGTAAATGAATGGGTTAGCGACAGACCCAAGGACCAcgc                                                                          >  2:135521/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGATAACATGCAGCGGCAGAGTGCGGTCCGCCATACCATAGATAAAGCAGTTTACCGTCGCATCGTCGAAGTTCGCGGTAAATGAGGTGGTACGCAACAGGCCGAAGGACCAcgc                                                                          >  2:386837/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgc                                                                          >  1:461190/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgc                                                                          >  1:172630/1‑149 (MQ=255)
                                                                           acTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgcg                                                                         >  1:160570/1‑149 (MQ=255)
                                                                             tGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTCGTTGGAGGTAAATGAAGTGGGTAGCCACAGGCCGAAGGACCACgcggg                                                                       >  2:919283/1‑149 (MQ=255)
                                                                                 aGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGggcggc                                                                   <  2:1081710/149‑1 (MQ=255)
                                                                                  gtccttCATCCAGAAACCGTCAAAGTGGAAAACATGCAGCGGCAGATTGCCTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCa                                                                  <  2:1090795/145‑1 (MQ=255)
                                                                                       ccATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATGGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGcc                                                             <  2:1020218/148‑1 (MQ=255)
                                                                                           ccAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGgac                                                         >  2:331414/1‑149 (MQ=255)
                                                                                           acaaaaaCAGTCAAAGTGGAAAACATGCAGCGGCAGATTGGGTTCCGCCATACCATCGATAAAGCCGTTTACCGTCGCTTCGGCTTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGgac                                                         <  2:574767/145‑1 (MQ=255)
                                                                                              gAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGGTAGCCACAGGCCGAAGGACCCCGCGGGCGGCAGCGCCGgacgac                                                      >  2:957679/1‑149 (MQ=255)
                                                                                               aaaCAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGAcc                                                     >  1:258772/1‑149 (MQ=255)
                                                                                                  cAGTCAACGTGGAAAACATGCAGCGCCAGATTGCGTTCCGCCGTACCATCGATAAAGCTGTTTACCGACGCTTCGTCGTAGATGGTGGTAAATGAAGTGGGTAGACACAGGCCGCAGGACCACGAGGGCGGCAGCGCCGGACGACCAGt                                                  >  2:289180/1‑149 (MQ=255)
                                                                                                  cAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGt                                                  <  1:864706/149‑1 (MQ=255)
                                                                                                   aGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTa                                                 >  2:1120839/1‑149 (MQ=255)
                                                                                                    gTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTaa                                                >  2:939234/1‑149 (MQ=255)
                                                                                                          gTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCggg                                          <  1:379039/149‑1 (MQ=255)
                                                                                                            ggAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAATTGGTTAGCCACAGGCCGAAGGGCCACGCGGGCGGCAGCGCCGGACGACCGGTAAGGCGGGta                                        >  2:504194/1‑149 (MQ=255)
                                                                                                               aaaCATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATaa                                     >  2:1125310/1‑149 (MQ=255)
                                                                                                               aaaCATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGCCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGCAAAGCGGGTATaa                                     <  2:587380/149‑1 (MQ=255)
                                                                                                                    tGCCGCGGCAGATTGCGTTCCGCCATACCCTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATc                                <  1:212377/149‑1 (MQ=255)
                                                                                                                      cAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCaa                              <  2:135985/149‑1 (MQ=255)
                                                                                                                      cAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCaa                              <  2:864645/149‑1 (MQ=255)
                                                                                                                            cAGATTGCGTTCCCCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGGGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCg                        <  2:734708/149‑1 (MQ=255)
                                                                                                                                 tGCGTTCCGCCATACTATCGATAAAGCTGGTTACCGTCGCTTCGTCGTAGTGGGTGGTAACTGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTc                   <  1:565643/149‑1 (MQ=255)
                                                                                                                                 tGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGGCGACCAGTAGAGCGGGTATAACGATCAAGTACCGCTTTc                   >  2:486275/1‑149 (MQ=255)
                                                                                                                                   cGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTTTTGGTAGGTAACTGCGAATGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGGATAACGATCAAGTACCGCTTTCgg                 >  1:558572/1‑149 (MQ=255)
                                                                                                                                    gTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGc                >  1:4551/1‑149 (MQ=255)
                                                                                                                                    gTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGc                <  1:421522/149‑1 (MQ=255)
                                                                                                                                     ttCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCg               >  1:735158/1‑149 (MQ=255)
                                                                                                                                         gCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCgg           <  2:948439/149‑1 (MQ=255)
                                                                                                                                         gCCATACCAACGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCgg           <  1:230646/149‑1 (MQ=255)
                                                                                                                                          ccATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGCCGACCAGTGAAGCGGGTATAACGATCAGGTGCCGCTTTCGGCGTTggg          >  1:67930/1‑149 (MQ=255)
                                                                                                                                           cATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGc         >  1:383404/1‑149 (MQ=255)
                                                                                                                                              aCCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGt      <  1:25554/149‑1 (MQ=255)
                                                                                                                                                 aTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTc     >  2:1120280/1‑147 (MQ=255)
                                                                                                                                                 aTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTc     <  1:1120280/147‑1 (MQ=255)
                                                                                                                                                  tCGATAAAGCTGCTTACCGCCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGAGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAt  <  2:28480/149‑1 (MQ=255)
                                                                                                                                                  tCGATAAAGCGGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATACCGCTCAAGTACCGCTTTCGGCGTCGGGCCGTCGAt  <  2:1111073/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGTGGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAT  >  NC_000913/3833519‑3833813

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: