Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,132 T→G 5.4% A75A (GCA→GCC) ‡ cutA ← divalent‑cation tolerance protein, copper sensitivity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,1320TG 5.4% 415.3 / 8.8 138A75A (GCA→GCC) ‡cutAdivalent‑cation tolerance protein, copper sensitivity
Reads supporting (aligned to +/- strand):  ref base T (76/54);  new base G (4/4);  total (80/58)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 138 total aligned reads displayed.

CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAA  >  NC_000913/4364984‑4365279
                                                                                                                                                    |                                                                                                                                                   
ccGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTg                                                                                                                                                    <  2:1648770/150‑1 (MQ=255)
ccGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTAGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGAAGTg                                                                                                                                                    >  1:1025153/1‑150 (MQ=255)
ccGGCAAAAACGGAAGTGCAGAAAAGTAGCAGGATCAGCGTAAAGATGCGGTGAGCCATGAGAGGTAACCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTg                                                                                                                                                    <  2:1834169/150‑1 (MQ=255)
   gCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCct                                                                                                                                                 >  1:1736919/1‑150 (MQ=255)
    cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCATCGTAAAGATGCGTTGAGCCATGAGACGTAATCTGTGTCTCCGTGTGGCACAGGTAAACCCAGAAGTTCCGGGGTGTGATATGGATGATGCGACTTCATGCATTCCAGCAGTGCctg                                                                                                                                                >  1:935877/1‑150 (MQ=255)
            taaGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatg                                                                                                                                        <  2:1876236/149‑5 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatg                                                                                                                                        <  1:3560/150‑5 (MQ=255)
                tgctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATAc                                                                                                                                    <  1:555060/150‑1 (MQ=255)
                 gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                   >  2:1452955/1‑150 (MQ=255)
                 gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                   >  2:1016558/1‑150 (MQ=255)
                 gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                   <  1:1570834/150‑1 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  1:1040484/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:844127/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:836856/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:409472/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:1728313/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTATGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:1029013/1‑150 (MQ=255)
                   tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACAAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                 >  2:1022408/1‑150 (MQ=255)
                    gcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTa                                                                                                                                >  1:775128/1‑150 (MQ=255)
                    gcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTa                                                                                                                                >  1:1298110/1‑150 (MQ=255)
                    gcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAAACAGAAGTTCCGGTGTTTGATATGCATGATGAGCCGTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTa                                                                                                                                >  2:1175800/1‑150 (MQ=255)
                     cAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAg                                                                                                                               >  2:1563585/1‑150 (MQ=255)
                     cAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGAAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAg                                                                                                                               <  1:662860/150‑1 (MQ=255)
                       aaGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtt                                                                                                                             >  2:862263/1‑150 (MQ=255)
                       aaGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGAGGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtt                                                                                                                             <  1:1498950/150‑1 (MQ=255)
                         gTAGCAGGCTCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtttt                                                                                                                           <  1:337696/150‑1 (MQ=255)
                         gTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtttt                                                                                                                           >  1:1503120/1‑150 (MQ=255)
                           aGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTa                                                                                                                         >  2:193907/1‑150 (MQ=255)
                             cAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaaa                                                                                                                       <  1:1405930/150‑1 (MQ=255)
                              aGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTATCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaaaa                                                                                                                      >  2:126341/1‑150 (MQ=255)
                               ggATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAat                                                                                                                     <  2:797750/150‑1 (MQ=255)
                                gATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatc                                                                                                                    >  1:1501745/1‑150 (MQ=255)
                                   cAGCGTAAATATGCGTTGAGCCATGAGAGGTAATCTGTGGCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatcatc                                                                                                                 <  1:1334267/150‑1 (MQ=255)
                                     gagTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg                                                                                                               <  2:1148556/148‑1 (MQ=255)
                                      cttAAAGATGCGTTGAGGCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTTATACCAGATGTTACGGGGTTTGATATGTATGATGAGACTTCAGGCATTCCAGCACTGGCTGCTGGTGAGATACGGTACTTTTTAAAATCATCTGc                                                                                                              <  1:1865799/148‑1 (MQ=255)
                                      cGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGc                                                                                                              <  2:1791880/150‑1 (MQ=255)
                                         aaaGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACt                                                                                                           <  2:502019/150‑1 (MQ=255)
                                           aGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTc                                                                                                         <  2:313660/150‑1 (MQ=255)
                                              tGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTa                                                                                                      >  2:911111/1‑150 (MQ=255)
                                                cGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTAtt                                                                                                    >  1:721532/1‑150 (MQ=255)
                                                 gTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTc                                                                                                   <  1:618601/150‑1 (MQ=255)
                                                       ccATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTc                                                                                             >  2:49925/1‑150 (MQ=255)
                                                         gtGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCa                                                                                           >  1:408916/2‑150 (MQ=255)
                                                          tGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAg                                                                                          <  1:1261051/150‑1 (MQ=255)
                                                           gagaGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGc                                                                                         <  2:223402/150‑1 (MQ=255)
                                                            agagGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCt                                                                                        <  1:1645628/150‑1 (MQ=255)
                                                                gTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTAcc                                                                                    >  2:405389/1‑150 (MQ=255)
                                                                   aTCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTc                                                                                 >  2:1090463/1‑150 (MQ=255)
                                                                     cTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGATACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTccc                                                                               <  1:1796306/150‑1 (MQ=255)
                                                                      tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                              <  1:501816/150‑1 (MQ=255)
                                                                      tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCATGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACATTACCa                                                                              >  2:1219458/1‑150 (MQ=255)
                                                                       gtgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAg                                                                             <  1:1777740/150‑1 (MQ=255)
                                                                        tgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGt                                                                            <  2:1618417/150‑1 (MQ=255)
                                                                           ctcCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAAt                                                                         >  1:1128316/1‑150 (MQ=255)
                                                                             ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                       >  1:937632/1‑150 (MQ=255)
                                                                               gtgtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagag                                                                     <  1:795502/150‑1 (MQ=255)
                                                                                 gtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagag                                                                   >  1:1808056/1‑150 (MQ=255)
                                                                                  tgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCATCTTACCTTCCCAGTAATagagaga                                                                  <  2:1754791/150‑1 (MQ=255)
                                                                                   gtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagagag                                                                 >  2:1721932/1‑150 (MQ=255)
                                                                                   gtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagagag                                                                 >  2:500765/1‑150 (MQ=255)
                                                                                    tatcAGGTAAAACCAGATGTTCCGGGGTTTGATATTGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATTCGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTGCAGATTACCTTCCCAGTAATAGATAGAgg                                                                <  1:803893/147‑1 (MQ=255)
                                                                                     aaCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCAAGTAATATATAGAGGt                                                               >  2:606646/1‑150 (MQ=255)
                                                                                      aCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTa                                                              >  1:1562878/1‑150 (MQ=255)
                                                                                      aCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTa                                                              >  1:1259041/1‑150 (MQ=255)
                                                                                      aCAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTa                                                              >  1:270211/1‑150 (MQ=255)
                                                                                       cAGGTAAAACCAGGAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAg                                                             >  2:1314940/1‑150 (MQ=255)
                                                                                        aGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgc                                                            >  1:960465/1‑150 (MQ=255)
                                                                                         ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                           >  1:1562107/1‑150 (MQ=255)
                                                                                              aaCCAGGAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                      >  2:819945/1‑150 (MQ=255)
                                                                                              aaCCAGAAGTGCCGGGGTGTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                      <  2:360665/150‑1 (MQ=255)
                                                                                               aCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACATCGTATTCTTGCTCCAGCATACCTTCCCAGTAATAGAGAGAGGTAGCGCCgggg                                                     >  1:1657050/1‑150 (MQ=255)
                                                                                                 cagcagTTCCGTGGTTTGATATGGATTATGAGACTTCAGGCATTCCAGCAGTGCCTGATGGTGAGATACGGTAGTTTTTAACATCATCTGCCCTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGTGGAt                                                   <  2:900309/150‑1 (MQ=255)
                                                                                                   gAAGTTCCTGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGAAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCa                                                 <  2:356357/150‑1 (MQ=255)
                                                                                                   gAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCa                                                 <  1:401949/150‑1 (MQ=255)
                                                                                                      gTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAgg                                              >  1:1860141/1‑150 (MQ=255)
                                                                                                       ttCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGt                                             <  2:1006472/150‑1 (MQ=255)
                                                                                                           ggggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcgcg                                         <  2:721532/150‑1 (MQ=255)
                                                                                                           ggggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcgcg                                         >  2:229934/1‑150 (MQ=255)
                                                                                                              gTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAg                                      >  1:1840337/1‑150 (MQ=255)
                                                                                                                ttGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgc                                    <  1:979125/150‑1 (MQ=255)
                                                                                                                 tGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgcc                                   >  1:478530/1‑150 (MQ=255)
                                                                                                                    aatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgcc                                <  1:388758/149‑1 (MQ=255)
                                                                                                                         atgatgAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGTGATCAAGGTCGCGCAGGCAGCCAGttt                           >  1:228330/1‑150 (MQ=255)
                                                                                                                             tgaGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTcc                       >  2:54817/1‑150 (MQ=255)
                                                                                                                             tgaGACTTAAGGCATTCCAGCAGGGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTcc                       <  2:246914/150‑1 (MQ=255)
                                                                                                                              gagaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCg                      >  2:1335502/1‑150 (MQ=255)
                                                                                                                              gagaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCGCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCg                      <  1:1631011/150‑1 (MQ=255)
                                                                                                                              gagaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTATAAAATCATCAGCACAGCGTATTCTTGCTCCAGCTTCCCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCg                      <  1:325964/150‑1 (MQ=255)
                                                                                                                                gaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATTGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                    >  1:1693998/1‑150 (MQ=255)
                                                                                                                                gaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                    <  1:746923/150‑1 (MQ=255)
                                                                                                                                  cTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAg                  >  1:1637825/1‑150 (MQ=255)
                                                                                                                                  cTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAg                  >  1:1908311/1‑150 (MQ=255)
                                                                                                                                     cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCAc               >  1:22311/1‑150 (MQ=255)
                                                                                                                                          tttCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg          >  1:441929/2‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAgg                                     >  1:1905621/4‑123 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAgg                                     <  2:1905621/120‑1 (MQ=255)
                                                                                                                                            accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg        >  2:1218391/2‑150 (MQ=255)
                                                                                                                                            accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg        >  2:347164/2‑150 (MQ=255)
                                                                                                                                              cagcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcggc      >  1:518568/1‑150 (MQ=255)
                                                                                                                                                  agTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGGTTTTCCGCCAGCACTTTGGCGGCTaaa  <  1:1712505/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAA  >  NC_000913/4364984‑4365279

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