Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,303,940 C→T 6.4% P178L (CCG→CTG)  cysS → cysteine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,303,9400CT 6.4% 216.8 / 2.6 74P178L (CCG→CTG) cysScysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base T (3/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGC  >  NC_002947/3303796‑3304088
                                                                                                                                                |                                                                                                                                                    
aCAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCg                                                                                                                                                 >  1:548045/1‑150 (MQ=255)
  aaaGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGAc                                                                                                                                               <  1:42642/150‑1 (MQ=255)
   aaGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACt                                                                                                                                              >  2:52446/1‑150 (MQ=255)
   aaGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACt                                                                                                                                              <  1:190130/150‑1 (MQ=255)
     gggTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTc                                                                                                                                            <  1:41069/150‑1 (MQ=255)
      ggTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCg                                                                                                                                           >  2:1052482/1‑150 (MQ=255)
        tACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                         <  1:1023115/150‑1 (MQ=255)
        tACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                         >  1:627259/1‑150 (MQ=255)
        tACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATAGAATTCCTGCGTATGTGTTCACTCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                         <  1:739870/150‑1 (MQ=255)
                     cAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCg                                                                                                                            <  2:758073/150‑1 (MQ=255)
                          aaTGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAg                                                                                                                       >  2:1101134/1‑150 (MQ=255)
                            tGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCCAGAGCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGcc                                                                                                                     >  2:114731/1‑150 (MQ=255)
                                   gtgtACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                         >  1:359639/1‑123 (MQ=255)
                                   gtgtACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTc                                                                                                                                         <  2:359639/123‑1 (MQ=255)
                                         tacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAg                                                                                                              <  2:237321/144‑1 (MQ=255)
                                         tacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAg                                                                                                              >  1:237321/1‑144 (MQ=255)
                                         tacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGc                                                                                                        >  2:1048763/1‑150 (MQ=255)
                                         tacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGc                                                                                                        >  2:791609/1‑150 (MQ=255)
                                           ccGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTg                                                                                                      <  2:548045/150‑1 (MQ=255)
                                           ccGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTg                                                                                                      <  1:607579/150‑1 (MQ=255)
                                                    cAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCTGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGcc                                                                                             >  1:270150/1‑150 (MQ=255)
                                                      aGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCTGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCAt                                                                                           <  2:270150/150‑1 (MQ=255)
                                                               gCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGaa                                                                                                                                 >  2:34047/1‑103 (MQ=255)
                                                               gCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGaa                                                                                                                                 <  1:34047/103‑1 (MQ=255)
                                                               gCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCgg                                                                                  >  1:214954/1‑150 (MQ=255)
                                                               gCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCgg                                                                                  >  2:1295199/1‑150 (MQ=255)
                                                                cTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCggg                                                                                 >  2:516101/1‑150 (MQ=255)
                                                                    ggCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCtgt                                                                                                                                     <  2:1247708/94‑1 (MQ=255)
                                                                    ggCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCtgt                                                                                                                                     >  1:1247708/1‑94 (MQ=255)
                                                                     gCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATc                                                                                                <  1:733880/130‑1 (MQ=255)
                                                                     gCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATc                                                                                                >  2:733880/1‑130 (MQ=255)
                                                                      cAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGccc                                                                           >  1:1132438/1‑150 (MQ=255)
                                                                         gCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCggg                                                                        <  2:214954/150‑1 (MQ=255)
                                                                            gTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                     <  1:1295199/150‑1 (MQ=255)
                                                                             gcgcgCAAGAGGATCGAAGACCTGCGCATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTgg                                                                    <  2:656614/149‑1 (MQ=255)
                                                                                   aaGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCggg                                                                                 >  2:226764/1‑131 (MQ=255)
                                                                                   aaGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCggg                                                                                 <  1:226764/131‑1 (MQ=255)
                                                                                     gagGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACAGCGa                                                            >  1:55597/1‑150 (MQ=255)
                                                                                       ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGc                                                                       <  1:584958/137‑1 (MQ=255)
                                                                                       ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGc                                                                       >  2:584958/1‑137 (MQ=255)
                                                                                        gATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTg                                                         >  1:823453/1‑150 (MQ=255)
                                                                                         aTCGAAGACCTGAGTATCGGTGTACTCATCGAGGTCGACGAAACCAAGCAGGACTTGCTCGACTTCGTCCTGTGGAAGGGCGTCAACCCGGGCGAGCCAATCTGGGAATCGCCATGGGTCCCGGGGGGCGCCGGCTGGCACTTCGAGTGc                                                        >  2:200465/1‑150 (MQ=255)
                                                                                          tCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGAc                                                                                                                                               >  2:415954/1‑62 (MQ=255)
                                                                                          tCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGAc                                                                                                                                               <  1:415954/62‑1 (MQ=255)
                                                                                             aaGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACAGCGAGTGCTCgg                                                    >  1:293652/1‑150 (MQ=255)
                                                                                                  cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCtgt                                                                                                                                     <  1:450627/64‑1 (MQ=255)
                                                                                                  cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCtgt                                                                                                                                     >  2:450627/1‑64 (MQ=255)
                                                                                                  cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGt                                                                                                                           >  1:907314/1‑74 (MQ=255)
                                                                                                  cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGt                                                                                                                           <  2:907314/74‑1 (MQ=255)
                                                                                                  cTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATg                                               >  1:588169/1‑150 (MQ=255)
                                                                                                   tGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                              >  2:98496/1‑150 (MQ=255)
                                                                                                   tGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                              <  1:1048763/150‑1 (MQ=255)
                                                                                                   tGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                              <  1:1101134/150‑1 (MQ=255)
                                                                                                   tGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                              <  1:98496/150‑1 (MQ=255)
                                                                                                                aCGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCt                                 >  1:819894/1‑150 (MQ=255)
                                                                                                                 cGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTg                                <  1:52446/150‑1 (MQ=255)
                                                                                                                    aTCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACAGCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGc                             <  1:114731/150‑1 (MQ=255)
                                                                                                                        aGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                         <  1:516101/150‑1 (MQ=255)
                                                                                                                          gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGcc                                                                                     <  1:322842/88‑1 (MQ=255)
                                                                                                                          gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGcc                                                                                     >  2:322842/1‑88 (MQ=255)
                                                                                                                          gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGc                       <  2:293652/150‑1 (MQ=255)
                                                                                                                           tCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCt                      >  1:741986/1‑150 (MQ=255)
                                                                                                                             gacgaAGTCGAGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTc                    >  2:1424897/1‑150 (MQ=255)
                                                                                                                             gacgaAGTCGAGCGGGTCCTGCTTGGCTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTc                    <  1:939183/150‑1 (MQ=255)
                                                                                                                             gacgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTc                    <  1:791609/150‑1 (MQ=255)
                                                                                                                               cgaAGCCAAGCAGTACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGa                  <  1:200465/150‑1 (MQ=255)
                                                                                                                               cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGa                  <  2:627259/150‑1 (MQ=255)
                                                                                                                                   gCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACAtt              <  1:1052482/150‑1 (MQ=255)
                                                                                                                                     cAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCa            <  2:588169/150‑1 (MQ=255)
                                                                                                                                      aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAc           >  1:1186830/1‑150 (MQ=255)
                                                                                                                                      aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAc           <  2:1132438/150‑1 (MQ=255)
                                                                                                                                      aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTAGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCAc           >  2:450672/1‑150 (MQ=255)
                                                                                                                                       aGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACg          >  1:1206604/1‑150 (MQ=255)
                                                                                                                                            gACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggc     <  2:1206604/150‑1 (MQ=255)
                                                                                                                                               ccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                            <  1:182189/92‑1 (MQ=255)
                                                                                                                                               ccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                            >  2:182189/1‑92 (MQ=255)
                                                                                                                                               ccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcggc  >  2:240862/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACAAAGGGTACGCCTACGCGCCAGGCAATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGC  >  NC_002947/3303796‑3304088

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