Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,610,547 T→C 6.2% E396G (GAA→GGA)  yhhS ← putative arabinose efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,610,5470TC 6.2% 281.8 / 7.3 97E396G (GAA→GGA) yhhSputative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (40/50);  new base C (2/4);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

TGGGTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGCGATAAATCCATAAATACG  >  NC_000913/3610400‑3610696
                                                                                                                                                   |                                                                                                                                                      
tGGGTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCg                                                                                                                                                     >  2:1652752/1‑150 (MQ=255)
 gggTAACGTGGGTAATCAACTGCTGCCCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCg                                                                                                                                                     >  2:1463662/1‑150 (MQ=255)
 gggTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCg                                                                                                                                                     <  2:707211/150‑1 (MQ=255)
  ggTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCgg                                                                                                                                                    <  1:514098/150‑1 (MQ=255)
  ggTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCAGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATGCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCgg                                                                                                                                                    <  1:289669/150‑1 (MQ=255)
   gTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGCGGCGGCCTCAGGGACGTGTTCCGGa                                                                                                                                                   >  1:154999/1‑150 (MQ=255)
   gTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGa                                                                                                                                                   <  1:1672905/150‑1 (MQ=255)
     aaCGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCCGGAACACGTCCCTGAgg                                                                                                                                                 <  2:165678/150‑1 (MQ=255)
      aCGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGACCGCTGACGATCAGCAAAAAGCGCAGACCAACGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCCGGAACACGTCCCTGAGGc                                                                                                                                                <  2:149398/150‑1 (MQ=255)
       cGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCCGGAACACGTCCCTGAGGcc                                                                                                                                               <  2:1155426/150‑1 (MQ=255)
       cGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGcc                                                                                                                                               <  2:1212377/150‑1 (MQ=255)
          gggTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGATTTAAGATGATGAGGCGGCCTCCGGAACACGTCCCTGAGGCCgcc                                                                                                                                            <  2:1252525/150‑3 (MQ=255)
           ggTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGttt                                                                                                                                          <  1:1473912/150‑1 (MQ=255)
            gTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGtttt                                                                                                                                         <  1:322257/150‑1 (MQ=255)
              aaTCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCATCAAAAATCGCAGACCACCGCAGAAAACATCATTAATACGCTTGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGtttttt                                                                                                                                       <  1:1510793/150‑1 (MQ=255)
              aaTCAACTGCTGACCAACCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGtttttt                                                                                                                                       <  1:150205/150‑1 (MQ=255)
               aTCAACTTCTGTCCATGCAAATTACGCTGCCCGCTGACGATCTGCTAAAAGCGCAGACCACCGCAGAAAACATCAATAATACGCTGGTTCTTCAGTAAGTTTAAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGttttttt                                                                                                                                      <  2:1381750/150‑1 (MQ=255)
               aTCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACAGCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGttttttt                                                                                                                                      <  1:681/150‑1 (MQ=255)
               aTCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCCCAGACCGCCGCAGAAACCATCATTAATACGCTGGATATTGAGTAAGTGTTAAGATGATGATGCGGACTCAGTGACGAGTTCCGTAGTCCGgtttttt                                                                                                                                      >  2:93418/1‑150 (MQ=255)
                 cAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTAAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTaa                                                                                                                                    >  2:1102172/1‑150 (MQ=255)
                   aCTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATc                                                                                                                                  >  1:435897/1‑150 (MQ=255)
                   aCTGCTGACCATGCAAATTACGCTGCCCGCTGACGAACAAAAAAAACCGCTCCCCAACGCAGAAAAAATCATAAATACGGTGGTTATTCTTTAAGTTTTAAGATGATGAGCCGGCATCAGGGACGTGTTCCGAAGGCCGTTTTTTTAATc                                                                                                                                  >  2:1328524/1‑150 (MQ=255)
                    ctgctgACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCg                                                                                                                                 <  1:127945/150‑1 (MQ=255)
                      gctgACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGcc                                                                                                                               <  1:752667/150‑1 (MQ=255)
                       ctgACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCa                                                                                                                              <  1:639200/150‑1 (MQ=255)
                       ctgACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCa                                                                                                                              <  1:991249/150‑1 (MQ=255)
                         gACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACg                                                                                                                            <  2:1631463/150‑1 (MQ=255)
                          aCCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGt                                                                                                                           >  2:1561000/1‑150 (MQ=255)
                           ccATGCAAATTACGCGGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTCGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTc                                                                                                                          <  1:1084689/150‑1 (MQ=255)
                            cATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCa                                                                                                                         >  2:1529623/1‑150 (MQ=255)
                             aTGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAg                                                                                                                        >  1:997797/1‑150 (MQ=255)
                              tGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGt                                                                                                                       >  1:584728/1‑150 (MQ=255)
                                cAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGtaa                                                                                                                     >  1:340072/1‑150 (MQ=255)
                                  aaTTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGtaata                                                                                                                   <  1:1542657/150‑1 (MQ=255)
                                    atACGCTGCCCGCTGACGATCTGCAAACAGCGCAGACCACCTCAGACAACATAATTAATTCGCTGGTTATTCAGTAAGTCTTAACATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAAc                                                                                                                 <  1:292060/149‑1 (MQ=255)
                                     tACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTACTACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAAcg                                                                                                                <  1:439213/150‑1 (MQ=255)
                                      aCGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTTAATCGCCACGTCAGTAATAAcg                                                                                                                >  2:710021/1‑150 (MQ=255)
                                         cTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGAt                                                                                                            <  1:222391/150‑1 (MQ=255)
                                         cTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGAt                                                                                                            >  1:761023/1‑150 (MQ=255)
                                          tGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGAtt                                                                                                           >  1:428617/1‑150 (MQ=255)
                                           gCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTg                                                                                                          >  2:1325422/1‑150 (MQ=255)
                                           gCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTg                                                                                                          <  1:1346059/150‑1 (MQ=255)
                                            cccGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGc                                                                                                         <  2:1731293/150‑1 (MQ=255)
                                             ccGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACCTCAGTAATAACGCGCTTGCg                                                                                                        >  1:1037833/1‑150 (MQ=255)
                                               gCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGAc                                                                                                      >  1:988383/1‑150 (MQ=255)
                                                cTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGAcc                                                                                                     <  2:781167/150‑1 (MQ=255)
                                                  gACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAg                                                                                                   <  2:55812/150‑1 (MQ=255)
                                                    cGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTc                                                                                                 >  2:967220/1‑150 (MQ=255)
                                                     gATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTcc                                                                                                <  1:700607/150‑1 (MQ=255)
                                                     gATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTcc                                                                                                >  2:477584/1‑150 (MQ=255)
                                                      aTCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACAAGTccc                                                                                               >  2:827723/1‑150 (MQ=255)
                                                        cAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCccgc                                                                                             >  1:533884/1‑150 (MQ=255)
                                                         aGCAAATTGCGCAGACCCCCGCAGAATCCATCAGTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGTGGCGGGTTGTTTAATCGCCCCGTCAGTAATAACGCGATTGCGACCATGCccgcc                                                                                            >  2:917117/1‑150 (MQ=255)
                                                         aGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGACACGTCAGTAATAACGCGATTGCGACCAGTCccgcc                                                                                            >  1:121672/1‑150 (MQ=255)
                                                          gCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCccgccg                                                                                           >  1:754940/1‑150 (MQ=255)
                                                          gCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCccgccg                                                                                           <  1:1118410/150‑1 (MQ=255)
                                                           cAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCccgccgc                                                                                          <  2:1583047/150‑1 (MQ=255)
                                                           cAAAAAGCGCAGAACACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCccgccgc                                                                                          <  2:1591605/150‑1 (MQ=255)
                                                                 gcgcAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAg                                                                                    <  2:340072/150‑1 (MQ=255)
                                                                 ccgcAGTCCACCGCACAAAACATCATTTATAGGCTGGTTATTCCGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAg                                                                                    <  2:1505162/149‑1 (MQ=255)
                                                                       accaccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGGGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAAt                                                                              <  1:1652487/150‑1 (MQ=255)
                                                                        ccaccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATc                                                                             >  1:1551145/1‑150 (MQ=255)
                                                                        ccaccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATc                                                                             >  1:1692024/1‑150 (MQ=255)
                                                                           ccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCAcc                                                                          >  2:1382040/1‑150 (MQ=255)
                                                                           ccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCAACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCAcc                                                                          <  1:1118912/150‑1 (MQ=255)
                                                                           ccGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCAGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCAcc                                                                          >  2:1194779/1‑150 (MQ=255)
                                                                            cGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCg                                                                         >  2:927557/1‑150 (MQ=255)
                                                                              cAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGt                                                                       <  1:1529623/150‑1 (MQ=255)
                                                                                gAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCTCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTAc                                                                     >  1:146054/1‑150 (MQ=255)
                                                                                gAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGCAGCCCGTTTTTTTAATAGCCACGTCAGTAATAACGCGATTGCGACCAGTCGCGCCGCCGCCAGATAAATCACCGGTAc                                                                     >  1:966636/1‑150 (MQ=255)
                                                                                   aaCATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTAcgc                                                                   <  1:710021/150‑1 (MQ=255)
                                                                                   aaCATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTAcgcc                                                                  <  2:328532/150‑1 (MQ=255)
                                                                                    aCATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTAcgccc                                                                 <  1:1715563/150‑1 (MQ=255)
                                                                                         atTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCa                                                            <  2:971538/150‑1 (MQ=255)
                                                                                            aaTACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCt                                                         <  2:420695/150‑1 (MQ=255)
                                                                                            aaTACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCt                                                         <  2:614054/150‑1 (MQ=255)
                                                                                                   tGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCAcc                                                  <  2:515655/150‑1 (MQ=255)
                                                                                                    ggTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCTGATAAATCACCGGTACGCCCGCCCAGCTCATCCcca                                                 >  1:386350/1‑150 (MQ=255)
                                                                                                        aTTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCAccagccc                                             >  2:1692665/1‑150 (MQ=255)
                                                                                                            aGTAAGTTTTTAGATGATGAGGCGGCCTCAGGGACGTGATCCGGAGGCCGTGTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCTGCTCATCACcagcccagcc                                         <  1:204652/150‑1 (MQ=255)
                                                                                                            aGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCCCGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCCCCGGTACGCCCGCCCAGCTCATCACccgccccgcc                                         >  2:1621676/1‑150 (MQ=255)
                                                                                                               taGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGt                                      <  1:1652752/149‑1 (MQ=255)
                                                                                                                 gTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTgg                                    <  1:1458062/150‑1 (MQ=255)
                                                                                                                 gTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTGCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTgg                                    <  2:1355442/150‑1 (MQ=255)
                                                                                                                  ttttAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGt                                   >  1:664998/1‑150 (MQ=255)
                                                                                                                   tttAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGGTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTc                                  <  2:1248449/150‑1 (MQ=255)
                                                                                                                    ttAAGATGATGTGGCGGCGTCAGGGACGTGTTCCGGAGGCAGTTTTTTAAATCGCCACGTGAGTAATAACGCGTTTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCAAGTGGTcc                                 <  2:1127967/150‑1 (MQ=255)
                                                                                                                      aaGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAg                               >  1:402729/1‑150 (MQ=255)
                                                                                                                        gatgatGAGGCGGCCTCAGGGACGTGTTCGGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACTACCGGCCATCTCATCACAAGCCCAGCCAGGGGAGCAGTc                             >  2:1051326/1‑150 (MQ=255)
                                                                                                                           gatgaGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACg                          <  1:11585/150‑1 (MQ=255)
                                                                                                                              gaGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCa                       <  1:597606/150‑1 (MQ=255)
                                                                                                                              aaaacggcCTCCGGAACACGTCCCTGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCa                       >  1:1724403/5‑150 (MQ=255)
                                                                                                                               aGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGCCCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCGCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCCCGCCac                      >  1:560142/1‑149 (MQ=255)
                                                                                                                                 acggcCTCCGGAACACGTCCCTGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGc                    >  2:1718908/2‑150 (MQ=255)
                                                                                                                                          gggACGTGCTCCGGAGGCCGTTTTTTTCATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGCGATAAATcc           <  1:333680/150‑1 (MQ=255)
                                                                                                                                              cGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGCGATAAATCCATaa       <  2:716180/150‑1 (MQ=255)
                                                                                                                                                 gtTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCCCCGGTACGCCCGCCCAGCTCATCAGCAGCCCAGCAAGTGGTCCAGTCAAGCCAAGCGATAAATCCATAAATa    >  2:281633/1‑150 (MQ=255)
                                                                                                                                                   tCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGCGATAAATCCATAAATACg  <  2:1445741/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGGGTAACGTGGGTAATCAACTGCTGACCATGCAAATTACGCTGCCCGCTGACGATCAGCAAAAAGCGCAGACCACCGCAGAAAACATCATTAATACGCTGGTTATTCAGTAAGTTTTAAGATGATGAGGCGGCCTCAGGGACGTGTTCCGGAGGCCGTTTTTTTAATCGCCACGTCAGTAATAACGCGATTGCGACCAGTCCCGCCGCCGCCAGATAAATCACCGGTACGCCCGCCCAGCTCATCACCAGCCCAGCCAGTGGTCCAGTCACGCCAAGCGATAAATCCATAAATACG  >  NC_000913/3610400‑3610696

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: