Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,633,811 C→A 7.3% A434S (GCA→TCA) ‡ uvrC ← UvrABC system protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,633,8110CA 7.3% 210.4 / 3.5 67A434S (GCA→TCA) ‡uvrCUvrABC system protein C
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (2/3);  total (32/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 = 1.00e+00

CCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCCAG  >  NC_002947/4633664‑4633958
                                                                                                                                                   |                                                                                                                                                   
ccAGCTCCTGCATCACATCGCTGGCCATGATCAGCTGGCCCTTGCCGCCGTCCACCAGCCGCACATCGCCCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGCCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCg                                                                                                                                                    >  5:91109/1‑149 (MQ=255)
ccAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCg                                                                                                                                                    <  2:98585/149‑1 (MQ=255)
  aGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGt                                                                                                                                                  <  2:26903/149‑1 (MQ=255)
       cTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGt                                                                                                                                             >  3:226431/1‑149 (MQ=255)
          cATCACCTCGCGGTCCATGTTCAGCTGGCCCTTGCCGCCGGCCACCAGCAGCACATCGGGCAGGTTGCCCTCGCCGTCCATGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATGGTCACCGCAGGTGACgatt                                                                                                                                           <  4:195990/148‑4 (MQ=255)
          cATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCtt                                                                                                                                          <  4:186018/149‑1 (MQ=255)
          cATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCGCAGGTGACgatt                                                                                                                                           >  3:195990/1‑145 (MQ=255)
           aTCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTc                                                                                                                                         <  7:82781/149‑1 (MQ=255)
              acaTCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCTTCACCTGCGGTGACGTCTTCGa                                                                                                                                       <  6:178583/148‑1 (MQ=255)
                aTCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGt                                                                                                                                    <  3:68743/149‑1 (MQ=255)
                 tCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGAGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGtt                                                                                                                                   >  4:191214/1‑149 (MQ=255)
                   gcgGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTg                                                                                                                                  <  1:138643/148‑1 (MQ=255)
                       gCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAgcg                                                                                                                             >  7:140966/1‑149 (MQ=255)
                            gTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGtt                                                                                                                                   <  6:61052/138‑1 (MQ=255)
                            gTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGtt                                                                                                                                   >  5:61052/1‑138 (MQ=255)
                             ttCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTa                                                                                                                       <  2:173963/149‑1 (MQ=255)
                                   tGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCgta                                                                                                                 >  8:7923/1‑147 (MQ=255)
                                     gCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGAc                                                                                                                <  6:199954/148‑1 (MQ=255)
                                     gCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACt                                                                                                               <  1:141281/149‑1 (MQ=255)
                                     gCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGGCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGACGCGGCGGTAGTCGGAc                                                                                                                <  6:30552/148‑1 (MQ=255)
                                       ccTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACtt                                                                                                              >  1:224709/1‑148 (MQ=255)
                                           gccgccGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATc                                                                                                          <  8:140966/148‑1 (MQ=255)
                                           gccgccGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACATCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATc                                                                                                          <  7:234654/148‑1 (MQ=255)
                                            ccgccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCg                                                                                                         <  5:55800/148‑1 (MQ=255)
                                               ccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGccc                                                                                                     <  3:191214/149‑1 (MQ=255)
                                                  tCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGCCGTCTTCGATGTTGAAGCGGCGGTAGTCTGACTTGATAGGCCCTTc                                                                                                  >  4:197131/1‑149 (MQ=255)
                                                       cagcagCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAgcggc                                                                                                                           <  6:29349/119‑1 (MQ=255)
                                                       cagcagCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAgcggc                                                                                                                           >  5:29349/1‑119 (MQ=255)
                                                       cagcagCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGc                                                                                             <  4:226431/149‑1 (MQ=255)
                                                                aTCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAAcacca                                                                                    >  6:31323/1‑149 (MQ=255)
                                                                    ggCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGc                                                                                 >  3:47983/1‑148 (MQ=255)
                                                                    ggCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGc                                                                                 >  8:208077/1‑148 (MQ=255)
                                                                    ggCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGc                                                                                 >  2:50395/1‑148 (MQ=255)
                                                                        gCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTg                                                                            >  2:210127/1‑149 (MQ=255)
                                                                         cTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTc                                                                                                                    >  7:167282/1‑108 (MQ=255)
                                                                         cTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTc                                                                                                                    <  8:167282/108‑1 (MQ=255)
                                                                          ttGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGc                                                                          >  3:315541/1‑149 (MQ=255)
                                                                               gtCGCGGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                               <  2:126947/122‑1 (MQ=255)
                                                                               cTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                               >  1:126947/1‑123 (MQ=255)
                                                                                          tGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGAGGTAGTCGTACTTTATCGGCCCTTCGGGGCCGAACACCACGCACCTGGCCACGGTAGCTTcacca                                                          >  1:99280/1‑149 (MQ=255)
                                                                                             tGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACgg                                                                     >  3:86199/1‑135 (MQ=255)
                                                                                             tGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACgg                                                                     <  4:86199/135‑1 (MQ=255)
                                                                                                   cATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACg                                                                                  <  6:317317/116‑1 (MQ=255)
                                                                                                   cATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACg                                                                                  >  5:317317/1‑116 (MQ=255)
                                                                                                    aTAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGt                                                   <  7:201456/146‑1 (MQ=255)
                                                                                                    aTAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGt                                                   >  8:201456/1‑146 (MQ=255)
                                                                                                      aGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCt                                               >  4:179927/1‑148 (MQ=255)
                                                                                                         ggcggGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATg                                           <  2:224709/149‑1 (MQ=255)
                                                                                                          gcggGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGt                                          >  4:196000/1‑149 (MQ=255)
                                                                                                                   ggCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAaca                                                                                       <  6:66282/95‑1 (MQ=255)
                                                                                                                   ggCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAaca                                                                                       >  5:66282/1‑95 (MQ=255)
                                                                                                                      cTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCa                              <  1:210127/149‑1 (MQ=255)
                                                                                                                        gATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAgg                            >  3:251211/1‑149 (MQ=255)
                                                                                                                         aTGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAg                             <  7:208077/147‑1 (MQ=255)
                                                                                                                         aTGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGc                           >  1:78954/1‑149 (MQ=255)
                                                                                                                             atcgCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCATCTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTg                       <  5:241510/146‑1 (MQ=255)
                                                                                                                                gCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGGGGTAGTCGGACTTGATCGGCCCTTCGGGGGCGAACACGACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGGATTCCAGGCGTTGCgg                    >  1:226384/1‑149 (MQ=255)
                                                                                                                                   gCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCAc                 <  4:315541/149‑1 (MQ=255)
                                                                                                                                    cATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGc                           <  8:288966/138‑1 (MQ=255)
                                                                                                                                    cATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGc                           >  7:288966/1‑138 (MQ=255)
                                                                                                                                      aagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCgg                    >  5:128424/5‑143 (MQ=255)
                                                                                                                                      aagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCgg                    <  6:128424/139‑1 (MQ=255)
                                                                                                                                         tCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGt                                                   >  7:133482/1‑109 (MQ=255)
                                                                                                                                         tCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGt                                                   <  8:133482/109‑1 (MQ=255)
                                                                                                                                           gTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATcc          <  4:47983/148‑1 (MQ=255)
                                                                                                                                              aCCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGc      >  1:198162/1‑149 (MQ=255)
                                                                                                                                               ccTGAGGTGAAGTAGTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGcc     <  5:287627/149‑1 (MQ=255)
                                                                                                                                                   cGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCCa   <  3:179927/147‑1 (MQ=255)
                                                                                                                                                   cGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCCAg  <  1:50395/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCCAG  >  NC_002947/4633664‑4633958

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