Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,272,426 T→A 7.9% intergenic (+77/‑81) chaB → / → chaC cation transport regulator/cation transport regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,4260TA 7.9% 199.5 / 3.5 71intergenic (+77/‑81)chaB/chaCcation transport regulator/cation transport regulator
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base A (3/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAC  >  NC_000913/1272282‑1272573
                                                                                                                                                |                                                                                                                                                   
tGGCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGttt                                                                                                                                                <  1:858437/150‑1 (MQ=255)
tGGCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGAGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGttt                                                                                                                                                >  2:671047/1‑150 (MQ=255)
 ggCCTGGGCTGCTGTGAAGCATGAATATGCCAAATTTGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTc                                                                                                                                               <  1:1073515/150‑1 (MQ=255)
 ggCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTc                                                                                                                                               <  2:925430/150‑1 (MQ=255)
  gCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCCAGTGGTCCGTGGATTGCATATTTTCCCGTTAGAGGTTTCa                                                                                                                                              >  1:1219895/1‑150 (MQ=255)
  gCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTAGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGAGGTTTCa                                                                                                                                              >  1:1214239/1‑150 (MQ=255)
   ccTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGACTTGCAAGAGGTCCGTGGATTGCATAATGTACAGTTAGAGGTTTCaa                                                                                                                                             >  2:362833/1‑150 (MQ=255)
    cTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCaaa                                                                                                                                            >  2:710740/1‑150 (MQ=255)
     tGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCaaaa                                                                                                                                           <  1:1176296/150‑1 (MQ=255)
      gggCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAAt                                                                                                                                          >  1:1165293/1‑150 (MQ=255)
          tgctgTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGc                                                                                                                                      <  1:1078705/150‑1 (MQ=255)
          tgctgTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGc                                                                                                                                      <  2:641887/150‑1 (MQ=255)
          tgctgTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGAATAAAAAATCGTAAAACCGGTGGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGc                                                                                                                                      <  2:1160221/150‑1 (MQ=255)
           gctgTGAAGCATGAATATGCCACAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCa                                                                                                                                     >  2:333109/1‑150 (MQ=255)
           gctgTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCa                                                                                                                                     >  2:328124/1‑150 (MQ=255)
              gtgAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGCGCAGTa                                                                                                                                  >  1:669557/1‑150 (MQ=255)
               tgAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAACTATTCGTTCGGTGTTGCCTTGCAAGTGGTCCGTGTATTGCATATTGTCACGTTAGTGGTTTCAAAATGAGCAGTaa                                                                                                                                 >  1:314669/1‑150 (MQ=255)
                   tcATGAATATGCCAAAGTGGATGTTGATAAATGGCATAAAAAATAGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGGGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAAt                                                                                                                             <  1:740562/149‑1 (MQ=255)
                   gCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGTTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTTCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAAt                                                                                                                             >  1:974698/1‑150 (MQ=255)
                           atGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATTTCCGGaa                                                                                                                     >  1:564951/1‑150 (MQ=255)
                            tGCCAAAGGGGATGATGAAAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGGTTCAAAATGAGCAGTAAAAATGTCCGGAAg                                                                                                                    <  2:122100/150‑1 (MQ=255)
                              ccAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGac                                                                                                                  >  1:591017/1‑150 (MQ=255)
                                 aaGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTTAAAATGTCCGGAAGACAcc                                                                                                               <  1:121650/150‑1 (MQ=255)
                                  aGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCa                                                                                                              <  2:987913/150‑1 (MQ=255)
                                     ggATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTACCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTACGGAAGACACCaaaa                                                                                                           >  2:1142578/1‑150 (MQ=255)
                                       atgatgatAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGTATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAg                                                                                                         <  1:807496/150‑1 (MQ=255)
                                       atgatgatAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAg                                                                                                         <  1:544151/150‑1 (MQ=255)
                                           tgatAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGt                                                                                                     >  2:896923/1‑150 (MQ=255)
                                              tAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGc                                                                                                  <  1:333109/150‑1 (MQ=255)
                                              tAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGc                                                                                                  <  1:328124/150‑1 (MQ=255)
                                                aaTGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCGGGGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAg                                                                                                <  2:495924/150‑1 (MQ=255)
                                                         aaaaaTCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATg                                                                                       <  1:312244/150‑1 (MQ=255)
                                                            aaTCGTAAAACCGGTCGCTTAGTTAAAGCTTTTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTTGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAg                                                                                    <  1:115688/150‑1 (MQ=255)
                                                                gTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGc                                                                                >  1:1182105/1‑150 (MQ=255)
                                                                 tAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCg                                                                               >  1:596080/1‑150 (MQ=255)
                                                                  aaaaCCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGAATGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTGCGGAAGCCACCAAAAAGTTGTCGAAGGGAAGTATCCAGTGGCgg                                                                              >  2:672866/1‑150 (MQ=255)
                                                                      ccGGTCGCATAGTTAAATCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTTGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAgt                                                                          <  1:809318/150‑1 (MQ=255)
                                                                               tAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTTGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGa                                                                 >  1:1184679/1‑150 (MQ=255)
                                                                                       gCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGt                                                         >  1:1154669/1‑150 (MQ=255)
                                                                                        cTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGACCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTg                                                        <  1:1099810/150‑1 (MQ=255)
                                                                                         tATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGa                                                       >  1:875980/1‑150 (MQ=255)
                                                                                           ttCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGAtt                                                     <  2:539190/150‑1 (MQ=255)
                                                                                              gTGCGGGGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTTTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCt                                                  <  2:1346458/150‑1 (MQ=255)
                                                                                                     gtTGCCTTGCAAGTGGTCAGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGaa                                           <  1:896923/150‑1 (MQ=255)
                                                                                                         ccTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGcc                                       <  2:1165293/150‑1 (MQ=255)
                                                                                                          cTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGCCACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCg                                      >  2:976781/1‑150 (MQ=255)
                                                                                                             gCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGAtt                                   >  2:1056756/1‑150 (MQ=255)
                                                                                                             gCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGAtt                                   >  2:464218/1‑150 (MQ=255)
                                                                                                              cAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTg                                  >  2:1128232/1‑150 (MQ=255)
                                                                                                              cAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGAAGTGGGGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATTAATGCCGATTg                                  >  2:828684/1‑150 (MQ=255)
                                                                                                              cAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGAGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTg                                  <  2:918741/150‑1 (MQ=255)
                                                                                                                  tGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTaaa                              >  1:658747/1‑150 (MQ=255)
                                                                                                                   ggTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTaaaa                             >  1:1299705/1‑150 (MQ=255)
                                                                                                                    gTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAACTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCTGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAAc                            <  1:667578/150‑1 (MQ=255)
                                                                                                                      ccGTGGATTGCATATTGTCCCGTTTGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACgg                          >  2:692049/1‑150 (MQ=255)
                                                                                                                       cGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATACCGATTGTAAAACGGc                         >  2:718525/1‑150 (MQ=255)
                                                                                                                        gTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCa                        >  1:380145/1‑150 (MQ=255)
                                                                                                                            aTTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTg                    >  2:48808/1‑150 (MQ=255)
                                                                                                                            aTTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTg                    >  1:269169/1‑150 (MQ=255)
                                                                                                                               gaatatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTg                 <  1:979319/148‑1 (MQ=255)
                                                                                                                               gCATATTGGCCCGATCGTGGTCTCAAAATGAGCAGTAAAAATCTCCGGAAGACACCAAAAAGTTCTCGCAGGGCAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGAGTTCTTGATGAATGCCGTTTTACAATCGGCAttcatca                 <  2:1300723/150‑6 (MQ=255)
                                                                                                                                  tatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCa              >  1:290188/1‑150 (MQ=255)
                                                                                                                                   atTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCAt             >  2:773655/1‑150 (MQ=255)
                                                                                                                                   atTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGACGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGAGTGTAAAACGGCATATGGTGCCAt             >  2:706216/1‑150 (MQ=255)
                                                                                                                                    ttGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCAtt            <  1:27789/150‑1 (MQ=255)
                                                                                                                                     tGTCCCGTTCGAGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGAAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTg           <  2:1214239/150‑1 (MQ=255)
                                                                                                                                     tGTCCCGTTAGTGGTTTCAAAATGCGCAGAAAAAATGTCCGGAAGACTCCAAAAAGTTGTCTCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTTATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTg           <  2:361080/150‑1 (MQ=255)
                                                                                                                                     tGTCCCGTTAGAGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTg           <  2:1219895/150‑1 (MQ=255)
                                                                                                                                      gTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTga          >  2:1246001/1‑150 (MQ=255)
                                                                                                                                            ttAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGCATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATc    <  2:669557/150‑1 (MQ=255)
                                                                                                                                              aGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAc  >  2:569957/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGCCTGGGCTGCTGTGAAGCATGAATATGCCAAAGGGGATGATGATAAATGGCATAAAAAATCGTAAAACCGGTCGCTTAGTTAAAGCTATTCGTGCGGTGTTGCCTTGCAAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAC  >  NC_000913/1272282‑1272573

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