Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 388,725 T→A 7.2% L198Q (CTG→CAG)  soxB → sarcosine oxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947388,7250TA 7.2% 377.7 / 13.3 122L198Q (CTG→CAG) soxBsarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (60/53);  new base A (4/5);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGG  >  NC_002947/388576‑388874
                                                                                                                                                     |                                                                                                                                                      
gAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCt                                                                                                                                                        >  2:110832/1‑150 (MQ=255)
  acacTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCGTGCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                      <  1:1966983/150‑1 (MQ=255)
  acacTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                      <  1:2086429/150‑1 (MQ=255)
   cacTGCCCAGGTTGCGGCCGAACTCACGTAGATGTACTCCTCTAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTggg                                                                                                                                                     <  2:1469230/150‑1 (MQ=255)
         ccAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTAGTGGCGTTGCCCGTCACTATACCGTTGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCgg                                                                                                                                               <  1:2242904/150‑2 (MQ=255)
          cAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGAc                                                                                                                                              <  1:110832/150‑1 (MQ=255)
           aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCggca                                                                                                                                             <  2:1647902/150‑4 (MQ=255)
           aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCggca                                                                                                                                             <  2:1246769/150‑4 (MQ=255)
            ctcaacGGCAGAAATCCCGTACCTGGCCTGCTCCAACAACACCCGCTACCCGATCCTCGGCGCCCCCGTACAGCGCCGTCGTGGCGTTGCCCGCCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCCGCCGACCCCCTGGGCGTCGACCt                                                                                                                                            >  2:532682/6‑150 (MQ=255)
                gCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGAGCTGTGGTGGAGGTGCCCGCAAAGAGGTCGTTGCCTGGGGCTCCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATc                                                                                                                                        <  1:1501240/150‑1 (MQ=255)
                    ccGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGc                                                                                                                                    >  1:946636/1‑150 (MQ=255)
                       aaaTCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaaa                                                                                                                                <  2:946636/150‑1 (MQ=255)
                          taccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCg                                                                                                                             >  1:2278967/3‑150 (MQ=255)
                            ccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                           >  1:1302862/1‑150 (MQ=255)
                                aCCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGa                                                                                                                       >  1:1756964/1‑150 (MQ=255)
                                 ccTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGAt                                                                                                                      >  2:311466/1‑150 (MQ=255)
                                 ccTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGAt                                                                                                                      >  2:538464/1‑150 (MQ=255)
                                      aCTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCt                                                                                                                 >  1:219751/1‑150 (MQ=255)
                                           tCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGc                                                                                                            >  2:2136304/1‑150 (MQ=255)
                                           tCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGc                                                                                                            >  1:741848/1‑150 (MQ=255)
                                             gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaaa                                                                                                                                >  1:160351/1‑128 (MQ=255)
                                             gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCaa                                                                                                          >  2:9350/1‑150 (MQ=255)
                                             gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGACCATCAAAACCGAAGTGATCGGCTTCCGCaa                                                                                                          <  1:1531209/150‑1 (MQ=255)
                                             gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCAGCGGACACCCTATCCACTCCGTCGACCTGATCCAGCAAACCGAAGTGACCGGCTTCCGCaa                                                                                                          >  2:1270584/1‑150 (MQ=255)
                                                gaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGa                                                                                                       <  1:1093120/150‑1 (MQ=255)
                                                gaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGa                                                                                                       >  2:414160/1‑150 (MQ=255)
                                                 aaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATcc                                                                                                                                       >  2:2205572/1‑118 (MQ=255)
                                                      ccGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTg                                                                                                                                           >  2:1059406/1‑109 (MQ=255)
                                                      ccGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTg                                                                                                                                           <  1:1059406/109‑1 (MQ=255)
                                                      ccGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGACGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                                 >  1:255904/1‑150 (MQ=255)
                                                       cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACGGAAGTGATCGGCTTCCGCAAGGAAAACGgc                                                                                                >  1:912790/1‑150 (MQ=255)
                                                       cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGa                                                                                                                       >  1:2408674/1‑127 (MQ=255)
                                                       cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGa                                                                                                                       <  2:2408674/127‑1 (MQ=255)
                                                          tACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAg                                                                                                         <  2:1633969/138‑1 (MQ=255)
                                                          tACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAg                                                                                                         >  1:1633969/1‑138 (MQ=255)
                                                            cccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                                           >  2:1733170/1‑150 (MQ=255)
                                                            cccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                                           >  1:1565467/1‑150 (MQ=255)
                                                             ccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTc                                                                                          >  2:1453569/1‑150 (MQ=255)
                                                             ccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTc                                                                                          >  2:357649/1‑150 (MQ=255)
                                                             ccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCGGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTc                                                                                          >  1:2104804/1‑150 (MQ=255)
                                                               gATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCAt                                                                                        >  2:2041946/1‑150 (MQ=255)
                                                                 tCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCg                                                                                      <  1:462383/150‑1 (MQ=255)
                                                                   cTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                    <  2:1976904/150‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAAc                                                                                                                               >  1:352024/1‑106 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAAc                                                                                                                               <  2:352024/106‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCg                                                                                      >  2:2171339/1‑147 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCg                                                                                      <  1:2171339/147‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtg                                                                                   <  1:1733170/150‑1 (MQ=255)
                                                                     cGGCGCCACGTTACAGCGCCGTGGTGGCTTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                  <  2:910812/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                  <  2:2410380/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                  <  2:161680/150‑1 (MQ=255)
                                                                           cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGa                                                                                                                                          <  2:134305/89‑1 (MQ=255)
                                                                           cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGa                                                                                                                                          >  1:134305/1‑89 (MQ=255)
                                                                           cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                    >  1:1620027/1‑142 (MQ=255)
                                                                           cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                    <  2:1620027/142‑1 (MQ=255)
                                                                            aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                      >  1:78385/1‑76 (MQ=255)
                                                                            aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                      <  2:78385/76‑1 (MQ=255)
                                                                                 aCAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACcaacaa                                                                      <  1:491186/150‑1 (MQ=255)
                                                                                     cgcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                                  <  1:311466/150‑1 (MQ=255)
                                                                                      gcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCt                                                                 <  1:275433/150‑1 (MQ=255)
                                                                                          tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCAt                                                             <  1:414160/150‑1 (MQ=255)
                                                                                             tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGtgtg                                                                                 >  1:61956/1‑127 (MQ=255)
                                                                                                 gTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCGCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGTTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGcc                                                      >  1:771156/1‑150 (MQ=255)
                                                                                                 gTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGcc                                                      >  2:2059125/1‑150 (MQ=255)
                                                                                                  ttGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCa                                                     <  2:1670246/150‑1 (MQ=255)
                                                                                                       cGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAACCGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAgcgc                                                <  2:219751/150‑1 (MQ=255)
                                                                                                        gTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCg                                                                                                                             >  2:821441/1‑72 (MQ=255)
                                                                                                        gTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCg                                                                                                                             <  1:821441/72‑1 (MQ=255)
                                                                                                        gTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCg                                                                                                                    <  2:1055056/81‑1 (MQ=255)
                                                                                                        gTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCg                                                                                                                    >  1:1055056/1‑81 (MQ=255)
                                                                                                        gTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGcgcg                                               >  1:2328507/1‑150 (MQ=255)
                                                                                                           aCGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCg                                            >  2:1038471/1‑150 (MQ=255)
                                                                                                             gATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGc                                          >  1:2383145/1‑150 (MQ=255)
                                                                                                                gCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTg                                       <  2:1302862/150‑1 (MQ=255)
                                                                                                                        cTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                    <  2:1517001/97‑1 (MQ=255)
                                                                                                                        cTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTAGGCGTGGTCAccgc                               >  1:2176963/1‑150 (MQ=255)
                                                                                                                         tGGGGCTTCGCCCGCGCTGCCGTCGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgcc                              >  2:760623/1‑150 (MQ=255)
                                                                                                                          ggggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgccg                             >  2:1748357/1‑150 (MQ=255)
                                                                                                                              cTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAgg                                                                   >  2:484201/1‑108 (MQ=255)
                                                                                                                              cTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAgg                                                                   <  1:484201/108‑1 (MQ=255)
                                                                                                                               ttCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAg                                                                    <  2:1969429/106‑1 (MQ=255)
                                                                                                                               ttCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAg                                                                    >  1:1969429/1‑106 (MQ=255)
                                                                                                                                   cccGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACaaa                                                                     <  1:1160033/101‑1 (MQ=255)
                                                                                                                                   cccGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACaaa                                                                     >  2:1160033/1‑101 (MQ=255)
                                                                                                                                     cgcgCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                             >  2:616162/1‑91 (MQ=255)
                                                                                                                                     cgcgCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                             <  1:616162/91‑1 (MQ=255)
                                                                                                                                         aggtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGcaca              >  2:376629/5‑150 (MQ=255)
                                                                                                                                          tGCCGACGCCCTGGGCGTCGACCTGTTCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACAt             <  2:1756964/150‑1 (MQ=255)
                                                                                                                                          tGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACAt             >  2:153490/1‑150 (MQ=255)
                                                                                                                                          tGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTAGACACAAGCAAAGGCTTCATCTGCGCCAAGCGCGTCGGCGTGGTCACCGCCCGTAGCTCCGGGCACAt             >  1:1490045/1‑150 (MQ=255)
                                                                                                                                           gtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                             >  1:958585/3‑85 (MQ=255)
                                                                                                                                           gtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                             <  2:958585/83‑1 (MQ=255)
                                                                                                                                            ccGACGCCCTGGGCGTCGACCTGATCCAGCAAACCTAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATgg           <  1:538464/150‑1 (MQ=255)
                                                                                                                                              gACGCCCTGGGCGTCGACCTGATCCAGNAAACCGAAGTGANCGGCTTNCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGcc         >  1:1984423/1‑150 (MQ=255)
                                                                                                                                                cGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCaa       <  1:9350/150‑1 (MQ=255)
                                                                                                                                                cGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCaa       <  2:2104804/150‑1 (MQ=255)
                                                                                                                                                   ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAGACCAACAAAGGCTTAATCGCCGCAAAGCGCGGCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCt    >  1:542693/1‑150 (MQ=255)
                                                                                                                                                   ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgccg                             >  2:698488/1‑125 (MQ=255)
                                                                                                                                                   ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCt    >  1:1816931/1‑150 (MQ=255)
                                                                                                                                                     tGGGCGACGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTgg  <  2:255904/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGG  >  NC_002947/388576‑388874

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: