Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,134 C→A 7.7% A75S (GCA→TCA) ‡ cutA ← copper binding protein CutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,1340CA 7.7% 1053.5 / 46.3 316A75S (GCA→TCA) ‡cutAcopper binding protein CutA
Reads supporting (aligned to +/- strand):  ref base C (154/135);  new base A (17/8);  total (172/144)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

CGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAATC  >  NC_000913/4364985‑4365281
                                                                                                                                                     |                                                                                                                                                   
cGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCCGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGc                                                                                                                                                     >  1:1003076/1‑150 (MQ=255)
 ggCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGcc                                                                                                                                                    >  1:759522/1‑150 (MQ=255)
   cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctg                                                                                                                                                  >  2:1631328/1‑150 (MQ=255)
   cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctg                                                                                                                                                  >  1:257185/1‑150 (MQ=255)
   cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctg                                                                                                                                                  >  1:184731/1‑150 (MQ=255)
   cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCATAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctg                                                                                                                                                  >  2:1959434/1‑150 (MQ=255)
    aaaaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgc                                                                                                                                                 <  2:305440/150‑1 (MQ=255)
         cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAATCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                            <  2:1670692/150‑1 (MQ=255)
         cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                            <  2:776418/150‑1 (MQ=255)
           gAAGTGCTGCAAAGTAGCAGGATCAGCGTAACGATGCGGTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgag                                                                                                                                          <  2:286818/150‑1 (MQ=255)
           gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgag                                                                                                                                          <  2:485370/150‑1 (MQ=255)
            aaGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATAATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatgc                                                                                                                                         <  2:756597/150‑6 (MQ=255)
              gTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATa                                                                                                                                       <  2:190225/150‑1 (MQ=255)
                gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                     >  2:1094155/1‑150 (MQ=255)
                gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                     >  2:1292776/1‑150 (MQ=255)
                  tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGt                                                                                                                                   >  1:1877743/1‑150 (MQ=255)
                        gTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtttt                                                                                                                             >  2:1042034/1‑150 (MQ=255)
                         tAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTTATATGTATGATGAGACTTCAGGCATTCCATCAGTGCGTTCTGGTGAGATAAGGTAGttttt                                                                                                                            <  2:831837/150‑1 (MQ=255)
                              ggAACAGCGTAAAGATACGTTGAGCCATGACAGCAAATCTCTGACTCCGTGTGAAACAGGTAAAACAACAAGTTCCGGGGTTTCAAATGCAACATGAGACTTAAGGCATTCCAGCAGTACATGCTGGTGAGATACGGTAGTTTTTAAAat                                                                                                                       >  2:1390586/1‑150 (MQ=255)
                                   aGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCt                                                                                                                  <  1:302200/150‑1 (MQ=255)
                                    gCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg                                                                                                                 >  1:1490260/1‑150 (MQ=255)
                                      gTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCa                                                                                                               >  1:1275103/1‑150 (MQ=255)
                                        aaaGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACt                                                                                                             >  1:561351/1‑150 (MQ=255)
                                         aaGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACtt                                                                                                            >  2:1182342/1‑150 (MQ=255)
                                           gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCg                                                                                                          >  2:1953400/1‑150 (MQ=255)
                                           gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCg                                                                                                          <  1:1182342/150‑1 (MQ=255)
                                           gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGAATTGGAGCAGTGCCTGCTGGTGAGATACTGTAGTTTTTAAAATCATCTGCACTTCg                                                                                                          <  2:1859150/150‑1 (MQ=255)
                                               cGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTAtt                                                                                                      >  2:1112444/1‑150 (MQ=255)
                                                gTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTc                                                                                                     >  2:556451/1‑150 (MQ=255)
                                                gTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTc                                                                                                     >  2:1722079/1‑150 (MQ=255)
                                                 ttGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCt                                                                                                    <  1:791199/150‑1 (MQ=255)
                                                   gAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTg                                                                                                  <  1:1457578/150‑1 (MQ=255)
                                                     gCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCt                                                                                                >  1:1765808/1‑150 (MQ=255)
                                                       cATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTcc                                                                                              >  2:874775/1‑150 (MQ=255)
                                                       cATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTcc                                                                                              <  1:27786/150‑1 (MQ=255)
                                                       cATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTcc                                                                                              >  2:1525760/1‑150 (MQ=255)
                                                        aTGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGTCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCa                                                                                             >  1:330605/1‑150 (MQ=255)
                                                         tGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAg                                                                                            >  1:76970/1‑150 (MQ=255)
                                                               gTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTAcc                                                                                      <  2:1490260/150‑1 (MQ=255)
                                                               gTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTAcc                                                                                      >  1:600931/1‑150 (MQ=255)
                                                                tAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCt                                                                                     >  1:1506765/1‑150 (MQ=255)
                                                                 aaTCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCtt                                                                                    >  1:1057617/1‑150 (MQ=255)
                                                                   tCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTcc                                                                                  <  2:903660/150‑1 (MQ=255)
                                                                    cTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTccc                                                                                 <  1:1631328/150‑1 (MQ=255)
                                                                     tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                                <  1:1292776/150‑1 (MQ=255)
                                                                     tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACATCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                                >  2:421952/1‑150 (MQ=255)
                                                                      gtgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAg                                                                               <  2:257185/150‑1 (MQ=255)
                                                                      gtgtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAg                                                                               >  2:1969968/1‑150 (MQ=255)
                                                                        gtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTa                                                                             <  1:946413/150‑1 (MQ=255)
                                                                           tcCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATa                                                                          >  2:1331022/1‑150 (MQ=255)
                                                                           tcCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATa                                                                          >  1:1966703/1‑150 (MQ=255)
                                                                            ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                         >  2:415848/1‑150 (MQ=255)
                                                                            ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                         >  2:1237498/1‑150 (MQ=255)
                                                                                      cAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAg                                                               >  2:1998906/1‑150 (MQ=255)
                                                                                      cAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAg                                                               >  1:359606/1‑150 (MQ=255)
                                                                                        ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                             >  2:1411229/1‑150 (MQ=255)
                                                                                         gTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcgc                                                            >  1:1160540/1‑150 (MQ=255)
                                                                                          tAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGcc                                                           >  1:1543432/1‑150 (MQ=255)
                                                                                           aaaaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCg                                                          >  1:84771/1‑150 (MQ=255)
                                                                                            aaaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgg                                                         >  1:1010853/1‑150 (MQ=255)
                                                                                              aCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgggg                                                       >  2:105623/1‑150 (MQ=255)
                                                                                              aCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgggg                                                       >  1:614154/1‑150 (MQ=255)
                                                                                                cAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGAt                                                     <  1:385207/150‑1 (MQ=255)
                                                                                                cAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGAt                                                     <  1:1026698/150‑1 (MQ=255)
                                                                                                  gAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCa                                                   <  2:1160540/150‑1 (MQ=255)
                                                                                                  caaGTTCCGTCTTATGATATAGATTATTATACTTCAGGCATTAAAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCa                                                   >  2:528136/2‑150 (MQ=255)
                                                                                                    aGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAg                                                 <  1:417905/150‑1 (MQ=255)
                                                                                                     gTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAgg                                                <  1:1378756/150‑1 (MQ=255)
                                                                                                       tCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTc                                              >  1:1321785/1‑150 (MQ=255)
                                                                                                        ccGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcg                                             >  2:1655422/1‑150 (MQ=255)
                                                                                                            ggTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCa                                         >  2:226693/1‑150 (MQ=255)
                                                                                                               ttGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgc                                      >  2:1543945/1‑150 (MQ=255)
                                                                                                                tGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGATGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgcc                                     <  1:730871/150‑1 (MQ=255)
                                                                                                                   tatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgcc                                  >  2:958423/1‑150 (MQ=255)
                                                                                                                     tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTTGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                                <  1:90217/150‑1 (MQ=255)
                                                                                                                      ggAAGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGt                               <  1:1366447/150‑1 (MQ=255)
                                                                                                                        atgatgAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGttt                             >  1:1567428/1‑150 (MQ=255)
                                                                                                                          gatgaGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGTTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCACAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGttttt                           <  2:1275103/150‑1 (MQ=255)
                                                                                                                            tgaGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTcc                         <  2:1765808/150‑1 (MQ=255)
                                                                                                                             gagaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGAGACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCg                        <  1:635360/150‑1 (MQ=255)
                                                                                                                               gaCTTCAGGCATTCCAGCAGTGCCTGTTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                      <  1:1115041/150‑1 (MQ=255)
                                                                                                                                  ttCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGc                   <  1:1650185/150‑1 (MQ=255)
                                                                                                                                  ttCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGc                   <  2:1003076/150‑1 (MQ=255)
                                                                                                                                    cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCa                     <  1:498210/146‑1 (MQ=255)
                                                                                                                                    cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCAc                 >  1:478884/1‑150 (MQ=255)
                                                                                                                                      ggCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACtt               >  2:1807276/1‑150 (MQ=255)
                                                                                                                                        ctcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTg             >  1:525101/5‑150 (MQ=255)
                                                                                                                                          caccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggc           >  1:1957522/3‑150 (MQ=255)
                                                                                                                                           tCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg          <  2:623026/150‑1 (MQ=255)
                                                                                                                                           tCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg          >  2:664853/1‑150 (MQ=255)
                                                                                                                                           accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg          >  1:447789/2‑150 (MQ=255)
                                                                                                                                           accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcg          >  1:121089/2‑150 (MQ=255)
                                                                                                                                           accAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGACGCCAGTTTTTCCGCCAGCACTTTggcg          >  2:1353705/2‑150 (MQ=255)
                                                                                                                                            ccAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcgg         <  1:1962733/150‑1 (MQ=255)
                                                                                                                                            ccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcgg         <  2:150339/150‑1 (MQ=255)
                                                                                                                                             cagcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                        >  2:743080/1‑102 (MQ=255)
                                                                                                                                              agcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCt       <  2:359606/150‑1 (MQ=255)
                                                                                                                                               gcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTa      >  1:215255/1‑150 (MQ=255)
                                                                                                                                               gcagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACATTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTa      <  2:1131018/150‑1 (MQ=255)
                                                                                                                                                cagTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTaa     <  1:1751430/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAATC  >  NC_000913/4364985‑4365281

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: