Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,028,713 T→A 5.6% Y31N (TAC→AAC)  PP_3553 → AMP‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,028,7130TA 5.6% 219.7 / 2.5 71Y31N (TAC→AAC) PP_3553AMP‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base A (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CCATCGGGGGCCAGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGATGCTGCCCAACACC  >  NC_002947/4028570‑4028856
                                                                                                                                                 |                                                                                                                                               
ccATCGGGGGCCAGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTAAACgg                                                                                                                                             >  6:197453/1‑148 (MQ=255)
 cATCGGGGGCCAGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTAAACGGCa                                                                                                                                           <  5:197453/149‑1 (MQ=255)
          gctggTGCACAGGGTGCGTCTGCCTCCCTTTCAGGATTACGGCTTGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACGATGTCGCCCTCACCCCGCTGAGCTTCATGGAACGGACCGCTTCCGTTTACGGCAACTACCCg                                                                                                                                   <  2:191713/145‑1 (MQ=255)
           cAGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGc                                                                                                                                 >  2:286911/1‑149 (MQ=255)
            aGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGcc                                                                                                                                >  7:27115/1‑149 (MQ=255)
            aGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGcc                                                                                                                                <  2:241254/149‑1 (MQ=255)
              ctGCACAGGCTGCGTCTGCCTCCCTTACAGGATAACCGCATGTAGATCTATGCCCAGGGCCTGATGCCAGAAGCCGTCAACCATGTCGCACTCACCCCCCTGAGCTTCATAGAACGCACCGCTCCCGTTTACGGCAACTACCCg                                                                                                                                   >  1:191713/2‑144 (MQ=255)
                       cTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACgg                                                                                                                     <  3:300945/149‑1 (MQ=255)
                         gCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGgcg                                                                                                                   <  1:286911/149‑1 (MQ=255)
                               gCCTCCCTTACAGGAATACCGCATGTCGATCTATGACCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGCGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATcc                                                                                                             >  8:116822/1‑149 (MQ=255)
                                  tCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCAc                                                                                                                       <  2:85035/136‑1 (MQ=255)
                                  tCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCAc                                                                                                                       >  1:85035/1‑136 (MQ=255)
                                   cccTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATccgccg                                                                                                         >  8:131052/1‑149 (MQ=255)
                                        aCAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCaa                                                                                                      <  6:169383/147‑1 (MQ=255)
                                        aCAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCaa                                                                                                      >  5:169383/1‑147 (MQ=255)
                                          aGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTgg                                                                                                  <  8:27115/149‑1 (MQ=255)
                                           ggATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCGTCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGc                                                                                                 <  8:71708/149‑1 (MQ=255)
                                            gATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATcc                                                                                                                         >  7:155825/1‑124 (MQ=255)
                                            gATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATcc                                                                                                                         <  8:155825/124‑1 (MQ=255)
                                             aTTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCaa                                                                                               <  3:92303/149‑1 (MQ=255)
                                              ttACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGACCGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGc                                                                                                 >  3:55874/1‑146 (MQ=255)
                                                     aTGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTa                                                                                       <  7:116822/149‑1 (MQ=255)
                                                        tCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTAcc                                                                                     <  2:163844/148‑1 (MQ=255)
                                                          gATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATccgc                                                                                                           >  2:277554/1‑124 (MQ=255)
                                                          gATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATccgc                                                                                                           <  1:277554/124‑1 (MQ=255)
                                                          gATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCACTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAAc                                                                                  >  8:125003/1‑149 (MQ=255)
                                                           aTCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACg                                                                                 >  5:182049/1‑149 (MQ=255)
                                                                tGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTc                                                                            >  1:60712/1‑149 (MQ=255)
                                                                tGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTc                                                                            >  5:42953/1‑149 (MQ=255)
                                                                  cccAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCgg                                                                          >  1:15989/1‑149 (MQ=255)
                                                                   ccAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCCCCCCGCTGAGTCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCgg                                                                          >  1:284547/1‑149 (MQ=255)
                                                                     aGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGgcgc                                                                       >  5:194505/1‑149 (MQ=255)
                                                                       ggCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGcc                                                                                                                 >  6:286896/1‑105 (MQ=255)
                                                                         ccTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTgg                                                                   >  2:296753/1‑149 (MQ=255)
                                                                               gCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAgcg                                                             >  3:183336/1‑149 (MQ=255)
                                                                                 ccGCTGCCGTCAACCATGTCGCCCTCACGCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTAc                                                                                      >  1:215876/1‑123 (MQ=255)
                                                                                 ccGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTAc                                                                                      <  2:215876/122‑1 (MQ=255)
                                                                                    cTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTAAACGGCag                                                                                                                                          <  2:316099/67‑2 (MQ=255)
                                                                                    cTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTAAACGGCag                                                                                                                                          >  1:316099/1‑66 (MQ=255)
                                                                                        cGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTggccg                                                    >  8:344554/1‑149 (MQ=255)
                                                                                         gTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGTAACTACCCgg                                                                                                                                  >  6:27686/1‑70 (MQ=255)
                                                                                         gTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGTAACTACCCgg                                                                                                                                  <  5:27686/70‑1 (MQ=255)
                                                                                         gTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCCGGCCAGGGCGCTggccag                                                   >  5:133997/1‑147 (MQ=255)
                                                                                          tCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTggccggc                                                  >  5:64928/1‑149 (MQ=255)
                                                                                            aaCCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCCGGCCAGGGCGCTggccaggc                                                 <  6:133997/148‑5 (MQ=255)
                                                                                             aCCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGt                                               <  4:183336/149‑1 (MQ=255)
                                                                                             aCCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGt                                               <  6:371037/149‑1 (MQ=255)
                                                                                              ccATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTggc                                                      >  5:1537/1‑141 (MQ=255)
                                                                                              ccATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTggc                                                      <  6:1537/141‑1 (MQ=255)
                                                                                              ccATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTg                                              <  3:157250/149‑1 (MQ=255)
                                                                                                 tGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCa                                                                >  6:187446/1‑128 (MQ=255)
                                                                                                 tGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCa                                                                <  5:187446/128‑1 (MQ=255)
                                                                                                           aCCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGt                                   <  2:60712/147‑1 (MQ=255)
                                                                                                            ccccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGc                                                                  <  6:64679/115‑1 (MQ=255)
                                                                                                            ccccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGc                                                                  >  5:64679/1‑115 (MQ=255)
                                                                                                            ccccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGcc                                                           >  6:278329/1‑122 (MQ=255)
                                                                                                            ccccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGcc                                                           <  5:278329/122‑1 (MQ=255)
                                                                                                            ccccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGc                                <  6:64928/149‑1 (MQ=255)
                                                                                                               ccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAAcc                                                                                         <  8:178403/90‑1 (MQ=255)
                                                                                                               ccGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAAcc                                                                                         >  7:178403/1‑90 (MQ=255)
                                                                                                                     aGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCg                         <  2:284547/148‑1 (MQ=255)
                                                                                                                         ttCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTgg                      >  3:331684/1‑148 (MQ=255)
                                                                                                                          tCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTgg                      >  7:27840/1‑147 (MQ=255)
                                                                                                                            aTCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGt                  >  8:288429/1‑149 (MQ=255)
                                                                                                                             tCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTg                 <  7:344554/149‑1 (MQ=255)
                                                                                                                                   gCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGAtgctg           >  2:362035/1‑149 (MQ=255)
                                                                                                                                     aCCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCTGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGATGCTGcc         <  7:125003/149‑1 (MQ=255)
                                                                                                                                     aCCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGc                                >  2:158115/1‑126 (MQ=255)
                                                                                                                                     aCCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGc                                <  1:158115/126‑1 (MQ=255)
                                                                                                                                     aCCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGATGCTGcc         <  7:131052/149‑1 (MQ=255)
                                                                                                                                            ccGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGATGCTGCCCAACAcc  <  8:27840/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCATCGGGGGCCAGGTGCACAGGCTGCGTCTGCCTCCCTTACAGGATTACCGCATGTCGATCTATGCCCAGGGCCTGATGCCCGCTGCCGTCAACCATGTCGCCCTCACCCCGCTGAGCTTCATCGAACGCACCGCTGCCGTTTACGGCAACTACCCGGCCGTCATCCACGGCGCCATCCGCCGCAACTGGCAAGAAACCTACCAACGCTGTCGGCGCCTGGCCAGCGCCCTGGCCGGCCGTGGCATCGGCCGTGGCGACACCGTGGCAGTGATGCTGCCCAACACC  >  NC_002947/4028570‑4028856

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: