Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,133 G→C 8.0% A75G (GCA→GGA) ‡ cutA ← copper binding protein CutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,1330GC 8.0% 415.8 / 21.3 125A75G (GCA→GGA) ‡cutAcopper binding protein CutA
Reads supporting (aligned to +/- strand):  ref base G (59/56);  new base C (6/4);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 127 total aligned reads displayed.

CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAATCC  >  NC_000913/4364984‑4365282
                                                                                                                                                     |                                                                                                                                                     
   gCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGCTGCGTTGAGCCATGAGTGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAGCCAGAAGTTCGGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGcgt                                                                                                                                                    >  2:893810/1‑148 (MQ=255)
   gCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCct                                                                                                                                                    >  2:394029/1‑150 (MQ=255)
    cAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctg                                                                                                                                                   >  1:417382/1‑150 (MQ=255)
     aaaaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgc                                                                                                                                                  <  2:1336004/150‑1 (MQ=255)
      aaaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgct                                                                                                                                                 <  2:1639674/150‑1 (MQ=255)
        aaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTCCAGGCACTGCTGGAATGCctgaagt                                                                                                                                               <  2:609360/150‑5 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                             <  2:519956/150‑1 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                             <  2:703458/150‑1 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                             >  2:39723/1‑150 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                             <  2:581288/150‑1 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaa                                                                                                                                             <  2:1312245/150‑3 (MQ=255)
          cGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGGTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTg                                                                                                                                             <  1:194512/150‑1 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgag                                                                                                                                           <  1:779135/150‑1 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgag                                                                                                                                           >  1:663858/1‑150 (MQ=255)
            gAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGGCACTGCTggaatg                                                                                                                                           <  1:1254522/150‑5 (MQ=255)
               gTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATa                                                                                                                                        >  2:1204291/1‑150 (MQ=255)
                 gctgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACg                                                                                                                                      >  1:1389456/1‑150 (MQ=255)
                     cAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAg                                                                                                                                  <  1:1315742/150‑1 (MQ=255)
                         gTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtttt                                                                                                                              <  1:167113/150‑1 (MQ=255)
                          tAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGttttt                                                                                                                             <  2:1072192/150‑1 (MQ=255)
                            gCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaa                                                                                                                           <  2:453435/150‑1 (MQ=255)
                            gCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaa                                                                                                                           >  1:1318711/1‑150 (MQ=255)
                            gCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaa                                                                                                                           <  2:766894/150‑1 (MQ=255)
                              aGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTaaaa                                                                                                                         <  2:940265/150‑1 (MQ=255)
                                 aTCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatca                                                                                                                      <  2:1110331/150‑1 (MQ=255)
                                   cAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAatcatc                                                                                                                    <  2:137202/150‑1 (MQ=255)
                                    aGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCt                                                                                                                   >  1:636017/1‑150 (MQ=255)
                                    aGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCt                                                                                                                   <  2:1478677/150‑1 (MQ=255)
                                     gCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg                                                                                                                  <  1:1448656/150‑1 (MQ=255)
                                     gCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg                                                                                                                  >  1:978529/1‑150 (MQ=255)
                                      cGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGc                                                                                                                 <  2:434588/150‑1 (MQ=255)
                                      cGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGc                                                                                                                 <  1:823797/150‑1 (MQ=255)
                                            gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCg                                                                                                           <  1:1301952/150‑1 (MQ=255)
                                              tGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTa                                                                                                         >  2:1008361/1‑150 (MQ=255)
                                                cGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTAtt                                                                                                       <  1:1562751/150‑1 (MQ=255)
                                                      gCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCt                                                                                                 <  2:663858/150‑1 (MQ=255)
                                                        cATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTcc                                                                                               <  1:1182491/150‑1 (MQ=255)
                                                         aTGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCa                                                                                              >  1:478120/1‑150 (MQ=255)
                                                               ggTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTAc                                                                                        <  2:1195986/150‑1 (MQ=255)
                                                                 tAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCt                                                                                      >  2:79108/1‑150 (MQ=255)
                                                                  aaTCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTACGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCtt                                                                                     >  2:1238885/1‑150 (MQ=255)
                                                                   aTCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTc                                                                                    <  1:910386/150‑1 (MQ=255)
                                                                    tCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTcc                                                                                   >  1:965391/1‑150 (MQ=255)
                                                                     cTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTccc                                                                                  <  1:79108/150‑1 (MQ=255)
                                                                      tgtgTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAAATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCa                                                                                 >  1:836538/1‑150 (MQ=255)
                                                                         gtCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTa                                                                              >  1:971175/1‑150 (MQ=255)
                                                                          tctcCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTaa                                                                             >  1:583713/1‑150 (MQ=255)
                                                                             ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                          <  2:517919/150‑1 (MQ=255)
                                                                             ccGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTAGTATTCTTGCTCCAGCTTACCTTCCCAGTAATag                                                                          <  2:395433/150‑1 (MQ=255)
                                                                                tgtgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagaga                                                                       <  2:257703/150‑1 (MQ=255)
                                                                                  tgtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagaga                                                                     >  2:310232/1‑150 (MQ=255)
                                                                                   gtAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATagagagag                                                                    >  2:677822/1‑150 (MQ=255)
                                                                                    tAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAgg                                                                   >  2:1614985/1‑150 (MQ=255)
                                                                                    tAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAgg                                                                   >  2:505599/1‑150 (MQ=255)
                                                                                    tAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAgg                                                                   >  1:1481805/1‑150 (MQ=255)
                                                                                        aGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgc                                                               <  1:28483/150‑1 (MQ=255)
                                                                                         ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                              >  2:1322047/1‑150 (MQ=255)
                                                                                         ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                              >  2:419975/1‑150 (MQ=255)
                                                                                         ggTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGACATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                              >  1:740260/1‑150 (MQ=255)
                                                                                           tAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGcc                                                            >  1:1230984/1‑150 (MQ=255)
                                                                                              aaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                         <  1:1301284/150‑1 (MQ=255)
                                                                                              aaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGAAGCGCCggg                                                         >  2:816827/1‑150 (MQ=255)
                                                                                              aaCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTACGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCggg                                                         <  2:1274037/150‑1 (MQ=255)
                                                                                               aCCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCgggg                                                        >  1:739629/1‑150 (MQ=255)
                                                                                                 cAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGAt                                                      >  2:1641437/1‑150 (MQ=255)
                                                                                                   gAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCa                                                    <  2:1459571/150‑1 (MQ=255)
                                                                                                      gTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                              <  1:262272/137‑1 (MQ=255)
                                                                                                      gTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAgcg                                                              >  2:262272/1‑137 (MQ=255)
                                                                                                      gTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAgg                                                 <  2:1484294/150‑1 (MQ=255)
                                                                                                       ttCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGt                                                >  1:1135311/1‑150 (MQ=255)
                                                                                                        tCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTc                                               >  2:268250/1‑150 (MQ=255)
                                                                                                         ccGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcg                                              <  2:636017/150‑1 (MQ=255)
                                                                                                         ccGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcg                                              >  2:924188/1‑150 (MQ=255)
                                                                                                         ccGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcg                                              >  2:1166082/1‑150 (MQ=255)
                                                                                                          cGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTcgc                                             >  1:283584/1‑150 (MQ=255)
                                                                                                                 tGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgcc                                      >  1:502590/1‑150 (MQ=255)
                                                                                                                  gATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccg                                     <  2:288150/150‑1 (MQ=255)
                                                                                                                  gATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccg                                     <  2:1518254/150‑1 (MQ=255)
                                                                                                                   atatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCCTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgc                                    <  1:148498/150‑1 (MQ=255)
                                                                                                                    tatGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCATATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGgccgcc                                   >  1:192868/1‑150 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                                 >  1:1089800/1‑150 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                                 <  1:1362134/150‑1 (MQ=255)
                                                                                                                      tGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAg                                 >  1:826912/1‑150 (MQ=255)
                                                                                                                       ggATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGt                                <  1:1204291/150‑1 (MQ=255)
                                                                                                                        gatgatGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGtt                               <  1:1611978/150‑1 (MQ=255)
                                                                                                                              gagaCTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCg                         >  2:1292561/1‑150 (MQ=255)
                                                                                                                                  cTTAAGACATTCCAGCAGTGACTGATGGTGAGTTACTGTAGTTTTTAAAATAATCTGAAATTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAg                     >  2:918696/1‑150 (MQ=255)
                                                                                                                                   ttCAGGCATTCCAGCAGTGCCTGCTTGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGGGCAGGCCGCCAGTTTTTCCGCCAGc                    <  1:508179/150‑1 (MQ=255)
                                                                                                                                   ttCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGc                    <  1:1008361/150‑1 (MQ=255)
                                                                                                                                     cAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCAc                  >  1:966001/1‑150 (MQ=255)
                                                                                                                                       ggCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                       <  1:419648/143‑1 (MQ=255)
                                                                                                                                       ggCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGcc                       >  2:419648/1‑143 (MQ=255)
                                                                                                                                         ctcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTg              >  1:494651/5‑150 (MQ=255)
                                                                                                                                         cATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTg              >  2:438174/1‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg             >  1:422667/4‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg             >  1:1199843/4‑150 (MQ=255)
                                                                                                                                          tcaccAGCAGGCACTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg             >  2:1263436/4‑150 (MQ=255)
                                                                                                                                          aTTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTgg             >  1:885263/1‑150 (MQ=255)
                                                                                                                                             ccAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTggcgg          <  1:1366837/150‑1 (MQ=255)
                                                                                                                                                     gCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAATcc  >  2:643408/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTCGTATTCTTGCTCCAGCTTACCTTCCCAGTAATAGAGAGAGGTAGCGCCGGGGATCAAGGTCGCGCAGGCCGCCAGTTTTTCCGCCAGCACTTTGGCGGCTAAATCC  >  NC_000913/4364984‑4365282

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: