Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,713,554 G→T 5.8% P152P (CCG→CCT)  tsaB → endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,713,5540GT 5.8% 279.4 / 4.9 86P152P (CCG→CCT) tsaBendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (30/51);  new base T (2/3);  total (32/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GCATTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGA  >  NC_002947/1713405‑1713703
                                                                                                                                                     |                                                                                                                                                     
gCATTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCg                                                                                                                                                       <  2:846497/150‑1 (MQ=255)
gCATTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCg                                                                                                                                                       >  1:2067869/1‑150 (MQ=255)
 cATTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCgg                                                                                                                                                      >  1:1720025/1‑150 (MQ=255)
  aTTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCCGGTGGCGCTGCCggc                                                                                                                                                     >  1:88217/1‑150 (MQ=255)
   ttGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAGAGTGTAAGGACACCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCggca                                                                                                                                                    >  2:1268953/1‑150 (MQ=255)
     gcgcGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCggcagg                                                                                                                                                  <  1:2155352/150‑1 (MQ=255)
        tgAGCATGGCGGGCAGCCGGTGGCCGTTGCTATCAATGCGCGCATGGATGAGGTTTACTGGGGCTGTTACCAGGCCATAGCGGGCGAAATGCGCCTGTTGGGCTGTGAAGCGGTGTTGCCGCCCGAGCAGTTGGCGCTGCCGgaaggcag                                                                                                                                               <  2:450689/149‑1 (MQ=255)
         gAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGgcagc                                                                                                                                              >  2:1409896/1‑150 (MQ=255)
         gAGCAGGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCTGCCGgcagc                                                                                                                                              <  1:93414/150‑1 (MQ=255)
              tGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATgg                                                                                                                                         <  1:1238562/150‑1 (MQ=255)
               ggCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCTGCCGgcagcgccacc                                                                                                                                        <  1:676852/150‑7 (MQ=255)
               ggCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGc                                                                                                                                        <  2:1030033/150‑1 (MQ=255)
                 cGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGa                                                                                                                                      <  2:460979/150‑1 (MQ=255)
                  gTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCACCAATgcccag                                                                                                                                     <  2:32665/150‑6 (MQ=255)
                    gcagcaGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTg                                                                                                                                   <  1:754220/150‑1 (MQ=255)
                    gcagcaGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTg                                                                                                                                   <  2:2067869/150‑1 (MQ=255)
                      agcagGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGt                                                                                                                                 >  1:1428001/1‑150 (MQ=255)
                      agcagGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGt                                                                                                                                 <  2:1446921/150‑1 (MQ=255)
                         gggTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCg                                                                                                                              <  2:1800856/149‑1 (MQ=255)
                           ctgaccGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGGCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGa                                                                                                                                      >  2:1399665/5‑140 (MQ=255)
                           ctgaccGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGGCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGa                                                                                                                                      <  1:1399665/136‑1 (MQ=255)
                             ggCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGc                                                                                                                          <  2:246738/150‑1 (MQ=255)
                              gCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGcc                                                                                                                         <  1:800538/150‑1 (MQ=255)
                              gCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGCGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGcc                                                                                                                         <  1:2294319/150‑1 (MQ=255)
                                 gCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGt                                                                                                                      >  2:1162523/1‑150 (MQ=255)
                                   tGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTAc                                                                                                                    <  2:1285035/150‑1 (MQ=255)
                                   tGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTAc                                                                                                                    <  1:1409896/150‑1 (MQ=255)
                                     ccATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCg                                                                                                                  >  2:1767125/1‑150 (MQ=255)
                                      cATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCgg                                                                                                                 <  2:1428001/150‑1 (MQ=255)
                                       aTCGATGCGCGCGTGGATAAAGACTACTGGGGATGCGACCAGGCCACAGCGGACGAAATACGACTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGAGGGTCGGTGCCGGTACCggg                                                                                                                >  2:70022/1‑150 (MQ=255)
                                       aTCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCggg                                                                                                                >  2:2128536/1‑150 (MQ=255)
                                         cGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTg                                                                                                              <  1:21286/150‑1 (MQ=255)
                                          gATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTgg                                                                                                             <  1:1165223/150‑1 (MQ=255)
                                          gATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTgg                                                                                                             <  1:2321821/150‑1 (MQ=255)
                                                  cATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGc                                                                                                     <  2:1934866/150‑1 (MQ=255)
                                                     ggATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGAAATGGCGAGTGGTTCGGTGCCGGTTCCTGGTGGGGCTATGCCgc                                                                                                  <  2:594003/150‑2 (MQ=255)
                                                       aTGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGGTCGGTGCCGGTACCGGGTGGGGCTATGTCGCTc                                                                                                >  2:509851/1‑150 (MQ=255)
                                                         gAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGc                                                                                              >  1:769429/1‑150 (MQ=255)
                                                             tttACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGTGCTGCCGGCAGGCAGAAATGGCGAGTGGTTCGGTGCCGGTACCGGGGGGGGCTATGCCGATCGCCTgg                                                                                          >  2:1387531/1‑150 (MQ=255)
                                                               tACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCg                                                                                        <  2:769429/150‑1 (MQ=255)
                                                                 cTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGt                                                                                      >  1:512160/1‑150 (MQ=255)
                                                                  tGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCTGCCGgcagcgcca                                                                                                                                          <  2:2018015/97‑5 (MQ=255)
                                                                  tGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCTGCCGgcagcgcca                                                                                                                                          >  1:2018015/1‑93 (MQ=255)
                                                                  tGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTa                                                                                     <  2:2148563/150‑1 (MQ=255)
                                                                    gggCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACa                                                                                   <  1:2128536/150‑1 (MQ=255)
                                                                    gggCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGCGGCACGTCCGGCAGGCAGCAATGGCGAGAGTTGCGGATCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACa                                                                                   >  1:1740559/1‑150 (MQ=255)
                                                                     ggCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAg                                                                                  <  2:1937684/150‑1 (MQ=255)
                                                                      gctgctACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTCCAgg                                                                                 >  1:228936/1‑150 (MQ=255)
                                                                      gctgctACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAgg                                                                                 >  2:1330002/1‑150 (MQ=255)
                                                                            aCCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTa                                                                                                        >  1:812425/1‑121 (MQ=255)
                                                                            aCCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTa                                                                                                        <  2:812425/121‑1 (MQ=255)
                                                                             ccAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCgg                                                                                                                 <  1:1249970/111‑1 (MQ=255)
                                                                             ccAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCgg                                                                                                                 >  2:1249970/1‑111 (MQ=255)
                                                                              cAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAgcggcg                                                                         >  2:548481/1‑150 (MQ=255)
                                                                               aGGNCACAGCGGGCGAAATCAGCCGTGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGNANCGGGTGGGGCTATGCNGATCGCCTGGCGGTACAGGTAGCGGCGa                                                                        <  1:1330002/150‑1 (MQ=255)
                                                                               aGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGa                                                                        <  2:2005346/150‑1 (MQ=255)
                                                                               aGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTAAGCCGAGCGCCTGGCGGTACAGGTAGCGGCGa                                                                        <  1:1268953/150‑1 (MQ=255)
                                                                                ggCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTCACGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAg                                                                       <  1:305557/150‑1 (MQ=255)
                                                                                   cacaGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCAGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCaa                                                                    <  1:1266564/150‑1 (MQ=255)
                                                                                          ggCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCt                                                             <  1:182064/150‑1 (MQ=255)
                                                                                             gAAATGCACCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGcc                                                          >  1:1100603/1‑150 (MQ=255)
                                                                                              aaaTGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGCTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCt                                                         <  2:1740559/150‑1 (MQ=255)
                                                                                                     ccTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCaa                                                  <  2:2218197/150‑1 (MQ=255)
                                                                                                      cTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAAt                                                 >  1:589613/1‑150 (MQ=255)
                                                                                                       tGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATg                                                <  2:2306797/150‑1 (MQ=255)
                                                                                                           cGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGctct                                            <  1:724311/150‑1 (MQ=255)
                                                                                                            ggCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTg                                           >  1:1192969/1‑150 (MQ=255)
                                                                                                             gCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTgg                                          >  2:2106465/1‑150 (MQ=255)
                                                                                                              ccGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGa                                         <  1:1208056/150‑1 (MQ=255)
                                                                                                              ccGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGa                                         <  1:468191/150‑1 (MQ=255)
                                                                                                               cGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGAc                                        >  2:1771936/1‑150 (MQ=255)
                                                                                                                     gCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATtctc                                  >  2:770141/1‑150 (MQ=255)
                                                                                                                      ccGTGTTGCCGCCCGAGCAGGTGGCGCTGCCggca                                                                                                                                                    <  2:22973/35‑1 (MQ=255)
                                                                                                                      ccGTGTTGCCGCCCGAGCAGGTGGCGCTGCCggca                                                                                                                                                    >  1:22973/1‑35 (MQ=255)
                                                                                                                        gtgtTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGc                               <  1:1162523/150‑1 (MQ=255)
                                                                                                                               cgccCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTAc                                                                                    >  2:1609971/1‑90 (MQ=255)
                                                                                                                               cgccCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTAc                                                                                    <  1:1609971/90‑1 (MQ=255)
                                                                                                                                gccCGAGCAGGTGGCGCTGCCGGCAGGCAGGAATGGCGGGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCttggccc                       <  2:88217/150‑7 (MQ=255)
                                                                                                                                  ccGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCt                     >  1:72470/1‑150 (MQ=255)
                                                                                                                                   cGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCtt                    >  1:1006288/1‑150 (MQ=255)
                                                                                                                                    gAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCttt                   >  1:710427/1‑150 (MQ=255)
                                                                                                                                      gCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGc                 <  1:204039/150‑1 (MQ=255)
                                                                                                                                      gCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGc                 <  1:1767125/150‑1 (MQ=255)
                                                                                                                                        aGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCt               <  2:512160/150‑1 (MQ=255)
                                                                                                                                          gTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTgg                                                                                          >  2:954548/1‑73 (MQ=255)
                                                                                                                                          gTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTgg                                                                                          <  1:954548/73‑1 (MQ=255)
                                                                                                                                          gTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCTGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTgg             <  1:770141/150‑1 (MQ=255)
                                                                                                                                          gTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTgg             <  1:2106465/150‑1 (MQ=255)
                                                                                                                                            attgcTGCCTGCCGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGc           >  1:2180551/4‑150 (MQ=255)
                                                                                                                                                cTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACa                                       >  1:1763908/1‑118 (MQ=255)
                                                                                                                                                cTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACa                                       <  2:1763908/118‑1 (MQ=255)
                                                                                                                                                  gCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTgg     <  2:72470/150‑1 (MQ=255)
                                                                                                                                                   ccGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGc    >  1:1286010/1‑150 (MQ=255)
                                                                                                                                                   ccGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGc    <  1:854204/150‑1 (MQ=255)
                                                                                                                                                     ggcaggcaGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTTGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGATTATTTCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTcgcgc  <  2:369796/150‑2 (MQ=255)
                                                                                                                                                     ggcaggcaGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGa  <  2:1100603/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATTGCGCGAGCATGGCGTGCAGCAGGTGGCCGCTGCCATCGATGCGCGCATGGATGAGGTTTACTGGGGCTGCTACCAGGCCACAGCGGGCGAAATGCGCCTGGTCGGCCGTGAAGCCGTGTTGCCGCCCGAGCAGGTGGCGCTGCCGGCAGGCAGCAATGGCGAGTGGTTCGGTGCCGGTACCGGGTGGGGCTATGCCGATCGCCTGGCGGTACAGGTAGCGGCGAGCAACCCGGCTGCCTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGA  >  NC_002947/1713405‑1713703

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: ‑