Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,671,701 T→C 5.3% V22A (GTG→GCG) ‡ 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,7010TC 5.3% 667.4 / 15.3 204V22A (GTG→GCG) ‡nuoNNADH‑quinone oxidoreductase subunit N
Reads supporting (aligned to +/- strand):  ref base T (99/94);  new base C (3/8);  total (102/102)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 221 total aligned reads displayed.

GACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAGT  >  NC_002947/4671552‑4671850
                                                                                                                                                     |                                                                                                                                                     
gACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggt                                                                                                                                                       <  2:603668/150‑1 (MQ=255)
  cacCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtgg                                                                                                                                                     <  2:789861/150‑1 (MQ=255)
    ccTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtg                                                                                                                                                   <  1:3436465/150‑1 (MQ=255)
     ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtgg                                                                                                                                                  <  1:2883832/150‑1 (MQ=255)
     ctctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGGAAGAGCGCTATGGAATTCACCACTCACCACTTCATCGCCGTGGCGCCGATGCTGAGCACCACCATCACCACggtggggg                                                                                                                                                  >  2:2522489/1‑150 (MQ=255)
      tctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACggtggtggt                                                                                                                                                 >  1:1391894/1‑150 (MQ=255)
       ctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGgtggtggtg                                                                                                                                                >  2:913573/1‑150 (MQ=255)
       ctGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGgtggtggtg                                                                                                                                                >  1:99124/1‑150 (MQ=255)
        tGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGa                                                                                                                                               <  1:2929442/150‑1 (MQ=255)
          ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCGTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATg                                                                                                                                             >  2:2346376/1‑150 (MQ=255)
          ccgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATg                                                                                                                                             <  1:1979006/150‑1 (MQ=255)
           cgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                            >  2:2323584/1‑150 (MQ=255)
           cgccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGc                                                                                                                                            <  1:2952956/150‑1 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                           <  1:529385/150‑1 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                           <  2:4088480/150‑1 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                           >  2:3435142/1‑150 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCt                                                                                                                                           >  2:2146723/1‑150 (MQ=255)
            gccACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATggt                                                                                                                                           <  1:4249194/150‑3 (MQ=255)
             ccaccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCAGCATCACCACCACCGTGGTGATGGTg                                                                                                                                          <  2:2823293/150‑1 (MQ=255)
             ccaccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACCACCGTGGTGATGGTg                                                                                                                                          <  1:3229676/150‑1 (MQ=255)
              caccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTgg                                                                                                                                         <  1:2545577/150‑1 (MQ=255)
              caccaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTgg                                                                                                                                         <  1:825423/150‑1 (MQ=255)
               accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGc                                                                                                                                        <  1:2579867/150‑1 (MQ=255)
               accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGc                                                                                                                                        >  2:2232804/1‑150 (MQ=255)
               accaTGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGc                                                                                                                                        >  2:1105871/1‑150 (MQ=255)
                     aGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGc                                                                                                                                  <  1:3802746/150‑1 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCCCCATCACCACGGTGGTGGTGATGCTGGTGATCGCCt                                                                                                                                <  1:1478634/150‑1 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGGTGGTGATcagca                                                                                                                                <  1:1713981/150‑5 (MQ=255)
                       tggtgTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCt                                                                                                                                >  2:1533285/1‑150 (MQ=255)
                          tgtgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGa                                                                                                                             >  2:2658001/1‑150 (MQ=255)
                           gtgCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGaa                                                                                                                            <  2:2752812/150‑1 (MQ=255)
                             gcggcagTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAgc                                                                                                                          <  1:720901/147‑1 (MQ=255)
                             gcagcaGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAgc                                                                                                                          <  1:3685084/150‑1 (MQ=255)
                                     ggCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCAc                                                                                                                  <  2:4210991/150‑1 (MQ=255)
                                      gCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACt                                                                                                                 >  1:4451660/1‑150 (MQ=255)
                                        tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                               >  2:4481963/1‑150 (MQ=255)
                                        tCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCg                                                                                                               >  2:3038425/1‑150 (MQ=255)
                                          ggATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCa                                                                                                             <  1:3835541/150‑1 (MQ=255)
                                            aTCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGa                                                                                                           <  1:1435189/150‑1 (MQ=255)
                                            aTCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGa                                                                                                           >  1:2153840/1‑150 (MQ=255)
                                             tCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAc                                                                                                          >  1:1204099/1‑150 (MQ=255)
                                             tCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAc                                                                                                          >  2:4320221/1‑150 (MQ=255)
                                              ccGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGAcc                                                                                                         >  1:4082591/1‑150 (MQ=255)
                                                 cTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTc                                                                                                      <  2:1738561/150‑1 (MQ=255)
                                                  tttCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTcc                                                                                                     >  1:677976/1‑150 (MQ=255)
                                                  tttCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTcc                                                                                                     >  2:3499522/1‑150 (MQ=255)
                                                       cTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                                >  1:1783354/1‑150 (MQ=255)
                                                       cTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCctgctg                                                                                                >  2:310312/1‑150 (MQ=255)
                                                                 ttCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTgg                                                                                      <  1:1258711/150‑1 (MQ=255)
                                                                   cGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGc                                                                                    <  2:900284/150‑1 (MQ=255)
                                                                        cGGTAAGAGCGCTATGGAAGTCACCACTCACCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGaa                                                                               <  2:149370/150‑1 (MQ=255)
                                                                         ggTAAGAGCGCTATGGAATTCACCACTCAACACTTCCTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAc                                                                              <  2:3223130/150‑1 (MQ=255)
                                                                          gTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAAcc                                                                             >  1:4355291/1‑150 (MQ=255)
                                                                               agCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGcc                                                                        <  1:28393/150‑1 (MQ=255)
                                                                                  gcTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctg                                                                     <  1:1499982/150‑1 (MQ=255)
                                                                                   cTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgc                                                                    <  2:4283134/150‑1 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                   <  2:4302780/150‑1 (MQ=255)
                                                                                    tATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCctgct                                                                   >  1:1855038/1‑150 (MQ=255)
                                                                                      tGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGt                                                                 >  1:4056766/1‑150 (MQ=255)
                                                                                      tGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGt                                                                 >  1:806299/1‑150 (MQ=255)
                                                                                             caccacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCt                                                          <  1:1899123/150‑1 (MQ=255)
                                                                                                cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGcc                                                       >  2:3330261/1‑150 (MQ=255)
                                                                                                cacTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGcc                                                       >  2:3945800/1‑150 (MQ=255)
                                                                                                  cTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCgg                                                     <  2:872281/150‑1 (MQ=255)
                                                                                                     aaCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcgc                                                  <  2:4129526/150‑1 (MQ=255)
                                                                                                     aaCACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcgc                                                  >  2:1383742/1‑150 (MQ=255)
                                                                                                        acTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGgcccgga                                               >  1:3940466/1‑145 (MQ=255)
                                                                                                         cTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                              <  2:876047/150‑1 (MQ=255)
                                                                                                          ttCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAg                                             <  1:1698468/150‑1 (MQ=255)
                                                                                                           tcctCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGCGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAgg                                            <  2:3678909/147‑1 (MQ=255)
                                                                                                             aTCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTc                                          >  1:3664122/1‑150 (MQ=255)
                                                                                                              tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGa                                                                                                                                               >  2:4203639/1‑48 (MQ=255)
                                                                                                              tCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcg                                         >  1:3138437/1‑150 (MQ=255)
                                                                                                                gCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGTAGGTcgcg                                       >  2:2419212/1‑150 (MQ=255)
                                                                                                                gCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTcgcg                                       >  2:1437786/1‑150 (MQ=255)
                                                                                                                gCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGAGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGACGGTcgcg                                       >  2:2674013/1‑150 (MQ=255)
                                                                                                                  aTTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGcgcc                                     >  2:2991708/1‑150 (MQ=255)
                                                                                                                        gcCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGc                               >  2:1572951/1‑150 (MQ=255)
                                                                                                                         ccGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCg                              <  1:2279789/150‑1 (MQ=255)
                                                                                                                             tGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTca                          <  1:3555344/150‑1 (MQ=255)
                                                                                                                                tGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACAt                       <  2:709873/150‑1 (MQ=255)
                                                                                                                                   tCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGc                    >  1:2521963/1‑150 (MQ=255)
                                                                                                                                     gccagcaTCACCACCACCGTGGTGATGGTGGTGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTg                  >  1:4187100/2‑150 (MQ=39)
                                                                                                                                     accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTg                  >  1:1084255/1‑150 (MQ=255)
                                                                                                                                     accaccATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTg                  >  1:854757/1‑150 (MQ=255)
                                                                                                                                       cagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                              >  1:664373/4‑120 (MQ=38)
                                                                                                                                       cagcaTCACCACCACCGTGGTGATGGTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGaa                                              <  2:664373/117‑1 (MQ=38)
                                                                                                                                       caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTgct                >  2:291878/1‑148 (MQ=255)
                                                                                                                                       caccaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGAt                >  1:925037/1‑150 (MQ=255)
                                                                                                                                        accaTCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATc               >  1:3271679/1‑150 (MQ=255)
                                                                                                                                         ccaTCACCACGCTGGGGGCGATGCGGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCa              <  1:3035659/150‑1 (MQ=255)
                                                                                                                                            tCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCa           <  2:1391894/150‑1 (MQ=255)
                                                                                                                                              accacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGACGGCGCTGAAGGTCGCGCCGCTGGCGGTCCCATCGCTGATCACCATc         >  2:25788/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa       <  1:3341658/150‑1 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa       >  1:1523321/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa       >  2:2385744/1‑150 (MQ=255)
                                                                                                                                                cacGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGa       >  1:1923731/1‑150 (MQ=255)
                                                                                                                                                 acGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGAc      >  2:3864906/1‑150 (MQ=255)
                                                                                                                                                    gtggtggtgATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAg   >  2:720882/1‑150 (MQ=255)
                                                                                                                                                     tggtggtgATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAGt  <  1:4001486/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACACCTCTGCCGCCACCATGAGTGGTGTGCAGCAGTGGCTCGGATCCGCTTTCACTCAACTCGCTTCGGCCCGGTAAGAGCGCTATGGAATTCACCACTCAACACTTCATCGCATTGGCGCCGATGCTGATCACCACCATCACCACGGTGGTGGTGATGCTGGCGATCGCCTGGAAGCGCAACCACTCGCAGACCTTCCTGCTGTCCACCGTGGGCCTGAACCTGGCCCTGCTGTCGCTCCTGCCGGCGCTGAAGGTCGCGCCGCTGGCGGTCACATCGCTGATCACCATCGACAAGT  >  NC_002947/4671552‑4671850

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: