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→ATA) yicIputative 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

GTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAT  >  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)
                                    ggTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTtggt‑g                                                                                                                >  1:183485/1‑149 (MQ=255)
                                    ggTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTACGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTtggt‑g                                                                                                                <  1:244625/149‑1 (MQ=255)
                                           gAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATg                                                                                                         <  1:776002/149‑1 (MQ=255)
                                                tCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTg                                                                                                    >  1:948439/1‑149 (MQ=255)
                                                tCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTg                                                                                                    <  1:1126215/149‑1 (MQ=255)
                                                     gggTTCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAg                                                                                               <  1:141829/149‑1 (MQ=255)
                                                      ggTCCCACTCAAAGTCGCACCACTGGAAGGCTTTCGTCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGc                                                                                              <  2:564865/149‑1 (MQ=255)
                                                      ggTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGACAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGGGGTTAGc                                                                                              >  2:45964/1‑149 (MQ=255)
                                                        tCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATCAAGCTGGGGACCGTCGCTTCGTCGAAGTTGGT‑GGTAAATGAAGTGGTTAGCca                                                                                            >  1:264724/1‑149 (MQ=255)
                                                            aCTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAgg                                                                                        >  2:581294/1‑149 (MQ=255)
                                                              tCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGGGGTTAGCCACAGGcc                                                                                      >  2:510055/1‑149 (MQ=255)
                                                               cAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCg                                                                                     <  2:837418/149‑1 (MQ=255)
                                                                aaaaTCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGa                                                                                    >  1:864645/1‑149 (MQ=255)
                                                                aaaaTCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGa                                                                                    >  1:135985/1‑149 (MQ=255)
                                                                        accacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGGCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAc                                                                            >  2:236019/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAAGAGTCCAAGTGGAAAAAATGCACCGGCAGATTGCGTTCCGCCATAAAAAAGATAAAGCTGGTGAACGTAGCTAAGGGGTAGCTGGT‑GGTAAATGAATGGGTTAGCGACAGACCCAAGGACCAcgc                                                                          >  2:135521/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGATAACATGCAGCGGCAGAGTGCGGTCCGCCATACCATAGATAAAGCAGTTTACCGTCGCATCGTCGAAGTTCGC‑GGTAAATGAGGTGGTACGCAACAGGCCGAAGGACCAcgc                                                                          >  2:386837/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgc                                                                          >  1:461190/1‑149 (MQ=255)
                                                                          cacTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgc                                                                          >  1:172630/1‑149 (MQ=255)
                                                                           acTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCAcgcg                                                                         >  1:160570/1‑149 (MQ=255)
                                                                             tGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTCGTTGGA‑GGTAAATGAAGTGGGTAGCCACAGGCCGAAGGACCACgcggg                                                                       >  2:919283/1‑149 (MQ=255)
                                                                                 aGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGggcggc                                                                   <  2:1081710/149‑1 (MQ=255)
                                                                                  gtccttCATCCAGAAACCGTCAAAGTGGAAAACATGCAGCGGCAGATTGCCTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCa                                                                  <  2:1090795/145‑1 (MQ=255)
                                                                                       ccATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATGGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGcc                                                             <  2:1020218/148‑1 (MQ=255)
                                                                                           ccAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGgac                                                         >  2:331414/1‑149 (MQ=255)
                                                                                           acaaaaaCAGTCAAAGTGGAAAACATGCAGCGGCAGATTGGGTTCCGCCATACCATCGATAAAGCCGTTTACCGTCGCTTCGGCTTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGgac                                                         <  2:574767/145‑1 (MQ=255)
                                                                                              gAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGGTAGCCACAGGCCGAAGGACCCCGCGGGCGGCAGCGCCGgacgac                                                      >  2:957679/1‑149 (MQ=255)
                                                                                               aaaCAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGAcc                                                     >  1:258772/1‑149 (MQ=255)
                                                                                                  cAGTCAACGTGGAAAACATGCAGCGCCAGATTGCGTTCCGCCGTACCATCGATAAAGCTGTTTACCGACGCTTCGTCGTAGATGGT‑GGTAAATGAAGTGGGTAGACACAGGCCGCAGGACCACGAGGGCGGCAGCGCCGGACGACCAGt                                                  >  2:289180/1‑149 (MQ=255)
                                                                                                  cAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGt                                                  <  1:864706/149‑1 (MQ=255)
                                                                                                   aGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTa                                                 >  2:1120839/1‑149 (MQ=255)
                                                                                                    gTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTaa                                                >  2:939234/1‑149 (MQ=255)
                                                                                                          gTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCggg                                          <  1:379039/149‑1 (MQ=255)
                                                                                                            ggAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAATTGGTTAGCCACAGGCCGAAGGGCCACGCGGGCGGCAGCGCCGGACGACCGGTAAGGCGGGta                                        >  2:504194/1‑149 (MQ=255)
                                                                                                               aaaCATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATaa                                     >  2:1125310/1‑149 (MQ=255)
                                                                                                               aaaCATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGCCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGCAAAGCGGGTATaa                                     <  2:587380/149‑1 (MQ=255)
                                                                                                                    tGCCGCGGCAGATTGCGTTCCGCCATACCCTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATc                                <  1:212377/149‑1 (MQ=255)
                                                                                                                      cAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCaa                              <  2:135985/149‑1 (MQ=255)
                                                                                                                      cAGCGGCAGATTGCGTTCCGCCATACCATCTATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCaa                              <  2:864645/149‑1 (MQ=255)
                                                                                                                            cAGATTGCGTTCCCCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGG‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCg                        <  2:734708/149‑1 (MQ=255)
                                                                                                                                 tGCGTTCCGCCATACTATCGATAAAGCTGGTTACCGTCGCTTCGTCGTAGTGGGT‑GGTAACTGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTc                   <  1:565643/149‑1 (MQ=255)
                                                                                                                                 tGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGGCGACCAGTAGAGCGGGTATAACGATCAAGTACCGCTTTc                   >  2:486275/1‑149 (MQ=255)
                                                                                                                                   cGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGT‑TTTGGTAGGTAACTGCGAATGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGGATAACGATCAAGTACCGCTTTCgg                 >  1:558572/1‑149 (MQ=255)
                                                                                                                                    gTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGc                >  1:4551/1‑149 (MQ=255)
                                                                                                                                    gTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGc                <  1:421522/149‑1 (MQ=255)
                                                                                                                                     ttCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCg               >  1:735158/1‑149 (MQ=255)
                                                                                                                                         gCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCgg           <  2:948439/149‑1 (MQ=255)
                                                                                                                                         gCCATACCAACGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCgg           <  1:230646/149‑1 (MQ=255)
                                                                                                                                          ccATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGCCGACCAGTGAAGCGGGTATAACGATCAGGTGCCGCTTTCGGCGTTggg          >  1:67930/1‑149 (MQ=255)
                                                                                                                                           cATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGc         >  1:383404/1‑149 (MQ=255)
                                                                                                                                              aCCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGt      <  1:25554/149‑1 (MQ=255)
                                                                                                                                                 aTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTc     >  2:1120280/1‑147 (MQ=255)
                                                                                                                                                 aTCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTc     <  1:1120280/147‑1 (MQ=255)
                                                                                                                                                  tCGATAAAGCTGCTTACCGCCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGAGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAt  <  2:28480/149‑1 (MQ=255)
                                                                                                                                                  tCGATAAAGCGGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATACCGCTCAAGTACCGCTTTCGGCGTCGGGCCGTCGAt  <  2:1111073/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCCTTTCGCTTTCAGGCGGCGGATCATCCCTTCCGGGTCAGGGAAAGTCAGCGGGTCCCACTCAAAATCGCACCACTGGAAGGCTTTCATCCAGAAACAGTCAAAGTGGAAAACATGCAGCGGCAGATTGCGTTCCGCCATACCATCGATAAAGCTGTTTACCGTCGCTTCGTCGTAGTTGGT‑GGTAAATGAAGTGGTTAGCCACAGGCCGAAGGACCACGCGGGCGGCAGCGCCGGACGACCAGTAAAGCGGGTATAACGATCAAGTACCGCTTTCGGCGTCGGGCCGTCGAT  >  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: ‑