Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,570,754 A→T 8.0% C578S (TGC→AGC)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,570,7540AT 8.0% 1055.9 / 55.7 326C578S (TGC→AGC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (141/159);  new base T (11/15);  total (152/174)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 335 total aligned reads displayed.

AAAGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5570604‑5570903
                                                                                                                                                      |                                                                                                                                                     
aaaGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATACGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCa                                                                                                                                                       >  1:6462227/1‑150 (MQ=255)
 aaGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCa                                                                                                                                                       >  2:134932/1‑150 (MQ=255)
    tttGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGc                                                                                                                                                    <  1:4696907/150‑1 (MQ=255)
      tGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCa                                                                                                                                                  <  2:413327/150‑1 (MQ=255)
        tCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                <  1:6643961/150‑1 (MQ=255)
        tCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                >  2:5203550/1‑150 (MQ=255)
        tCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                >  2:6281554/1‑150 (MQ=255)
          cAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGt                                                                                                                                              <  2:4473749/150‑1 (MQ=255)
               cGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCgt                                                                                                                                         <  1:5726675/150‑2 (MQ=255)
               cGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGa                                                                                                                                         <  2:4838270/150‑1 (MQ=255)
               cGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCATGGCTTGCCACGGTGCCgt                                                                                                                                         <  2:5442497/150‑2 (MQ=255)
                 tCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagc                                                                                                                                       >  1:3991451/1‑150 (MQ=255)
                 tCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagc                                                                                                                                       <  2:6782625/150‑1 (MQ=255)
                   cAGCCCCATCTGGGTCCCGGTAGGCCGACCATGGTTGCACTGGCCGCTGCGCGCGGGGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCTCCGTGGCAAGCCATGGTGCCGagcag                                                                                                                                     >  2:6196141/1‑150 (MQ=255)
                   cAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCGtgcgc                                                                                                                                     <  1:70242/150‑3 (MQ=255)
                      ccccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTc                                                                                                                                  <  1:2329415/150‑1 (MQ=255)
                      ccccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTc                                                                                                                                  <  1:5795684/150‑1 (MQ=255)
                        ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                >  2:3569517/1‑150 (MQ=255)
                        ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTGTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCag                                                                                                                                >  1:1197934/1‑149 (MQ=255)
                          aTCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTg                                                                                                                              >  1:6231611/1‑150 (MQ=255)
                           tCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGa                                                                                                                             >  1:6402984/1‑150 (MQ=255)
                           tCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGa                                                                                                                             >  2:6322653/1‑150 (MQ=255)
                              gggTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGt                                                                                                                          <  2:3038198/150‑1 (MQ=255)
                              gggTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGt                                                                                                                          >  2:1568011/1‑150 (MQ=255)
                                gTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                               <  2:513096/143‑1 (MQ=255)
                                gTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgc                                                                                                                        <  1:6348681/150‑1 (MQ=255)
                                  ccAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgcgc                                                                                                                      >  1:5936869/1‑150 (MQ=255)
                                  ccAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgcgc                                                                                                                      <  1:2211407/150‑1 (MQ=255)
                                     ggTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTg                                                                                                                   >  2:2850075/1‑150 (MQ=255)
                                        aGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAAt                                                                                                                <  1:2797930/150‑1 (MQ=255)
                                         ggACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATg                                                                                                               <  1:5508621/150‑1 (MQ=255)
                                          gtcgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGACAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGc                                                                                                              <  2:6499813/148‑1 (MQ=255)
                                          gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGc                                                                                                              >  2:3512700/1‑150 (MQ=255)
                                           acgacCATGGGTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGGTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCg                                                                                                             <  1:1969448/150‑1 (MQ=255)
                                            cgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGa                                                                                                            <  2:2646021/150‑1 (MQ=255)
                                              acCATGGTTGCACTGGCCGCTGCTCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCTAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATc                                                                                                          <  1:6831670/150‑1 (MQ=255)
                                               ccATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCa                                                                                                         >  1:3409567/1‑150 (MQ=255)
                                               ccATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCa                                                                                                         >  1:3298024/1‑150 (MQ=255)
                                                 aTGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACt                                                                                                       <  1:387436/150‑1 (MQ=255)
                                                 aTGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGCCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACt                                                                                                       >  2:3765661/1‑150 (MQ=255)
                                                     ttGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTg                                                                                                   <  1:4521058/150‑1 (MQ=255)
                                                      tGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGc                                                                                                  >  2:1863815/1‑150 (MQ=255)
                                                       gCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                 <  1:1538693/150‑1 (MQ=255)
                                                       gCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                 <  1:6789382/150‑1 (MQ=255)
                                                       gCACTGGCCGCTGCGCTCGGTGTTTGCCATGTCGCGCAGCAGGGCGTTCATCGCGGGAAGTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                 <  1:6694214/150‑1 (MQ=255)
                                                            ggCCGCTGCGCTCGGTGTTTTCCCTGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACt                                                                                            <  2:4645243/150‑1 (MQ=255)
                                                            ggCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACt                                                                                            <  2:2932930/150‑1 (MQ=255)
                                                              ccGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTcc                                                                                          <  2:2011740/150‑1 (MQ=255)
                                                                gCTGCGCTCGGTGTTTTCCATGTCGCACAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCAt                                                                                        <  1:2537931/150‑1 (MQ=255)
                                                                  tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCGATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCAtta                                                                                      <  1:5154459/150‑3 (MQ=255)
                                                                  tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                      <  1:3685020/150‑1 (MQ=255)
                                                                  tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                      >  1:429891/1‑150 (MQ=255)
                                                                  tGCGCTCGGGGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGAACTCCATaa                                                                                      >  2:4682287/1‑150 (MQ=255)
                                                                      cTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTc                                                                                  <  1:318510/150‑1 (MQ=255)
                                                                      cTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTc                                                                                  <  1:973631/150‑1 (MQ=255)
                                                                       tCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCg                                                                                 >  1:777232/1‑150 (MQ=255)
                                                                       tCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCg                                                                                 <  1:2589159/150‑1 (MQ=255)
                                                                       tCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCg                                                                                 >  1:4208860/1‑150 (MQ=255)
                                                                        cGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCgg                                                                                <  1:844563/150‑1 (MQ=255)
                                                                        cGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCgg                                                                                <  1:3919178/150‑1 (MQ=255)
                                                                         ggTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                               >  2:4521710/1‑150 (MQ=255)
                                                                         ggTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                               >  2:5611467/1‑150 (MQ=255)
                                                                           tgtTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCaa                                                                             >  2:6552281/1‑150 (MQ=255)
                                                                            gtTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAg                                                                            <  2:4247867/150‑1 (MQ=255)
                                                                             ttttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGc                                                                           <  1:2907014/150‑1 (MQ=255)
                                                                              tttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGTCGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGTATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAgac                                                                          <  2:6856219/150‑3 (MQ=255)
                                                                               ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCAc                                                                         >  2:2558207/1‑150 (MQ=255)
                                                                                 ccATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGt                                                                       >  2:1802936/1‑150 (MQ=255)
                                                                                  cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTc                                                                      >  2:1223587/1‑150 (MQ=255)
                                                                                     gTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTg                                                                   <  1:1863815/150‑1 (MQ=255)
                                                                                        gcgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCAc                                                                >  2:2346841/1‑150 (MQ=255)
                                                                                               agGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCgg                                                         >  2:757309/1‑150 (MQ=255)
                                                                                                 ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                       <  1:6574592/150‑1 (MQ=255)
                                                                                                 ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                       <  2:3709891/150‑1 (MQ=255)
                                                                                                 ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                       <  2:2198914/150‑1 (MQ=255)
                                                                                                  gCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTg                                                      <  2:831890/150‑1 (MQ=255)
                                                                                                    gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGc                                                    <  1:2383297/150‑1 (MQ=255)
                                                                                                       cATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTc                                                 <  2:3930591/150‑1 (MQ=255)
                                                                                                          ctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGc                                              <  1:4974975/150‑1 (MQ=255)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGc                                                                                                                                            >  1:5310091/1‑55 (MQ=18)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGc                                                                                                                                            <  2:5310091/55‑1 (MQ=18)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  1:4948016/1‑150 (MQ=255)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  1:4393533/1‑150 (MQ=255)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  1:3612670/1‑150 (MQ=255)
                                                                                                           tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  2:4100249/1‑150 (MQ=255)
                                                                                                             gggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTg                                           >  1:194459/1‑150 (MQ=255)
                                                                                                             gggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAGTGAGGTGCGCCTGACTGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTg                                           >  2:5833431/1‑150 (MQ=255)
                                                                                                              ggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGc                                          >  2:4629190/1‑150 (MQ=255)
                                                                                                               gAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCt                                         <  1:2832374/150‑1 (MQ=255)
                                                                                                                aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                        <  2:4324096/150‑1 (MQ=255)
                                                                                                                aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                        <  1:5108525/150‑1 (MQ=255)
                                                                                                                aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                        >  2:5067893/1‑150 (MQ=255)
                                                                                                                                   gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTgg                     <  2:6461677/150‑1 (MQ=255)
                                                                                                                                   gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCCTCAGCCTGCTTGAGCAAGGCCGGGATCTgg                     <  2:6611401/150‑1 (MQ=255)
                                                                                                                                      tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                  >  1:5265007/5‑150 (MQ=255)
                                                                                                                                      cgcACGGCACCGTGGCAAGCCATGGTGCCGAGTAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                  >  1:1035947/1‑150 (MQ=255)
                                                                                                                                       gcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGACCTGGCGaa                 <  2:3376282/150‑1 (MQ=255)
                                                                                                                                           ggCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATTCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCATGCACCAGGCGGTTGGCTTCAGACTGCTTGAGAAAGGCCGGGATCTGGCGAATCGc             >  1:2992008/1‑150 (MQ=255)
                                                                                                                                            gCACCATGGCTTGCCACGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGcc            >  2:3205753/1‑150 (MQ=255)
                                                                                                                                             cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa           <  1:4274083/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAAGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5570604‑5570903

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