Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,213,954 A→C 5.2% I1002S (ATC→AGC) ‡ PP_1061 ← ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,213,9540AC 5.2% 248.3 / 3.8 77I1002S (ATC→AGC) ‡PP_1061ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base C (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGGTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACAGTGCCAG  >  NC_002947/1213807‑1214096
                                                                                                                                                      |                                                                                                                                              
cGGTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGAt                                                                                                                                               >  5:162917/1‑149 (MQ=255)
 ggTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAATTGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATc                                                                                                                                              <  3:252699/149‑1 (MQ=255)
  gTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCt                                                                                                                                             >  5:14481/1‑149 (MQ=255)
  gTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCt                                                                                                                                             >  6:40781/1‑149 (MQ=255)
     gATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGTAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTgc                                                                                                                                           >  1:318725/1‑149 (MQ=255)
             ccAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAtggcc                                                                                                                                  >  2:179398/1‑147 (MQ=255)
                      ggCGTTCGCACAATTGCAGGTCATTTGCACCTGGGCTGAAATGCCCGCCCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGGTACCGGTTTGCCAAAAGGCGCACCGATCTGCGCCAGGGACTGGGGGCCa                                                                                                                         >  2:315232/1‑149 (MQ=255)
                             gCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGc                                                                                                                  <  6:162917/149‑1 (MQ=255)
                               aCCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCt                                                                                                                 >  3:237694/1‑148 (MQ=255)
                                  aTTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTCCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAgg                                                                                                             <  2:355580/149‑1 (MQ=255)
                                      cAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGAGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcgg                                                                                                          <  1:297575/149‑1 (MQ=255)
                                         gTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcgg                                                                                                          >  7:38412/1‑145 (MQ=255)
                                         gTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcgg                                                                                                          >  8:178956/1‑145 (MQ=255)
                                         gTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcgg                                                                                                          <  7:178956/145‑1 (MQ=255)
                                         gTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcgg                                                                                                          <  8:38412/145‑1 (MQ=255)
                                             gTTGCACCTGGGCTGAAATGCCCGCGCAGCACATGAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGgcgcg                                                                                                     >  1:55179/1‑147 (MQ=255)
                                             gTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCgcag                                                                                                  <  2:318725/149‑1 (MQ=255)
                                             gTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGATCGCTGCGCCa                                                                                                                                       <  5:383877/112‑1 (MQ=255)
                                             gTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGATCGCTGCGCCa                                                                                                                                       >  6:383877/1‑112 (MQ=255)
                                                 cACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCt                                                                                                                                             <  6:168476/102‑1 (MQ=255)
                                                 cACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCt                                                                                                                                             >  5:168476/1‑102 (MQ=255)
                                                  aCCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGtt                                                                                             >  5:54109/1‑149 (MQ=255)
                                                  aCCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGtt                                                                                             >  2:208359/1‑149 (MQ=255)
                                                   ccTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTc                                                                                            >  2:117708/1‑149 (MQ=255)
                                                        gCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTgggg                                                                                                                             >  5:38267/1‑111 (MQ=255)
                                                        gCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTgggg                                                                                                                             <  6:38267/111‑1 (MQ=255)
                                                        gCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGATCGCTGCGCCa                                                                                                                                       <  7:339536/101‑1 (MQ=255)
                                                        gCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGATCGCTGCGCCa                                                                                                                                       >  8:339536/1‑101 (MQ=255)
                                                         cTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAg                                                                                      >  7:44213/1‑149 (MQ=255)
                                                          tGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGt                                                                                     >  3:267778/1‑149 (MQ=255)
                                                               aTGCCCGCGCAGCTCATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGt                                                                                 <  4:237694/149‑1 (MQ=255)
                                                                    cgcgCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCg                                                                            <  4:267778/149‑1 (MQ=255)
                                                                          gcaCATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGTTATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTGAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCg                                                                      <  1:315232/149‑1 (MQ=255)
                                                                            acaTAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGAt                                                                    <  8:44213/149‑1 (MQ=255)
                                                                              aTAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATgg                                                                  >  3:225245/1‑149 (MQ=255)
                                                                               tAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATggg                                                                 >  2:19362/1‑149 (MQ=255)
                                                                                 aGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTg                                                                                                                                >  4:370974/1‑85 (MQ=255)
                                                                                 aGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGc                                                                >  2:176895/1‑149 (MQ=255)
                                                                                   ccTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcggtg                                                                                                        >  5:113871/1‑105 (MQ=255)
                                                                                   ccTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAggcggtg                                                                                                        <  6:113871/107‑3 (MQ=255)
                                                                                   ccTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTg                                                              <  4:257990/149‑1 (MQ=255)
                                                                                   ccTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTg                                                              <  5:40781/149‑1 (MQ=255)
                                                                                    cTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTgg                                                             >  2:348075/1‑149 (MQ=255)
                                                                                    cTTCGGCTTACAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGATCAGTGCGTGGTCGGGCTCGATGGGCTgg                                                             >  4:406646/1‑149 (MQ=255)
                                                                                         gCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGaa                                                        <  2:55179/149‑1 (MQ=255)
                                                                                               aGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGc                                                  >  2:150873/1‑149 (MQ=255)
                                                                                                      ccAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGt                                           <  6:14481/149‑1 (MQ=255)
                                                                                                        agagCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTc                                         <  7:333306/149‑1 (MQ=255)
                                                                                                         gagCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGt                                                                              >  8:359396/1‑111 (MQ=255)
                                                                                                         gagCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGt                                                                              <  7:359396/111‑1 (MQ=255)
                                                                                                           gCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGa                                       <  1:348075/148‑1 (MQ=255)
                                                                                                                cTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCt                                  <  1:117708/148‑1 (MQ=255)
                                                                                                                    atatCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTc                             <  1:176895/149‑1 (MQ=255)
                                                                                                                    atatCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTc                             >  7:340466/1‑149 (MQ=255)
                                                                                                                       tCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCt                          <  6:54109/149‑1 (MQ=255)
                                                                                                                       tCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCt                          <  1:208359/149‑1 (MQ=255)
                                                                                                                        cGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGTCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGtt                                          <  7:267296/132‑1 (MQ=255)
                                                                                                                        cGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGTCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGtt                                          >  8:267296/1‑132 (MQ=255)
                                                                                                                         ggTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCtgt                        <  4:35184/149‑1 (MQ=255)
                                                                                                                           ttACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTggg                                                                                                                              >  1:246749/1‑45 (MQ=255)
                                                                                                                           ttACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTggg                                                                                                                              <  2:246749/45‑1 (MQ=255)
                                                                                                                                 gTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCtgtg                >  7:317420/1‑149 (MQ=255)
                                                                                                                                 gTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCtgtg                >  5:1876/1‑149 (MQ=255)
                                                                                                                                      cccAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGt                                                            <  7:265963/100‑1 (MQ=255)
                                                                                                                                      cccAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGt                                                            >  8:265963/1‑100 (MQ=255)
                                                                                                                                       ccAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTa           >  6:219687/1‑148 (MQ=255)
                                                                                                                                       ccAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTa           >  1:256525/1‑148 (MQ=255)
                                                                                                                                        cAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTc                                                                       <  5:267512/87‑1 (MQ=255)
                                                                                                                                        cAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTc                                                                       >  6:267512/1‑87 (MQ=255)
                                                                                                                                        cAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTa           >  8:92724/1‑147 (MQ=255)
                                                                                                                                        cAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACa         <  7:92724/149‑1 (MQ=255)
                                                                                                                                        cAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACa         <  1:19362/149‑1 (MQ=255)
                                                                                                                                         aaTGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACAg        >  7:369485/1‑149 (MQ=255)
                                                                                                                                         aaTGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACAg        >  7:364055/1‑149 (MQ=255)
                                                                                                                                          ctGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGaa                                                        <  3:359691/99‑1 (MQ=255)
                                                                                                                                          ctGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGaa                                                        >  4:359691/2‑100 (MQ=255)
                                                                                                                                                gcAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACAGTGCcag  >  4:9695/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CGGTGGATGCGTGCCAGCAAGTGGCGTTCGCACCATTGCAGGTCAGTTGCACCTGGGCTGAAATGCCCGCGCAGCACATAGCCTTCGGCTTCCAGGCGGGCCAGAGCCTGCTCGATATCGGTTACCGGTTTGCCCAATGGCGCAGCGATCTGCGCCAGGGTCTGGGGGCCATGGCCGCTCAGGCGGGCGCGCAGCAGTTCGCTCAGTGCGTGGTCGGGCTCGATGGGCTGGTCGAAGCCAGCGAGCCGTTCGAGCGCTGGCTCCCTGTGCGCTGTGGGGTACAGTGCCAG  >  NC_002947/1213807‑1214096

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