Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,121 G→A 6.8% L81L (CTG→CTA) ‡ PP_4263 → cytochrome biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,850,1210GA 6.8% 604.4 / 22.8 191L81L (CTG→CTA) ‡PP_4263cytochrome biogenesis protein
Reads supporting (aligned to +/- strand):  ref base G (91/87);  new base A (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.97e-01

Only 100 of 208 total aligned reads displayed.

GCACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTCCAG  >  NC_002947/4849972‑4850270
                                                                                                                                                     |                                                                                                                                                     
gcACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCAGTGTCAGCGATGCCTGCGCCGGCCGGGTGCGGGGCCTGGCCGGCGGGGCCGTGGCCAGCAGCCCGGCAgcgctg                                                                                                                                                       >  2:590539/1‑150 (MQ=255)
    tAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                   >  2:1078170/1‑150 (MQ=255)
    tAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                   >  2:3219138/1‑150 (MQ=255)
    tAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                   >  1:3910044/1‑150 (MQ=255)
    tAACATTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgc                                                                                                                                                   >  2:2760797/1‑150 (MQ=255)
      aCCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAgcgctggcgctg                                                                                                                                                 <  2:3726153/150‑1 (MQ=255)
       ccTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGc                                                                                                                                                <  1:2878899/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:3802795/150‑1 (MQ=255)
         ttGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCCAGCGCTGccg                                                                                                                                              <  2:109085/150‑3 (MQ=255)
         ttGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGACCGGCAGCGCCAGCGCTGccg                                                                                                                                              <  1:199501/150‑3 (MQ=255)
             ccATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGggtgg                                                                                                                                          <  2:2487277/150‑1 (MQ=255)
                  ccccccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTcgcg                                                                                                                                     <  1:2473275/150‑1 (MQ=255)
                   cccccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCgg                                                                                                                                    >  1:3084900/1‑150 (MQ=255)
                     cccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCa                                                                                                                                  >  1:1238717/1‑150 (MQ=255)
                     cccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCa                                                                                                                                  >  1:2857248/1‑150 (MQ=255)
                     cccGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCa                                                                                                                                  >  2:27909/1‑150 (MQ=255)
                        gAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActg                                                                                                                               >  1:1036978/1‑150 (MQ=255)
                            aacGCGGCTAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCTCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGTCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                           <  2:1013318/150‑1 (MQ=255)
                            aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                           <  2:3112875/150‑1 (MQ=255)
                            aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                           >  1:1759928/1‑150 (MQ=255)
                            aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                           <  2:1492010/150‑1 (MQ=255)
                            aacGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgc                                                                                                                           >  1:400644/1‑150 (MQ=255)
                             acGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                          >  2:4020870/1‑150 (MQ=255)
                              cgcgGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgctg                                                                                                                         >  2:1352642/1‑150 (MQ=255)
                               gcgGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGa                                                                                                                        >  1:3583806/1‑150 (MQ=255)
                                 ggCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGAtt                                                                                                                      <  2:3866416/150‑1 (MQ=255)
                                  gCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTg                                                                                                                     >  1:4141667/1‑150 (MQ=255)
                                    aaGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGcc                                                                                                                   >  1:2315424/1‑150 (MQ=255)
                                    aaGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGcc                                                                                                                   >  1:1455234/1‑150 (MQ=255)
                                       cgcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATg                                                                                                                <  1:3560395/150‑1 (MQ=255)
                                       cgcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATg                                                                                                                <  1:3560272/150‑1 (MQ=255)
                                        gcctgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATgg                                                                                                               <  1:1882525/150‑1 (MQ=255)
                                           tgcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGcc                                                                                                            >  2:1284254/1‑150 (MQ=255)
                                            gcgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCt                                                                                                           <  1:3306969/150‑1 (MQ=255)
                                             cgcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTg                                                                                                          <  1:1142834/150‑1 (MQ=255)
                                              gcctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGt                                                                                                         <  1:1102537/150‑1 (MQ=255)
                                                ctgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTAc                                                                                                       >  1:4345801/1‑150 (MQ=255)
                                                 tgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTAcc                                                                                                      >  1:647850/1‑150 (MQ=255)
                                                 tgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTAcc                                                                                                      >  1:708303/1‑150 (MQ=255)
                                                 tgTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTAcc                                                                                                      >  1:2978512/1‑150 (MQ=255)
                                                  gTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCt                                                                                                     >  2:3212911/1‑150 (MQ=255)
                                                     gCTGGCGTACAACCTGGGGCGGATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTTGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGc                                                                                                  <  1:724347/150‑1 (MQ=255)
                                                       tGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCg                                                                                                >  2:675503/1‑150 (MQ=255)
                                                       tGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCg                                                                                                >  2:4289358/1‑150 (MQ=255)
                                                       tGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCg                                                                                                <  1:478444/150‑1 (MQ=255)
                                                             aCAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCtggt                                                                                          >  1:2808103/1‑150 (MQ=255)
                                                             aCAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCtggt                                                                                          <  2:1161784/150‑1 (MQ=255)
                                                                 ccTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGc                                                                                                             <  2:3564616/127‑1 (MQ=255)
                                                                 ccTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGc                                                                                                             >  1:3564684/1‑127 (MQ=255)
                                                                    ggggCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCgg                                                                                   >  1:2636618/1‑150 (MQ=255)
                                                                        cgcATTTTCAGCTATGCCTGCGACGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCActgctgct                                                                                                                          >  2:935077/1‑107 (MQ=255)
                                                                         gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                              <  1:4246290/150‑1 (MQ=255)
                                                                         gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                              >  2:2184614/1‑150 (MQ=255)
                                                                         gcATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGa                                                                              <  2:3504001/150‑1 (MQ=255)
                                                                           aTTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGAcc                                                                            >  2:2959951/1‑150 (MQ=255)
                                                                                 aGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATc                                                                      >  1:278287/1‑150 (MQ=255)
                                                                                  gCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGTCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCg                                                                     <  1:209434/150‑1 (MQ=255)
                                                                                    tATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAg                                                                   <  2:2516519/150‑1 (MQ=255)
                                                                                     aTGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAgg                                                                  <  2:2565115/150‑1 (MQ=255)
                                                                                         cTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACt                                                              >  2:849496/1‑150 (MQ=255)
                                                                                         cTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACt                                                              >  2:4193753/1‑150 (MQ=255)
                                                                                           gcgcCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTgg                                                            <  1:590539/150‑1 (MQ=255)
                                                                                                  ccTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCgg                                                     >  2:1605786/1‑150 (MQ=255)
                                                                                                   cTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCggg                                                    >  1:1429099/1‑150 (MQ=255)
                                                                                                   cTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCggg                                                    <  2:3578865/150‑1 (MQ=255)
                                                                                                    tGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGc                                                   >  1:3440643/1‑150 (MQ=255)
                                                                                                     gTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCt                                                  >  1:618049/1‑150 (MQ=255)
                                                                                                      ttGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCtg                                                 <  2:2653018/150‑1 (MQ=255)
                                                                                                      ttGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCtg                                                 <  2:4474855/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:2959951/150‑1 (MQ=255)
                                                                                                              ccTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCa                                         >  2:2839627/1‑150 (MQ=255)
                                                                                                                 ggCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCatat                                      <  1:4342400/150‑1 (MQ=255)
                                                                                                                    cGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCa                                   <  2:1470931/150‑1 (MQ=255)
                                                                                                                      gCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGc                                 >  1:1698054/1‑150 (MQ=255)
                                                                                                                         gggCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCgg                              >  1:4506530/1‑150 (MQ=255)
                                                                                                                          ggCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGt                             <  1:4467184/150‑1 (MQ=255)
                                                                                                                          ggCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGt                             >  1:1543405/1‑150 (MQ=255)
                                                                                                                            ccGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTgg                           >  1:4022059/1‑150 (MQ=255)
                                                                                                                            ccGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGGGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTgg                           >  1:1423091/1‑150 (MQ=255)
                                                                                                                              gTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCg                         >  2:2089622/1‑150 (MQ=255)
                                                                                                                                   cagcagCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTcgcg                    >  1:3379843/1‑150 (MQ=255)
                                                                                                                                   cagcagCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGGATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTcgcg                    >  2:4079083/1‑150 (MQ=255)
                                                                                                                                       agcCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtg                <  2:540344/150‑1 (MQ=255)
                                                                                                                                         cccGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgct              >  2:2011701/1‑150 (MQ=255)
                                                                                                                                         cccGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgct              >  2:3084546/1‑150 (MQ=255)
                                                                                                                                           ccgCAGCGCCAGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGTCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc            >  2:2331633/3‑150 (MQ=255)
                                                                                                                                           ccgCAGCGCCAGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc            >  1:812471/3‑150 (MQ=255)
                                                                                                                                           ccgCAGCGCCAGCGCTGCCGGTGGTCGCGGCACTGCTGCTGATTGCCa                                                                                                                  >  2:2694511/3‑48 (MQ=21)
                                                                                                                                           ccgCAGCGCCAGCGCTGCCGGTGGTCGCGGCACTGCTGCTGATTGCCa                                                                                                                  <  1:2694474/46‑1 (MQ=21)
                                                                                                                                           ccgCAGCGCCAGCGCTGCCGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc            >  2:1405857/3‑150 (MQ=255)
                                                                                                                                           cGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATtgctgc            <  1:3959464/150‑1 (MQ=255)
                                                                                                                                               aGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAgt        <  2:708303/150‑1 (MQ=255)
                                                                                                                                               aGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAgt        <  2:1759928/150‑1 (MQ=255)
                                                                                                                                                   ctggcgctgCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTcc    <  1:2339992/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACTAACCTTGGCCATCCCCCCCGAACAACGCGGCAAGCGCCTGCGCCTGTTGCTGGCGTACAACCTGGGGCGCATTTTCAGCTATGCCTGCGCCGGCCTGTTGCTGGGCCTGGCCGGCTGGGCCGTGGCCAGCAGCCCGGCAGCGCTGGCGCTGCGGGTGGTCGCGGCACTGCTGCTGATTGCCATGGGCCTGTACCTGGCCGGCTGGTGGAGCGGGCTGACCCGTATCGAGGCACTGGGCCGCGGGCTGTGGCGGCATATCCAGCCGGTGGCGTCGCGATTGCTGCCAGTGTCCAG  >  NC_002947/4849972‑4850270

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: