Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,895,235 C→G 5.1% intergenic (+17/+37) yieF → / ← adeP chromate reductase/adenine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,895,2350CG 5.1% 617.8 / 16.6 178intergenic (+17/+37)yieF/adePchromate reductase/adenine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (93/76);  new base G (5/4);  total (98/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 185 total aligned reads displayed.

GGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGA  >  NC_000913/3895086‑3895384
                                                                                                                                                     |                                                                                                                                                     
ggAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGc                                                                                                                                                       <  2:227624/150‑1 (MQ=255)
ggAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGc                                                                                                                                                       >  2:171681/1‑150 (MQ=255)
ggAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGc                                                                                                                                                       >  2:101875/1‑150 (MQ=255)
              gcgTGATTCATAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCg                                                                                                                                         <  2:495041/150‑1 (MQ=255)
                   aTTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGCCCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGttt                                                                                                                                    >  2:1180163/1‑150 (MQ=255)
                    ttCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCATGCATGATTGCTGGCGGGtttt                                                                                                                                   <  2:738205/150‑1 (MQ=255)
                    ttCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCATGCATGATTGCTGGCGGGtttt                                                                                                                                   <  1:465833/150‑1 (MQ=255)
                    ttCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGtttt                                                                                                                                   <  1:55318/150‑1 (MQ=255)
                      cAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGc                                                                                                                                                       <  1:1053624/128‑1 (MQ=255)
                      cAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGc                                                                                                                                                       >  2:1053624/1‑128 (MQ=255)
                      cAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTa                                                                                                                                 >  2:681233/1‑150 (MQ=255)
                        gAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAc                                                                                                                               >  1:716015/1‑150 (MQ=255)
                         aaCAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcg                                                                                                                              <  1:1030140/150‑1 (MQ=255)
                         aaCAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcg                                                                                                                              >  1:346480/1‑150 (MQ=255)
                          aCAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgc                                                                                                                             >  1:162148/1‑150 (MQ=255)
                           cAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgcg                                                                                                                            <  1:424947/150‑1 (MQ=255)
                           cAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgcg                                                                                                                            <  2:716015/150‑1 (MQ=255)
                            aaaGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGATCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACgcgc                                                                                                                           <  2:842960/150‑1 (MQ=255)
                            aaaGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACgcgc                                                                                                                           <  1:77524/150‑1 (MQ=255)
                            aaaGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACgcgc                                                                                                                           <  2:251895/150‑1 (MQ=255)
                              aGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTa                                                                                                                         >  2:327660/1‑150 (MQ=255)
                                ttGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACAACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATc                                                                                                                       >  2:1163299/1‑150 (MQ=255)
                                    tCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAtt                                                                                                                   >  1:68811/1‑150 (MQ=255)
                                      cGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAtttt                                                                                                                 >  1:258787/1‑150 (MQ=255)
                                           aCCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGa                                                                                                            >  2:81847/1‑150 (MQ=255)
                                            ccGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAg                                                                                                           >  2:1366991/1‑150 (MQ=255)
                                              ggAGAAGTGTTTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCAAGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                         <  1:1318576/150‑1 (MQ=255)
                                              ggAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                         <  1:150907/150‑1 (MQ=255)
                                              ggAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                         >  2:1006554/1‑150 (MQ=255)
                                               gagaAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGt                                                                                                        >  1:282626/1‑150 (MQ=255)
                                                  aaGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGa                                                                                                     <  2:155136/150‑1 (MQ=255)
                                                     tgattgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATaa                                                                                                  <  1:171681/150‑1 (MQ=255)
                                                     tgattgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATaa                                                                                                  <  2:93197/150‑1 (MQ=255)
                                                      gattgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATaaa                                                                                                 <  2:59556/150‑1 (MQ=255)
                                                        ttgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGa                                                                                               <  2:593010/150‑1 (MQ=255)
                                                          gatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACa                                                                                             >  2:459374/1‑150 (MQ=255)
                                                                 gTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCa                                                                                      <  2:552272/150‑1 (MQ=255)
                                                                 gTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCa                                                                                      >  1:783141/1‑150 (MQ=255)
                                                                  tACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAg                                                                                     <  1:968759/150‑1 (MQ=255)
                                                                   aCGCTGGACCACCTGACCGGTCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAgg                                                                                    <  2:639583/150‑1 (MQ=255)
                                                                   aCGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAgg                                                                                    <  1:1126495/150‑1 (MQ=255)
                                                                    cGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGa                                                                                   >  1:538622/1‑150 (MQ=255)
                                                                    cGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCACCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCCGGTATTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGa                                                                                   >  2:736839/1‑150 (MQ=255)
                                                                         gACCACCTTACCTGGCAATTGAACGCATTTGGTGAGTTTATTCAGCGATTTAAGATCTAAATAAAAAACCCGCCATCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAc                                                                              <  1:604714/150‑1 (MQ=255)
                                                                         gACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAc                                                                              >  1:904701/1‑150 (MQ=255)
                                                                           ccaccTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAcag                                                                            >  1:1214605/1‑150 (MQ=255)
                                                                             accTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAcagca                                                                          >  1:120725/1‑150 (MQ=255)
                                                                                tGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAgcg                                                                       >  2:1416938/1‑150 (MQ=255)
                                                                                    agggCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAAc                                                                   >  1:403084/2‑150 (MQ=255)
                                                                                      ggCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACga                                                                 <  2:773072/150‑1 (MQ=255)
                                                                                      ggCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACga                                                                 >  2:17992/1‑150 (MQ=255)
                                                                                         aaTTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACgatga                                                              >  1:1127740/1‑150 (MQ=255)
                                                                                          aTTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACgatgat                                                             >  2:76671/1‑150 (MQ=255)
                                                                                            tGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTa                                                           <  2:54/150‑1 (MQ=255)
                                                                                              aCCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACg                                                         <  1:81847/150‑1 (MQ=255)
                                                                                              aCCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACg                                                         >  1:1316329/1‑150 (MQ=255)
                                                                                                  cATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACg                                                     <  1:7370/150‑1 (MQ=255)
                                                                                                    tttGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACggg                                                   <  1:1006554/150‑1 (MQ=255)
                                                                                                       ggTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCtt                                                >  2:1274540/1‑150 (MQ=255)
                                                                                                       ggTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCtt                                                >  1:1343667/1‑150 (MQ=255)
                                                                                                        gTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTa                                               >  2:1384606/1‑150 (MQ=255)
                                                                                                          gAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGa                                             >  1:1272372/1‑150 (MQ=255)
                                                                                                           aGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGAt                                            >  1:113915/1‑150 (MQ=255)
                                                                                                           aGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGAt                                            >  2:1123745/1‑150 (MQ=255)
                                                                                                              ttATTCAGAGAGTTAAGATCTAAATAAAAAACACGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCAc                                         >  1:942998/1‑150 (MQ=255)
                                                                                                               tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACg                                        >  1:1233360/1‑150 (MQ=255)
                                                                                                               tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACg                                        >  2:182737/1‑150 (MQ=255)
                                                                                                                aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGc                                       >  2:170185/1‑150 (MQ=255)
                                                                                                                aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGc                                       <  1:525313/150‑1 (MQ=255)
                                                                                                                aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGc                                       <  2:258787/150‑1 (MQ=255)
                                                                                                                   cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGa                                    >  1:449268/1‑150 (MQ=255)
                                                                                                                   cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGa                                    >  1:1252575/1‑150 (MQ=255)
                                                                                                                        aGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTcc                               >  1:1374032/1‑150 (MQ=255)
                                                                                                                         gTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCg                              >  1:62626/1‑150 (MQ=255)
                                                                                                                            aaGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTa                           >  1:1403633/1‑150 (MQ=255)
                                                                                                                              gATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTAcc                         >  1:501859/1‑150 (MQ=255)
                                                                                                                              gATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTAcc                         >  2:993944/1‑150 (MQ=255)
                                                                                                                                 cTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAt                      >  1:289152/1‑150 (MQ=255)
                                                                                                                                  taaataaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtt                     >  1:1013711/1‑150 (MQ=255)
                                                                                                                                  taaataaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtt                     >  1:1018883/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                >  1:6359/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                >  2:237938/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                <  1:97662/150‑1 (MQ=255)
                                                                                                                                        aaaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatc               >  1:748999/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatc               <  1:998098/150‑1 (MQ=255)
                                                                                                                                        aaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatc               <  2:1252575/150‑1 (MQ=255)
                                                                                                                                         aaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca              >  2:246886/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGAACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca              >  2:974796/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca              >  2:1080401/1‑150 (MQ=255)
                                                                                                                                          aaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCAc             <  2:68811/150‑1 (MQ=255)
                                                                                                                                            aCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGc           <  2:504580/150‑1 (MQ=255)
                                                                                                                                             cccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCa          >  1:821974/1‑150 (MQ=255)
                                                                                                                                             cccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCa          >  2:298343/1‑150 (MQ=255)
                                                                                                                                             cccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCa          <  1:1416938/150‑1 (MQ=255)
                                                                                                                                              ccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAg         >  2:787757/1‑150 (MQ=255)
                                                                                                                                               cGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGt        <  2:538622/150‑1 (MQ=255)
                                                                                                                                                   aGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGATCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGga    <  2:942998/150‑1 (MQ=255)
                                                                                                                                                   aGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGga    <  2:1316382/150‑1 (MQ=255)
                                                                                                                                                    gCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGgag   >  2:1192716/1‑150 (MQ=255)
                                                                                                                                                     cAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGgaga  >  1:1110180/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGA  >  NC_000913/3895086‑3895384

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