Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,671,700 G→A 5.5% V22M (GTG→ATG) ‡ nuoN → NADH‑quinone oxidoreductase subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,671,7000GA 5.5% 661.2 / 16.2 204V22M (GTG→ATG) ‡nuoNNADH‑quinone oxidoreductase subunit N
Reads supporting (aligned to +/- strand):  ref base G (99/93);  new base A (3/8);  total (102/102)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAG  >  NC_002947/4671552‑4671849
                                                                                                                                                    |                                                                                                                                                     
gACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggt                                                                                                                                                      <  2:603668/150‑1 (MQ=255)
    ccTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  1:3436465/150‑1 (MQ=255)
    ccTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                  <  1:710847/150‑1 (MQ=255)
     ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                                 <  1:2883832/150‑1 (MQ=255)
     ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                                 <  2:1563363/150‑1 (MQ=255)
      tctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtggt                                                                                                                                                >  1:1391894/1‑150 (MQ=255)
      tctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtggt                                                                                                                                                >  1:4367080/1‑150 (MQ=255)
       ctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGgtggtggtg                                                                                                                                               >  2:913573/1‑150 (MQ=255)
          ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCGTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATg                                                                                                                                            >  2:2346376/1‑150 (MQ=255)
          ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATg                                                                                                                                            <  1:1979006/150‑1 (MQ=255)
           cgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                           <  1:4321795/150‑1 (MQ=255)
           cgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                           >  2:2323584/1‑150 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                          <  2:4088480/150‑1 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                          >  2:3435142/1‑150 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCCTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                          <  2:1358491/150‑1 (MQ=255)
             ccaccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATGGTg                                                                                                                                         <  1:4465866/150‑1 (MQ=255)
             ccaccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATGGTg                                                                                                                                         <  1:3229676/150‑1 (MQ=255)
               accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCTCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGACGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGc                                                                                                                                       <  2:2681786/150‑1 (MQ=255)
                    gAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCAGCATCACCACCACCGTGGTGATGGTGGTGATca                                                                                                                                  <  1:2089380/150‑2 (MQ=38)
                     aGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCGTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGc                                                                                                                                 <  1:3802772/150‑1 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCCACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCACCACTTCATCGCATTGGCGCCGAGGCTGATCACCACCCTCACCACGGTGGTGGTGATGCTGGCGATCGCCt                                                                                                                               <  2:2399474/150‑1 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGGTGGTGATcagca                                                                                                                               <  1:1713981/150‑5 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCt                                                                                                                               <  1:664267/150‑1 (MQ=255)
                          tgtgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGa                                                                                                                            >  2:2658001/1‑150 (MQ=255)
                              cagcagTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAgcg                                                                                                                        >  2:28393/1‑150 (MQ=255)
                                gcagTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCa                                                                                                                      >  2:1009767/1‑150 (MQ=255)
                                     ggCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCAc                                                                                                                 <  2:4210991/150‑1 (MQ=255)
                                        tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                              >  1:854786/1‑150 (MQ=255)
                                          ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCa                                                                                                            <  2:1435391/150‑1 (MQ=255)
                                          ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCa                                                                                                            <  1:3835541/150‑1 (MQ=255)
                                           gATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGACGCGCAACCACTCGCAg                                                                                                           >  2:3166146/1‑150 (MQ=255)
                                            aTCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGa                                                                                                          <  1:1435189/150‑1 (MQ=255)
                                             tCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAc                                                                                                         >  1:1204099/1‑150 (MQ=255)
                                             tCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAc                                                                                                         >  2:4320221/1‑150 (MQ=255)
                                              ccGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAcc                                                                                                        >  1:4082591/1‑150 (MQ=255)
                                                gCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCtt                                                                                                      >  1:336300/1‑150 (MQ=255)
                                                 cTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTc                                                                                                     <  1:2546533/150‑1 (MQ=255)
                                                             tCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCAcc                                                                                         >  2:3285750/1‑150 (MQ=255)
                                                                cTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTg                                                                                      >  1:414594/1‑150 (MQ=255)
                                                                 ttCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTgg                                                                                     <  1:1258711/150‑1 (MQ=255)
                                                                   cGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGc                                                                                   <  2:900284/150‑1 (MQ=255)
                                                                   cGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATGGGCGCCGATGCTGATCACCACCATCACCCCGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGc                                                                                   <  2:3094541/150‑1 (MQ=255)
                                                                     gCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCt                                                                                 <  2:2058895/150‑1 (MQ=255)
                                                                         ggTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAc                                                                             <  1:2949259/150‑1 (MQ=255)
                                                                          gTAAGAGCGCTATGGAATTCACCACTCAACACTTCTTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAcc                                                                            >  1:1440581/1‑150 (MQ=255)
                                                                            aaGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTg                                                                          <  1:842722/150‑1 (MQ=255)
                                                                             agagCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTgg                                                                         <  2:138749/150‑1 (MQ=255)
                                                                               agCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGcc                                                                       <  1:28393/150‑1 (MQ=255)
                                                                                  gcTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                                    <  1:1499982/150‑1 (MQ=255)
                                                                                   cTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgc                                                                   >  1:959837/1‑150 (MQ=255)
                                                                                   cTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgc                                                                   <  2:4283134/150‑1 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                  <  2:4302780/150‑1 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                  >  2:2534992/1‑150 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                  >  1:1855038/1‑150 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                  <  1:1742631/150‑1 (MQ=255)
                                                                                      tGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGt                                                                >  1:4056615/1‑150 (MQ=255)
                                                                                             caccacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCt                                                         <  1:1899123/150‑1 (MQ=255)
                                                                                               ccacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGc                                                       >  2:1489408/1‑150 (MQ=255)
                                                                                                cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGcc                                                      >  2:3330261/1‑150 (MQ=255)
                                                                                                cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGcc                                                      >  2:3945800/1‑150 (MQ=255)
                                                                                                 acTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCg                                                     <  1:731879/150‑1 (MQ=255)
                                                                                                  cTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCgg                                                    <  2:872281/150‑1 (MQ=255)
                                                                                                   tCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgc                                                   >  2:2232104/1‑150 (MQ=255)
                                                                                                   tCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgc                                                   >  2:3681008/1‑150 (MQ=255)
                                                                                                     aaCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcgc                                                 <  2:4129526/150‑1 (MQ=255)
                                                                                                        acTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcccgga                                              >  1:3940466/1‑145 (MQ=255)
                                                                                                         cTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                             <  2:876047/150‑1 (MQ=255)
                                                                                                          ttCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAg                                            <  1:1698468/150‑1 (MQ=255)
                                                                                                             aTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTc                                         >  1:3664122/1‑150 (MQ=255)
                                                                                                              tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGa                                                                                                                                              <  1:4203639/48‑1 (MQ=255)
                                                                                                              tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcg                                        >  1:3138437/1‑150 (MQ=255)
                                                                                                              tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcg                                        >  2:4204172/1‑150 (MQ=255)
                                                                                                                 cATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCgcgc                                     >  1:2543910/1‑150 (MQ=255)
                                                                                                                  aTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGcgcc                                    >  2:2991708/1‑150 (MQ=255)
                                                                                                                      gcgcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTg                                >  1:3005244/1‑150 (MQ=255)
                                                                                                                         ccGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCg                             <  1:2279789/150‑1 (MQ=255)
                                                                                                                           gATGCTGATCACCACCATCCCCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGACGGTCGCGCCGCTGGCGGt                           <  1:33220/150‑1 (MQ=255)
                                                                                                                             tGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTca                         <  2:2500576/150‑1 (MQ=255)
                                                                                                                             tGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTca                         <  1:3555344/150‑1 (MQ=255)
                                                                                                                                tGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACAt                      <  2:709873/150‑1 (MQ=255)
                                                                                                                                  aTCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCg                    <  1:676852/150‑1 (MQ=255)
                                                                                                                                    caccacCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCt                  <  2:3609073/150‑1 (MQ=255)
                                                                                                                                    caccacCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCt                  <  2:1610391/150‑1 (MQ=255)
                                                                                                                                     accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTg                 >  1:1084255/1‑150 (MQ=255)
                                                                                                                                       cagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                             <  2:664373/117‑1 (MQ=38)
                                                                                                                                       caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTgct               >  2:291878/1‑148 (MQ=255)
                                                                                                                                       caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGAt               >  2:2821542/1‑150 (MQ=255)
                                                                                                                                         ccaTCACCACGCTGGGGGCGATGCGGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCa             <  1:3035659/150‑1 (MQ=255)
                                                                                                                                          caTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCAc            >  2:4505064/1‑150 (MQ=255)
                                                                                                                                             caccacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcg                                        >  2:806608/1‑119 (MQ=255)
                                                                                                                                             caccacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcg                                        <  1:806608/119‑1 (MQ=255)
                                                                                                                                              accacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGACGGCGCTGAAGGTCGCGCCGCTGGCGGTCCCATCGCTGATCACCATc        >  2:25788/1‑150 (MQ=255)
                                                                                                                                               ccacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCg       >  1:3725673/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAAGCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa      >  1:2763886/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTg                                >  2:2239326/1‑124 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa      >  1:1523321/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa      >  1:1923731/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa      <  2:1515013/150‑1 (MQ=255)
                                                                                                                                                 acGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGAc     >  2:3864906/1‑150 (MQ=255)
                                                                                                                                                    gtggtggtgATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAg  >  2:720882/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAG  >  NC_002947/4671552‑4671849

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