Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,120 T→C 6.8% L81P (CTG→CCG) ‡ PP_4263 → cytochrome biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,850,1200TC 6.8% 615.3 / 21.9 191L81P (CTG→CCG) ‡PP_4263cytochrome biogenesis protein
Reads supporting (aligned to +/- strand):  ref base T (91/87);  new base C (5/8);  total (96/95)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 208 total aligned reads displayed.

CGCACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTCCA  >  NC_002947/4849971‑4850269
                                                                                                                                                     |                                                                                                                                                     
cgcACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCGGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgct                                                                                                                                                       <  2:180815/150‑1 (MQ=255)
 gcACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCAGTGTCAGCGATGCCTGCGCCGGCCGGGTGCGGGGCCTGGCCGGCGGGGCCGTGGCCAGCAGCCCGGCAgcgctg                                                                                                                                                      >  2:590539/1‑150 (MQ=255)
     tAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                  >  2:1084861/1‑150 (MQ=255)
     tAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                  >  2:3219138/1‑150 (MQ=255)
     tAACATTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                  >  2:2760797/1‑150 (MQ=255)
       aCCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgctg                                                                                                                                                <  2:3726153/150‑1 (MQ=255)
        ccTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGc                                                                                                                                               >  2:3023234/1‑150 (MQ=255)
         cTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCCAGCGCTGcc                                                                                                                                              <  2:2903126/150‑2 (MQ=255)
          ttGGCCATCCCCCCCGAACAACGCGGCACGCGCCTGCGCCTGTTGCTGGCGTACACCCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCgg                                                                                                                                             <  2:748359/150‑1 (MQ=255)
          ttGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCgg                                                                                                                                             <  1:3802728/150‑1 (MQ=255)
          ttGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGACCGGCAGCGCCAGCGCTGccg                                                                                                                                             <  1:199501/150‑3 (MQ=255)
           tGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCCAGCGCTGCcgg                                                                                                                                            <  1:3535401/150‑1 (MQ=255)
              ccATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGggtgg                                                                                                                                         <  2:2487277/150‑1 (MQ=255)
                   ccccccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTcgcg                                                                                                                                    <  1:2473275/150‑1 (MQ=255)
                   ccccccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTcgcg                                                                                                                                    <  1:4116389/150‑1 (MQ=255)
                     ccccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGc                                                                                                                                  >  1:47309/1‑150 (MQ=255)
                      cccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCa                                                                                                                                 >  2:27909/1‑150 (MQ=255)
                      cccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCa                                                                                                                                 >  1:2857248/1‑150 (MQ=255)
                         gAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActg                                                                                                                              >  1:1036978/1‑150 (MQ=255)
                         gAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActg                                                                                                                              >  1:2874868/1‑150 (MQ=255)
                             aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                          >  1:1759928/1‑150 (MQ=255)
                             aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                          <  2:1013375/150‑1 (MQ=255)
                             aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                          <  2:3112875/150‑1 (MQ=255)
                              acGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                         >  2:4020870/1‑150 (MQ=255)
                              acGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                         <  2:1251285/150‑1 (MQ=255)
                              acGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGACCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                         <  1:739926/150‑1 (MQ=255)
                                gcgGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGa                                                                                                                       <  1:616467/150‑1 (MQ=255)
                                  ggCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGAtt                                                                                                                     <  2:3866416/150‑1 (MQ=255)
                                     aaGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGcc                                                                                                                  >  1:1455234/1‑150 (MQ=255)
                                        cgcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATg                                                                                                               <  1:3560374/150‑1 (MQ=255)
                                         gcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATgg                                                                                                              <  2:1590112/150‑1 (MQ=255)
                                         gcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATgg                                                                                                              <  1:1882525/150‑1 (MQ=255)
                                         gcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATgg                                                                                                              <  1:123848/150‑1 (MQ=255)
                                          cctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATggg                                                                                                             <  2:1442642/150‑1 (MQ=255)
                                            tgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGcc                                                                                                           >  2:1284277/1‑150 (MQ=255)
                                             gcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCt                                                                                                          >  1:1356624/1‑150 (MQ=255)
                                             gcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCt                                                                                                          <  1:2298093/150‑1 (MQ=255)
                                               gcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGt                                                                                                        <  1:1102537/150‑1 (MQ=255)
                                                 ctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTAc                                                                                                      >  1:4345697/1‑150 (MQ=255)
                                                   gTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCt                                                                                                    >  2:3212911/1‑150 (MQ=255)
                                                      gCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGc                                                                                                 <  1:724199/150‑1 (MQ=255)
                                                        tGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCg                                                                                               >  2:675503/1‑150 (MQ=255)
                                                        tGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCg                                                                                               >  2:4289358/1‑150 (MQ=255)
                                                            gTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCtg                                                                                           >  1:3895451/1‑150 (MQ=255)
                                                              aCAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCtggt                                                                                         >  1:2808103/1‑150 (MQ=255)
                                                              aCAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCtggt                                                                                         <  2:1161777/150‑1 (MQ=255)
                                                                 aCCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGa                                                                                      <  1:66344/150‑1 (MQ=255)
                                                                  ccTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGc                                                                                                            <  2:3564616/127‑1 (MQ=255)
                                                                  ccTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGc                                                                                                            >  1:3564616/1‑127 (MQ=255)
                                                                    tGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCg                                                                                   >  2:1660188/1‑150 (MQ=255)
                                                                     ggggCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCgg                                                                                  >  1:2636618/1‑150 (MQ=255)
                                                                      gggCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCggg                                                                                 >  1:2970716/1‑150 (MQ=255)
                                                                      gggCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCggg                                                                                 <  2:1391773/150‑1 (MQ=255)
                                                                         cgcATTTTCAGCTATGCCTGCGACGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                         >  2:935077/1‑107 (MQ=255)
                                                                         cgcATTTTCAGCTATGCCTGCGACGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                         <  1:935077/107‑1 (MQ=255)
                                                                          gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                             <  2:3504001/150‑1 (MQ=255)
                                                                          gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                             >  1:3179737/1‑150 (MQ=255)
                                                                          gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                             >  2:2184614/1‑150 (MQ=255)
                                                                            aTTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGAcc                                                                           >  2:2959951/1‑150 (MQ=255)
                                                                                tCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTa                                                                       <  1:2987832/150‑1 (MQ=255)
                                                                                  aGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATc                                                                     >  1:278287/1‑150 (MQ=255)
                                                                                   gCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGTCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCg                                                                    <  1:3403949/150‑1 (MQ=255)
                                                                                   gCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGTCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCg                                                                    <  1:209434/150‑1 (MQ=255)
                                                                                     tATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAg                                                                  <  1:3921436/150‑1 (MQ=255)
                                                                                     tATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAg                                                                  <  2:2516519/150‑1 (MQ=255)
                                                                                      aTGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAgg                                                                 <  2:2565115/150‑1 (MQ=255)
                                                                                          cTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACt                                                             >  2:4193753/1‑150 (MQ=255)
                                                                                          cTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACt                                                             >  2:849496/1‑150 (MQ=255)
                                                                                          cTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACt                                                             >  2:147682/1‑150 (MQ=255)
                                                                                            gcgcCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTgg                                                           <  1:4501659/150‑1 (MQ=255)
                                                                                            gcgcCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTgg                                                           <  1:3574981/150‑1 (MQ=255)
                                                                                                    cTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCggg                                                   >  1:1429099/1‑150 (MQ=255)
                                                                                                    cTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCggg                                                   <  2:3578865/150‑1 (MQ=255)
                                                                                                     tGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGc                                                  >  1:4262223/1‑150 (MQ=255)
                                                                                                       ttGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCtg                                                >  1:3339355/1‑150 (MQ=255)
                                                                                                       ttGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCtg                                                <  2:4474855/150‑1 (MQ=255)
                                                                                                       ttGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCtg                                                <  2:2653018/150‑1 (MQ=255)
                                                                                                           tGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTggc                                            >  2:3796027/1‑150 (MQ=255)
                                                                                                           ggggCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTggc                                            <  2:334957/149‑1 (MQ=255)
                                                                                                            gggCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTggcg                                           >  2:3465190/1‑150 (MQ=255)
                                                                                                             ggCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTggcgg                                          <  1:4170702/150‑1 (MQ=255)
                                                                                                               ccTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCa                                        >  2:2839627/1‑150 (MQ=255)
                                                                                                                       gCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGc                                >  1:1698054/1‑150 (MQ=255)
                                                                                                                           ggCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGt                            <  1:4418988/150‑1 (MQ=255)
                                                                                                                           ggCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGt                            >  1:1543405/1‑150 (MQ=255)
                                                                                                                               gTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCg                        >  2:2089622/1‑150 (MQ=255)
                                                                                                                                   ccAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTcgc                    >  1:2933934/1‑150 (MQ=255)
                                                                                                                                    cagcagCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGGATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTcgcg                   >  2:4079083/1‑150 (MQ=255)
                                                                                                                                     agcagcCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGa                  <  1:4020870/150‑1 (MQ=255)
                                                                                                                                        agcCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtg               <  2:540344/150‑1 (MQ=255)
                                                                                                                                          cccGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgct             >  2:2011701/1‑150 (MQ=255)
                                                                                                                                          cccGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgct             >  2:3084546/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctg            >  1:4949/1‑150 (MQ=255)
                                                                                                                                            ccgCAGCGCCAGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGTCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc           >  2:2331633/3‑150 (MQ=255)
                                                                                                                                            ccgCAGCGCCAGCGCTGCCGGTGGTCGCGGCACTGCTGCTGATTGCCa                                                                                                                 <  1:2694474/46‑1 (MQ=21)
                                                                                                                                            ccgCAGCGCCAGCGCTGCCGGTGGTCGCGGCACTGCTGCTGATTGCCa                                                                                                                 >  2:2694511/3‑48 (MQ=21)
                                                                                                                                            cGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc           <  1:3959561/150‑1 (MQ=255)
                                                                                                                                                aGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAgt       <  2:708303/150‑1 (MQ=255)
                                                                                                                                                  cgctggcgctgCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAgtgt     <  1:633884/150‑1 (MQ=255)
                                                                                                                                                    ctggcgctgCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTcc   <  1:2339992/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTCCA  >  NC_002947/4849971‑4850269

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: