Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,895,244 A→C 5.1% intergenic (+26/+28) 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,2440AC 5.1% 572.9 / 18.0 178intergenic (+26/+28)yieF/adePchromate reductase/adenine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (93/76);  new base C (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 191 total aligned reads displayed.

GCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCA  >  NC_000913/3895097‑3895393
                                                                                                                                                   |                                                                                                                                                     
   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)
           cAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTa                                                                                                                                          >  2:681233/1‑150 (MQ=255)
           cAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTAATCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTa                                                                                                                                          <  2:164566/150‑1 (MQ=255)
             caaCAAAGTTGATCCGCAAACCGGTGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCATGCATGATTGCTGGCGGGttttttat                                                                                                                                        <  1:1213439/149‑4 (MQ=38)
              aaCAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcg                                                                                                                                       <  1:1030140/150‑1 (MQ=255)
               aCAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgc                                                                                                                                      >  1:71270/1‑150 (MQ=255)
                cAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgcg                                                                                                                                     <  2:716015/150‑1 (MQ=255)
                cAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAAcgcg                                                                                                                                     <  1:424947/150‑1 (MQ=255)
                 aaaGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGATCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACgcgc                                                                                                                                    <  2:842960/150‑1 (MQ=255)
                 aaaGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACgcgc                                                                                                                                    <  2:251895/150‑1 (MQ=255)
                   aGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTa                                                                                                                                  >  1:1215521/1‑150 (MQ=255)
                     ttGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACAACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATc                                                                                                                                >  2:1163299/1‑150 (MQ=255)
                       gATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGa                                                                                                                              >  2:1367028/1‑150 (MQ=255)
                        aTCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAt                                                                                                                             <  1:1051422/150‑1 (MQ=255)
                         tCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAtt                                                                                                                            <  1:723067/150‑1 (MQ=255)
                           cGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAtttt                                                                                                                          >  1:313305/1‑150 (MQ=255)
                           cGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGAtttt                                                                                                                          <  1:76778/150‑1 (MQ=255)
                             cAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAg                                                                                                                        <  2:1139355/150‑1 (MQ=255)
                                aCCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGa                                                                                                                     >  2:81847/1‑150 (MQ=255)
                                   ggAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                                  >  2:1006554/1‑150 (MQ=255)
                                   ggAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                                  >  2:178532/1‑150 (MQ=255)
                                   ggAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCg                                                                                                                  <  1:150907/150‑1 (MQ=255)
                                       aaGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGa                                                                                                              <  2:155136/150‑1 (MQ=255)
                                         gTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATa                                                                                                            <  1:1329959/150‑1 (MQ=255)
                                          tgattgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATaa                                                                                                           <  2:93197/150‑1 (MQ=255)
                                           gattgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATaaa                                                                                                          <  2:59556/150‑1 (MQ=255)
                                            attgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAg                                                                                                         <  2:543628/150‑1 (MQ=255)
                                             ttgatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGa                                                                                                        <  2:593010/150‑1 (MQ=255)
                                               gatCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACa                                                                                                      >  2:459374/1‑150 (MQ=255)
                                                      gTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCa                                                                                               <  1:273351/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)
                                                                   ccTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAcagcag                                                                                  <  2:182718/150‑1 (MQ=255)
                                                                     tGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAgcg                                                                                >  2:1416938/1‑150 (MQ=255)
                                                                         agggCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAAc                                                                            >  1:403084/2‑150 (MQ=255)
                                                                           ggCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACga                                                                          <  2:773072/150‑1 (MQ=255)
                                                                            gCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACgat                                                                         >  1:1137856/1‑150 (MQ=255)
                                                                               aTTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACgatgat                                                                      >  2:76671/1‑150 (MQ=255)
                                                                                   aCCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACg                                                                  <  1:81847/150‑1 (MQ=255)
                                                                                    ccGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGc                                                                 <  2:47832/150‑1 (MQ=255)
                                                                                      gCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCAc                                                               >  2:1035912/1‑150 (MQ=255)
                                                                                         tttGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACggg                                                            <  1:1006554/150‑1 (MQ=255)
                                                                                            ggTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCtt                                                         <  2:707395/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)
                                                                                              tGAGTTTATTCAGCTAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACTGGCTTAg                                                       <  1:1163299/150‑1 (MQ=255)
                                                                                              tGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAg                                                       >  1:186307/1‑150 (MQ=255)
                                                                                               gAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGa                                                      >  1:283217/1‑150 (MQ=255)
                                                                                               gAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGa                                                      >  1:1272372/1‑150 (MQ=255)
                                                                                                     aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGc                                                <  1:525313/150‑1 (MQ=255)
                                                                                                       tCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTAAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACTATGATTACGCACGGGCTTAGATCACGCAg                                              <  2:563087/150‑1 (MQ=255)
                                                                                                        cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGa                                             <  2:219469/150‑1 (MQ=255)
                                                                                                        cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGa                                             >  1:1252575/1‑150 (MQ=255)
                                                                                                        cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGa                                             >  1:449268/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)
                                                                                                                       taaataaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtt                              >  1:1013711/1‑150 (MQ=255)
                                                                                                                       taaataaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtt                              >  1:1018883/1‑150 (MQ=255)
                                                                                                                         cgttaaaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtttt                            >  1:1247496/5‑150 (MQ=38)
                                                                                                                         aataaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGAtttt                            >  1:1316382/1‑150 (MQ=255)
                                                                                                                           taaaAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtca                          <  2:11088/150‑1 (MQ=255)
                                                                                                                            aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                         <  1:97662/150‑1 (MQ=255)
                                                                                                                            aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                         >  1:6359/1‑150 (MQ=255)
                                                                                                                            aaaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcat                         >  2:237938/1‑150 (MQ=255)
                                                                                                                             aaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatc                        <  2:1252575/150‑1 (MQ=255)
                                                                                                                             aaaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatc                        <  1:998098/150‑1 (MQ=255)
                                                                                                                              aaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca                       >  2:246886/1‑150 (MQ=255)
                                                                                                                              aaaaCCCGCCATGCATGATTGCTGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGAACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca                       >  2:974796/1‑150 (MQ=255)
                                                                                                                              aaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca                       <  1:82031/150‑1 (MQ=255)
                                                                                                                              aaaaCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca                       <  1:170185/150‑1 (MQ=255)
                                                                                                                                 aCCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGc                    <  2:504580/150‑1 (MQ=255)
                                                                                                                                  cccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTtcatca                       >  2:1015057/1‑146 (MQ=255)
                                                                                                                                  cccGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCa                   <  1:1416938/150‑1 (MQ=255)
                                                                                                                                    cGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGt                 <  2:538622/150‑1 (MQ=255)
                                                                                                                                       cAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAgg              <  2:313305/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           <  2:283217/150‑1 (MQ=255)
                                                                                                                                          cAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGgaga           >  1:1110180/1‑150 (MQ=255)
                                                                                                                                             tCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATaa        <  2:1343667/150‑1 (MQ=255)
                                                                                                                                             tCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATaa        <  1:847190/150‑1 (MQ=255)
                                                                                                                                              catgcatgGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGTGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATaaa       <  2:403084/150‑1 (MQ=255)
                                                                                                                                                tgcatgGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGc     <  1:231678/150‑1 (MQ=255)
                                                                                                                                                tgcatgGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGc     <  1:743680/150‑1 (MQ=255)
                                                                                                                                                  catgGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGccc   >  1:686899/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCA  >  NC_000913/3895097‑3895393

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