Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,162,068 T→C 5.6% L299P (CTG→CCG)  PP_2776 → homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,162,0680TC 5.6% 268.7 / 3.8 89L299P (CTG→CCG) PP_2776homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base C (2/3);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

AACGATACGCTGCTCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAGCGAAGTGCGC  >  NC_002947/3161923‑3162217
                                                                                                                                                 |                                                                                                                                                     
aaCGATACGCTGCTCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATctgcaa                                                                                                                                                   <  2:2077791/150‑1 (MQ=255)
    aTACGCTGCTCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATctgcaagcct                                                                                                                                               <  2:1642610/150‑1 (MQ=255)
         ctgctCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGCGCCCCTCCCCCCCCTCGACGCGCCTctgcacccctgcaag                                                                                                                                          >  2:1853620/1‑150 (MQ=255)
            ctCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAg                                                                                                                                       <  2:2022510/150‑1 (MQ=255)
             tCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgg                                                                                                                                      >  1:629193/1‑150 (MQ=255)
               gATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgggg                                                                                                                                    >  2:2026406/1‑150 (MQ=255)
                 tatTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGAt                                                                                                                                  >  2:185117/1‑150 (MQ=255)
                  atTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATg                                                                                                                                 <  1:185117/150‑1 (MQ=255)
                    tcttGCCGACCTGGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGa                                                                                                                               <  2:602738/147‑1 (MQ=255)
                        gCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATaa                                                                                                                           >  2:1229994/1‑150 (MQ=255)
                                  ggCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACaga                                                                                                                 <  1:1874953/150‑1 (MQ=255)
                                  ggCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACaga                                                                                                                 <  1:79305/150‑1 (MQ=255)
                                   gCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgg                                                                                                                                      <  1:651710/128‑1 (MQ=255)
                                   gCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgg                                                                                                                                      >  2:651710/1‑128 (MQ=255)
                                    cccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAgagg                                                                                                               >  1:1402173/1‑150 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAATGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAgagga                                                                                                              >  1:1496267/1‑150 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATaa                                                                                                                           <  1:1461337/137‑1 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATaa                                                                                                                           >  2:1461337/1‑137 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAAc                                                                                                                    >  1:1478908/1‑144 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAAc                                                                                                                    <  2:1478908/144‑1 (MQ=255)
                                     ccGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAgagga                                                                                                              >  2:167574/1‑150 (MQ=255)
                                         actgtaacTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGcc                                                                                                          >  2:803648/8‑150 (MQ=255)
                                            gTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCcgcg                                                                                                       >  2:2386803/1‑150 (MQ=255)
                                                 tGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAACCAGAGGAGGCCGCGCTTgg                                                                                                  <  2:2436325/150‑1 (MQ=255)
                                                 tGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTgg                                                                                                  >  2:1471543/1‑150 (MQ=255)
                                                  gCAGTGGGCGCAGCAGTGGGTCGAGGCCGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTNGGc                                                                                                 <  2:1496267/150‑1 (MQ=255)
                                                       gggCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAgg                                                                                            <  1:2026406/150‑1 (MQ=255)
                                                             agcagTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGttt                                                                                      >  2:289987/1‑150 (MQ=255)
                                                                   gggTCGAGGAAGTGGCTGACATTGGGGGCGGCTGCGGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGATGAGGCCGCGCTTGGCTAATTCTGTTTCTggcg                                                                                <  2:2333701/150‑1 (MQ=255)
                                                                    ggTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTggcgg                                                                               >  1:1089194/1‑150 (MQ=255)
                                                                     gTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAGGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTggcggc                                                                              >  1:160342/1‑150 (MQ=255)
                                                                      tCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCt                                                                             >  2:1354262/1‑150 (MQ=255)
                                                                      tCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCt                                                                             >  1:1508450/1‑150 (MQ=255)
                                                                        gAGGAAGGGGCTGACATTGTGGGCGCGCACTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTgg                                                                           >  1:1244827/1‑150 (MQ=255)
                                                                         aGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGa                                                                          >  2:1999441/1‑150 (MQ=255)
                                                                          ggAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAg                                                                         <  1:925344/150‑1 (MQ=255)
                                                                           gAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGc                                                                        >  2:186644/1‑150 (MQ=255)
                                                                             aGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCgg                                                                      <  2:505334/150‑1 (MQ=255)
                                                                              ggggCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCggg                                                                     <  2:374477/150‑1 (MQ=255)
                                                                              ggggCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCggg                                                                     <  1:2318208/150‑1 (MQ=255)
                                                                                ggCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCg                                                                   >  2:1260409/1‑150 (MQ=255)
                                                                                  cTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAGGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGtt                                                                 <  2:160342/150‑1 (MQ=255)
                                                                                  cTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGNGtt                                                                 <  2:1508450/150‑1 (MQ=255)
                                                                                   tGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGCAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTg                                                                >  1:659242/1‑150 (MQ=255)
                                                                                      cATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGt                                                             <  1:1443029/150‑1 (MQ=255)
                                                                                       aTTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTg                                                            >  2:696468/1‑150 (MQ=255)
                                                                                       aTTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTg                                                            >  2:1068055/1‑150 (MQ=255)
                                                                                         tgtgGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAg                                                          <  1:116471/150‑1 (MQ=255)
                                                                                           tgGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCgcgc                                                                                                      <  1:1782264/104‑1 (MQ=255)
                                                                                           tgGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCgcgc                                                                                                      >  2:1782264/1‑104 (MQ=255)
                                                                                                ggcTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCg                                                   <  1:803648/150‑1 (MQ=255)
                                                                                                ggcTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCg                                                   <  1:1999441/150‑1 (MQ=255)
                                                                                                  ctgctgCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGgt                                                 <  1:186644/150‑1 (MQ=255)
                                                                                                  ctgctgCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGgt                                                 <  1:1229994/150‑1 (MQ=255)
                                                                                                    gctgcGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGgtgt                                               <  1:2427343/150‑1 (MQ=255)
                                                                                                    gctgcGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGgtgt                                               <  2:2338298/150‑1 (MQ=255)
                                                                                                      tgcGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGa                                             <  1:1532787/150‑1 (MQ=255)
                                                                                                      tgcGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGAAGGGGATGCAATAATCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGa                                             <  1:2386803/150‑1 (MQ=255)
                                                                                                            aTCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACa                                       >  1:1450455/1‑150 (MQ=255)
                                                                                                             tCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgg                                                                                                                                      >  2:629716/1‑54 (MQ=255)
                                                                                                             tCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAgg                                                                                                                                      <  1:629716/54‑1 (MQ=255)
                                                                                                             tCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAg                                      >  2:2447134/1‑150 (MQ=255)
                                                                                                               ggCCCGGAGCACATCGCCCGACTTGACGAGCATCCGCAAGCCTGCAACCTGCGTATGTAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGcc                                    <  2:1244827/150‑1 (MQ=255)
                                                                                                                gCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCcg                                   >  2:1527482/1‑150 (MQ=255)
                                                                                                                gCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCcg                                   >  2:2316339/1‑150 (MQ=255)
                                                                                                                  ccGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCcgcg                                 >  1:978662/1‑150 (MQ=255)
                                                                                                                      aGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGTTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCgg                             <  1:1022821/150‑1 (MQ=255)
                                                                                                                       gCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCggg                            >  1:1564496/1‑150 (MQ=255)
                                                                                                                          caTCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCAAAATCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGa                         <  1:167574/150‑1 (MQ=255)
                                                                                                                           aTCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAg                        >  1:914376/1‑150 (MQ=255)
                                                                                                                            tCGCCCGACTTGACGAGCATCCGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGcc                                    >  2:1588066/1‑137 (MQ=255)
                                                                                                                            tCGCCCGACTTGACGAGCATCCGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGcc                                    >  2:386055/1‑137 (MQ=255)
                                                                                                                            tCGCCCGACTTGACGAGCATCCGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGcc                                    <  1:1588066/137‑1 (MQ=255)
                                                                                                                            tCGCCCGACTTGACGAGCATCCGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGcc                                    <  1:386055/137‑1 (MQ=255)
                                                                                                                               cccGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTCGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAg                    >  2:1686103/1‑150 (MQ=255)
                                                                                                                                ccGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGa                   >  2:1381557/1‑150 (MQ=255)
                                                                                                                                ccGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGa                   >  1:1538881/1‑150 (MQ=255)
                                                                                                                                ccGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGCCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGa                   <  1:872223/150‑1 (MQ=255)
                                                                                                                                  gACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGtt                                                                 <  2:1459820/102‑1 (MQ=255)
                                                                                                                                  gACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGtt                                                                 >  1:1459820/1‑102 (MQ=255)
                                                                                                                                  gACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGGTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACg                 >  2:295543/1‑150 (MQ=255)
                                                                                                                                    gttGGCGAGCATTTGGAAGGCTGGAAGCAGGGTATGGTATAAGCGAAACAGAGGAGGCCGCGGCTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACggg               <  1:1853620/149‑1 (MQ=255)
                                                                                                                                     ttGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGc              >  2:1536125/1‑150 (MQ=255)
                                                                                                                                      tGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACTGGCa             <  1:1319552/150‑1 (MQ=255)
                                                                                                                                      tGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCa             <  1:1260409/150‑1 (MQ=255)
                                                                                                                                      tGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCa             <  1:289987/150‑1 (MQ=255)
                                                                                                                                       gACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAg            <  2:1450455/150‑1 (MQ=255)
                                                                                                                                       gACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAg            <  1:1381557/150‑1 (MQ=255)
                                                                                                                                        aCGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTCATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAGc           <  2:457123/150‑1 (MQ=255)
                                                                                                                                             cATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGttt                                                                                      <  1:932191/70‑1 (MQ=255)
                                                                                                                                             cATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGttt                                                                                      >  2:932191/1‑70 (MQ=255)
                                                                                                                                              aTCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCAGGCGTTGCCTCACCTTGCCTGTCGAATCGACAGCCGCGTCGGGTGAGGCAGAGGGGCAGCGACGTg     >  2:392399/1‑150 (MQ=255)
                                                                                                                                                 tgcaagcctgcaagcAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAGCGAAGTgcgc  <  2:1847053/150‑1 (MQ=255)
                                                                                                                                                 tgcaagcctgcaagcAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAGCGAAGTgcgc  >  2:2058842/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AACGATACGCTGCTCGATATTCGTGCCGACCTAGGGCCCGAGGCGTACCTGCAGTGGGCGCAGCAGTGGGTCGAGGAAGGGGCTGACATTGTGGGCGGCTGCTGCGGCATCGGCCCGGAGCACATCGCCCGACTTGACGAGCATCTGCAAGCCTGCAAGCAGGGGATGGAATAAGCGAAACAGAGGAGGCCGCGCTTGGCGAAGGCTGTTTCTGGCGGCTGGAGCGGGCGTTGCGTGAGCTTGGCGGTGTGATCGACAGCCGCGTCGGGTGAGGCAGACGGGCAGCGAAGTGCGC  >  NC_002947/3161923‑3162217

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