Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 10,619 G→T 6.5% intergenic (‑125/+24) satP ← / ← yaaW succinate‑acetate transporter/UPF0174 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,6190GT 6.5% 302.6 / 7.9 101intergenic (‑125/+24)satP/yaaWsuccinate‑acetate transporter/UPF0174 family protein
Reads supporting (aligned to +/- strand):  ref base G (48/46);  new base T (2/5);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 106 total aligned reads displayed.

ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCC  >  NC_000913/10471‑10768
                                                                                                                                                     |                                                                                                                                                     
aTTAGCCAACTTAGTGTTGCCCATAATGCCGCAAAACTCATCATCGCATGAATGGTGAAATACTTTCCCTGAATAACGGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGCGGGgaaaatga                                                                                                                                                      <  1:1516894/150‑1 (MQ=255)
     ccAACTTATTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGGGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGACAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAAtt                                                                                                                                                 <  2:1292269/150‑1 (MQ=255)
         cTTAGTGTTGCCCATAATTCCTCAAAAATCCTCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTcc                                                                                                                                             <  1:1430295/150‑1 (MQ=255)
         cTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTcc                                                                                                                                             <  1:967183/150‑1 (MQ=255)
         attAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTcc                                                                                                                                             <  1:1423730/149‑1 (MQ=255)
           tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccc                                                                                                                                           <  2:681246/150‑1 (MQ=255)
           tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccc                                                                                                                                           <  2:1396346/150‑1 (MQ=255)
              tgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATATTTTCCATGATTAAGGTTAGGGTTTACGGGCGCGGCATTAATATACAGCCATTGGGCAATTTTGTACAAAATTTTGCGGGTAAATGAAAAATTTACCCCGtt                                                                                                                                       >  1:100533/1‑150 (MQ=255)
              tgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccctca                                                                                                                                        <  2:1012511/150‑4 (MQ=255)
                 tGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTc                                                                                                                                     >  1:691542/1‑150 (MQ=255)
                  gCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATTTTCATTTTcccctcaaaag                                                                                                                                    <  2:545374/150‑8 (MQ=255)
                  gCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTcc                                                                                                                                    >  2:1252906/1‑150 (MQ=255)
                  gCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTcc                                                                                                                                    >  2:812065/1‑150 (MQ=255)
                  gCCCATAATTCCTCAAAAATCATAATCGAATGAATGGTGAGATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATGTTTTGAGGGGTAAATGAATATTTTCACCGGTTTcc                                                                                                                                    >  1:646975/1‑150 (MQ=255)
                   cccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCg                                                                                                                                   >  1:1330947/1‑150 (MQ=255)
                    ccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCgg                                                                                                                                  <  1:286172/150‑1 (MQ=255)
                    ccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCgg                                                                                                                                  <  1:346757/150‑1 (MQ=255)
                        aaTTCCTCAAAAATCATCATCGAATGAATGCTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATc                                                                                                                              <  2:571074/150‑1 (MQ=255)
                         aTTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCa                                                                                                                             <  1:799477/150‑1 (MQ=255)
                            ccTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAc                                                                                                                          <  2:1645386/150‑1 (MQ=255)
                             cTCAAAAATGATCATCGAATGAATGGTGAAATAATTTCCCTGCAAAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAcc                                                                                                                         <  2:301228/150‑1 (MQ=255)
                              tCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCt                                                                                                                        <  1:108175/150‑1 (MQ=255)
                                 aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAg                                                                                                                     >  1:1341620/1‑150 (MQ=255)
                                      catcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGgcg                                                                                                                >  1:1236610/1‑150 (MQ=255)
                                       atcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCCGACCTGAGTGgcgc                                                                                                               >  2:1002003/1‑150 (MQ=255)
                                        tcatcGAATGTATGGTGAAATAATTTCCCTTAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCGCCGGTTTCCGGTATCAGACCTGAGTGGCGCt                                                                                                              <  2:1229695/150‑1 (MQ=255)
                                        tcatcGAATGAATGGTGAAATAATATCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTATCCGCGGTTTCCGGTATCAGACCTGAGTGGCGCt                                                                                                              >  1:300781/1‑150 (MQ=255)
                                         catcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTa                                                                                                             >  2:586148/1‑150 (MQ=255)
                                         catcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTa                                                                                                             >  1:1179614/1‑150 (MQ=255)
                                          atcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTaa                                                                                                            >  1:1616718/1‑150 (MQ=255)
                                           tcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAAc                                                                                                           <  1:672083/150‑1 (MQ=255)
                                             taatgaatgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTTTCAGACCTGAGTGGCGCTAAcc                                                                                                          >  1:390155/2‑149 (MQ=255)
                                             gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAAcc                                                                                                          <  2:390155/149‑1 (MQ=255)
                                             gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCa                                                                                                         >  2:1662355/1‑150 (MQ=255)
                                             gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCa                                                                                                         >  2:724780/1‑150 (MQ=255)
                                             caatgaaagGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGTGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCa                                                                                                         <  2:310199/149‑1 (MQ=255)
                                                     ggTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGgcgc                                                                                                 <  2:693535/150‑1 (MQ=255)
                                                       tGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcag                                                                                               <  1:206558/150‑1 (MQ=255)
                                                           ataatTTCCCTGAATAACTGTAGTGTTTTCAGGGCTCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcaggcag                                                                                           <  1:625595/150‑1 (MQ=255)
                                                            taatTTCCCTGAATAACTGTAGTGTTTTCGGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGcaggcagg                                                                                          >  2:1645463/1‑150 (MQ=255)
                                                            taatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGcaggcagg                                                                                          <  1:1718375/150‑1 (MQ=255)
                                                                 tCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAtt                                                                                     >  1:1128519/1‑150 (MQ=255)
                                                                 tCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAtt                                                                                     <  1:962782/150‑1 (MQ=255)
                                                                 tCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTACCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAtt                                                                                     <  2:859137/150‑1 (MQ=255)
                                                                    cTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGc                                                                                  >  2:756170/1‑150 (MQ=255)
                                                                            cTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGc                                                                          <  1:858183/150‑1 (MQ=255)
                                                                              gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTg                                                                        >  2:507102/1‑150 (MQ=255)
                                                                              gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTg                                                                        <  1:597002/150‑1 (MQ=255)
                                                                               tAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTgg                                                                       >  1:977652/1‑150 (MQ=255)
                                                                                   gtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATc                                                                   >  2:1599994/1‑150 (MQ=255)
                                                                                         aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCAcg                                                             <  1:756170/150‑1 (MQ=255)
                                                                                               cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                       >  1:1733914/1‑150 (MQ=255)
                                                                                               cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                       >  2:560246/1‑150 (MQ=255)
                                                                                               cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                       >  2:132820/1‑150 (MQ=255)
                                                                                                ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTgg                                                                       >  1:518540/1‑133 (MQ=255)
                                                                                                ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTgg                                                                       <  2:518540/133‑1 (MQ=255)
                                                                                                ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAgg                                                      <  2:1296350/150‑1 (MQ=255)
                                                                                                  cATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgcg                                                    <  1:1599994/150‑1 (MQ=255)
                                                                                                   ataattATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgcgc                                                   <  1:644636/150‑1 (MQ=255)
                                                                                                   ataataATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgcgc                                                   >  2:496742/1‑150 (MQ=255)
                                                                                                         aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGTCGATGccc                                              >  2:592826/1‑149 (MQ=255)
                                                                                                         aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCg                                              >  2:423988/1‑150 (MQ=255)
                                                                                                            cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTg                                           >  2:1068086/1‑150 (MQ=255)
                                                                                                            cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCAGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTg                                           >  1:255911/1‑150 (MQ=255)
                                                                                                             aGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCTGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGa                                          >  2:1249178/1‑150 (MQ=255)
                                                                                                              gCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGAc                                         <  1:1252906/150‑1 (MQ=255)
                                                                                                                 aGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGAATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGc                                      <  2:1330947/150‑1 (MQ=255)
                                                                                                                   tGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCtt                                    <  1:1068086/150‑1 (MQ=255)
                                                                                                                    ggggCAGTGTCTACGATCTTTTTACGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCGGCCGGGCGCTAACCATGCGGCGCATGCAGGCGTTTTGCAGTACGGCTGGAATCGTCACGCGATATGCGCTGCCGCTGACCGCttt                                   <  2:1202912/150‑1 (MQ=255)
                                                                                                                     gggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTa                                  >  1:1113277/1‑150 (MQ=255)
                                                                                                                      ggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGAc                                         >  1:418839/1‑142 (MQ=255)
                                                                                                                      ggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGAc                                         <  2:418839/142‑1 (MQ=255)
                                                                                                                      ggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTaa                                 >  1:1647215/1‑150 (MQ=255)
                                                                                                                      ggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTaa                                 >  2:898745/1‑150 (MQ=255)
                                                                                                                       gCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAc                                <  1:808565/150‑1 (MQ=255)
                                                                                                                        cAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcc                               <  2:379011/150‑1 (MQ=255)
                                                                                                                        cAGTGGCAACGATCTTTTGAGGGGAAAATTAAAATTTTCACCGGTTTCCGGTATCAGACCTGAGGGGCGCTAACAATCAGGAGAATCCAGGCGAATTCCAGTACGGCTGGAATCGGAACGCGATAGGCGCTTCCGCTGACCGCTTTAAcc                               <  2:1234546/150‑1 (MQ=255)
                                                                                                                          gtgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATGAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGATTCGTCACGCGATAGGCGCTGCAGCTGGACGCTTTACcccc                             >  2:169774/1‑150 (MQ=255)
                                                                                                                          gtgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcccc                             <  2:691542/150‑1 (MQ=255)
                                                                                                                          gtgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCACGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGAAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcccc                             >  1:1366419/1‑150 (MQ=255)
                                                                                                                            gtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAt                           <  2:1616718/150‑1 (MQ=255)
                                                                                                                                   aTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGTCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGc                    <  2:300781/150‑1 (MQ=255)
                                                                                                                                    tCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAAACATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGcc                   <  2:1220595/150‑1 (MQ=255)
                                                                                                                                      ttttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                 <  1:724780/150‑1 (MQ=255)
                                                                                                                                      ttttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                 <  1:1662355/150‑1 (MQ=255)
                                                                                                                                        ttGAGGGGAAAATGAAAATTTTCGCCGGTTTCCGGTATGAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGTCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAc               <  1:300040/150‑1 (MQ=255)
                                                                                                                                        aaccggggAAAATTTTCATTTTCCCCTGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAc               >  1:982864/5‑150 (MQ=255)
                                                                                                                                         tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAcc              <  2:1128519/150‑1 (MQ=255)
                                                                                                                                          gAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTGTACCCCCATTTAGTGCCGCACCt             >  2:245363/1‑150 (MQ=255)
                                                                                                                                           aGGGGAAAATGAAAATTTTCCCCGGTGTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCGCCTa            >  2:1480241/1‑150 (MQ=255)
                                                                                                                                            ggggAAAATTTTCATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc           >  1:1316000/1‑150 (MQ=255)
                                                                                                                                            ggggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc           >  2:960273/1‑150 (MQ=255)
                                                                                                                                            ggggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc           >  2:1305842/1‑150 (MQ=255)
                                                                                                                                             gggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTACCCCCATTTAGTGCCGCACCTAcc          >  2:858443/1‑149 (MQ=255)
                                                                                                                                             gggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACa          >  1:759635/1‑150 (MQ=255)
                                                                                                                                              ggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAg         <  1:586148/150‑1 (MQ=255)
                                                                                                                                              ggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGAGCAGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCTCCTACAg         >  1:1124160/1‑150 (MQ=255)
                                                                                                                                                aaaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAggg       >  1:1261911/1‑150 (MQ=255)
                                                                                                                                                 aaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGc      >  2:1208354/1‑150 (MQ=255)
                                                                                                                                                     gaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTcc  >  1:984698/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCC  >  NC_000913/10471‑10768

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: