Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,215 G→C 10.3% V63L (GTT→CTT)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,2150GC10.3% 228.7 / 14.9 78V63L (GTT→CTT) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base C (5/3);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

CTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTCGG  >  NC_002947/571067‑571363
                                                                                                                                                    |                                                                                                                                                     
cTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCAGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAgg                                                                                                                                                       >  5:184089/1‑149 (MQ=255)
   tgGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGt                                                                                                                                                      <  6:322122/147‑1 (MQ=255)
      ccAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGcc                                                                                                                                                  >  4:242205/1‑148 (MQ=255)
          gATCCCGAACTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaaca                                                                                                                                             <  6:213203/149‑3 (MQ=255)
            tCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaacagc                                                                                                                                           <  7:66183/149‑5 (MQ=255)
               cGAAGTCCGCTACCTGCCCAACGGTAACGCGGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTg                                                                                                                                        <  6:184089/149‑1 (MQ=255)
                gAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                       >  7:35574/1‑149 (MQ=255)
                gAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                       <  2:70434/149‑1 (MQ=255)
                gAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                       >  2:235395/1‑149 (MQ=255)
                      cGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAAt                                                                                                                                  <  2:370114/148‑1 (MQ=255)
                      cGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaacag                                                                                                                                            <  8:173458/138‑4 (MQ=255)
                      cGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaacag                                                                                                                                            >  7:173458/1‑135 (MQ=255)
                       gCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAc                                                                                                                                <  6:262655/149‑1 (MQ=255)
                           ccTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTgc                                                                                                                            <  7:103367/149‑1 (MQ=255)
                           ccTGCCCAACCCTACCGCCGTCACCAACCTGAGCCTTGCCACCAGCGATCAGTTGACTGACAAGCAGGCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGGTCGGCAAGGATGCCGAAATCGCTGGCGAATACCTgc                                                                                                                            <  4:229220/149‑1 (MQ=255)
                                 cAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAgg                                                                                                                      <  3:242205/149‑1 (MQ=255)
                                    cGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGtt                                                                                                                   <  2:112112/149‑1 (MQ=255)
                                      gTAACGCCGTGACCAACCTCATCCTGGCCAAAAGCGAGCAGTGAACTGACAAGCAGTCGCGCCAGAACCTCGAGCGTACCGAGTGGCACCGTGTGTCCCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAACGGATCg                                                                                                                 >  1:73373/1‑149 (MQ=255)
                                         aCGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAg                                                                                                              >  6:56242/1‑149 (MQ=255)
                                         aCGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAg                                                                                                              >  8:184920/1‑149 (MQ=255)
                                           gCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTg                                                                                                            <  8:69137/149‑1 (MQ=255)
                                            ccGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGc                                                                                                           >  3:312404/1‑149 (MQ=255)
                                             cGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGc                                                                                                           >  6:96899/1‑148 (MQ=255)
                                               tGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGt                                                                                                             <  2:66675/144‑1 (MQ=255)
                                               tGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGt                                                                                                             >  1:66675/1‑144 (MQ=255)
                                               tGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTAc                                                                                                       <  7:19158/149‑1 (MQ=255)
                                                gACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACa                                                                                                      >  2:289148/1‑149 (MQ=255)
                                                 aCCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACa                                                                                                      <  1:275385/148‑1 (MQ=255)
                                                     aCCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGaa                                                                                                 <  5:56242/149‑1 (MQ=255)
                                                     aCCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGaa                                                                                                 <  6:322894/149‑1 (MQ=255)
                                                     aCCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGaa                                                                                                 <  5:47942/149‑1 (MQ=255)
                                                         gAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGc                                                                                              >  1:207271/1‑148 (MQ=255)
                                                            ccTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGc                                                                                          >  2:184754/1‑149 (MQ=255)
                                                            ccTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGc                                                                                          <  1:184754/149‑1 (MQ=255)
                                                             cTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCt                                                                                         <  7:184920/149‑1 (MQ=255)
                                                                gCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCa                                                                                      <  5:96899/149‑1 (MQ=255)
                                                                   accaGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAc                                                                                   <  4:312404/149‑1 (MQ=255)
                                                                   accaGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAc                                                                                   >  6:97924/1‑149 (MQ=255)
                                                                     caGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAccc                                                                                 >  2:236222/1‑149 (MQ=255)
                                                                     caGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAccc                                                                                 >  8:23860/1‑149 (MQ=255)
                                                                            cAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTg                                                                          <  1:236222/149‑1 (MQ=255)
                                                                             aGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTTCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                         <  5:353943/149‑1 (MQ=255)
                                                                              gTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTggg                                                                        <  2:250142/149‑1 (MQ=255)
                                                                                 gACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGaa                                                                      >  2:75297/1‑148 (MQ=255)
                                                                                  aCTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGaaa                                                                     >  7:225578/1‑148 (MQ=255)
                                                                                  aCTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGaaa                                                                     >  7:231387/1‑148 (MQ=255)
                                                                                    tGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGaaa                                                                     >  7:120551/1‑146 (MQ=255)
                                                                                      aCAAGCAGTCGTGCCAGAAGGTCGAGCGTAGCGTGTGGCACCGTGTGTCGCTGTTCTGCAATGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTAGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCTAGTGGGAAAAAGAc                                                                <  6:95773/149‑1 (MQ=255)
                                                                                      aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGAc                                                                >  8:182106/1‑149 (MQ=255)
                                                                                       cAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGCCCCGCGAGTGGGAAAAAGACg                                                               >  6:170864/1‑149 (MQ=255)
                                                                                       cAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACg                                                               <  1:75297/149‑1 (MQ=255)
                                                                                           cAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCa                                                            >  4:205011/1‑148 (MQ=255)
                                                                                            aGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCTAGTGGGAAAAAGACGGCATc                                                          <  2:298352/149‑1 (MQ=255)
                                                                                            aGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCAt                                                           >  3:352764/1‑148 (MQ=255)
                                                                                            aGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATc                                                          <  1:289148/149‑1 (MQ=255)
                                                                                            aGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATc                                                          <  8:45168/149‑1 (MQ=255)
                                                                                                tggCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATc                                                          >  6:266790/2‑145 (MQ=255)
                                                                                                tggCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATc                                                          <  5:266790/144‑1 (MQ=255)
                                                                                                  gCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAgcgc                                                    >  6:207062/1‑149 (MQ=255)
                                                                                                    cAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTa                                                  <  3:205011/149‑1 (MQ=255)
                                                                                                    cAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTa                                                  >  6:93762/1‑149 (MQ=255)
                                                                                                     aGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTac                                                 <  7:97444/149‑1 (MQ=255)
                                                                                                       aaGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTacac                                               <  1:235395/149‑1 (MQ=255)
                                                                                                              aGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACaccacc                                           >  7:217812/1‑146 (MQ=255)
                                                                                                                    ccGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGatcat                                   <  7:182106/148‑1 (MQ=255)
                                                                                                                     cGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCAtcg                                 >  2:63479/1‑149 (MQ=255)
                                                                                                                      gAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCAtcgt                                >  5:42303/1‑149 (MQ=255)
                                                                                                                           gCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGata                           <  5:170864/149‑1 (MQ=255)
                                                                                                                                 tgtgtCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACg                     <  5:93762/149‑1 (MQ=255)
                                                                                                                                    cagcgatttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCa                  >  8:79616/8‑149 (MQ=255)
                                                                                                                                    cagcgatttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCa                  >  4:302174/8‑149 (MQ=255)
                                                                                                                                       gCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCGGACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCAc                 >  1:109067/1‑148 (MQ=255)
                                                                                                                                       gCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCAcc                <  2:207271/148‑1 (MQ=255)
                                                                                                                                         tGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTTGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATg             >  1:124645/1‑149 (MQ=255)
                                                                                                                                         atttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATg             >  1:344289/3‑149 (MQ=255)
                                                                                                                                         atttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATg             >  1:175527/3‑149 (MQ=255)
                                                                                                                                               gCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGCTCATCGTCGAGATCAGCGGCACCAggcggcgg       >  6:261612/1‑147 (MQ=255)
                                                                                                                                               gCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAgctg       <  8:217812/149‑1 (MQ=255)
                                                                                                                                                 aaGGGTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATAGAACGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACAACCGAGATCATCTTCGATATCAACGGCACCATGCAgctgct     >  5:259479/1‑149 (MQ=255)
                                                                                                                                                    gTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTcgg  <  8:35574/149‑1 (MQ=255)
                                                                                                                                                    gTTGACGAAATCGCTGGCGAATTCCTGCGCATGGGTTCGCTGTTCTACATCGAAGGCAAGCTGCAGACACGCGTGTGGGAAATAGACGGCATGAAGCGCTTGACCACCGAGATCATCGTCGATATCAACGGCACCATGCTGCTGCTcgg  <  2:73373/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTCGG  >  NC_002947/571067‑571363

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