Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,759 T→C 6.1% H567R (CAC→CGC)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,7590TC 6.1% 192.6 / 2.7 65H567R (CAC→CGC) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base C (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

AAATCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGG  >  NC_002947/6104612‑6104901
                                                                                                                                                   |                                                                                                                                               
aaaTCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCgtgt                                                                                                                                                 <  3:337602/148‑1 (MQ=255)
        aGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGa                                                                                                                                        >  2:119250/1‑149 (MQ=255)
        aGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGa                                                                                                                                        >  1:49745/1‑149 (MQ=255)
        aGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGAGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGa                                                                                                                                        >  3:226399/1‑149 (MQ=255)
           cTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                     <  1:119250/149‑1 (MQ=255)
              tttGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGc                                                                                                                                   <  7:56456/148‑1 (MQ=255)
              tttGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATCCGCGTGCACACGCGCATgca                                                                                                                                    <  6:305997/147‑3 (MQ=255)
                                       ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCca‑ccaccac                                                                                                        >  7:60423/1‑149 (MQ=255)
                                       ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCACTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGAt                                                                                                                                       <  6:62978/119‑1 (MQ=255)
                                       ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCACTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGAt                                                                                                                                       >  5:62978/1‑119 (MQ=255)
                                                ctcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGAGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTAc                                                                                               >  3:291223/1‑149 (MQ=255)
                                                 tcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATCCGCGTGCACACGCGCATg                                                                                                                                      <  2:138126/110‑1 (MQ=255)
                                                 tcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATCCGCGTGCACACGCGCATg                                                                                                                                      >  1:138126/1‑110 (MQ=255)
                                                   gggcgggcGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTAcc                                                                                              <  8:60423/147‑1 (MQ=255)
                                                   gggcgggcGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCAc                                                                                            >  8:289033/1‑149 (MQ=255)
                                                          cgAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAAt                                                                                      >  1:47594/1‑148 (MQ=255)
                                                              aTACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCAc                                                                                   >  3:174339/1‑147 (MQ=255)
                                                                  tttGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCa                                                                                             <  2:19326/133‑1 (MQ=255)
                                                                  tttGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCa                                                                                             >  1:19326/1‑133 (MQ=255)
                                                                    tGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCagcagc                                                                            <  4:348681/148‑1 (MQ=255)
                                                                     gCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCaa                                                                          >  5:180893/1‑149 (MQ=255)
                                                                          gggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGc                                                                     >  2:74255/1‑149 (MQ=255)
                                                                            gCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtg                                                                   <  1:306690/149‑1 (MQ=255)
                                                                             cAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtg                                                                   <  1:74255/149‑1 (MQ=255)
                                                                             cAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtgg                                                                  >  7:191094/1‑149 (MQ=255)
                                                                              aGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggt                                                                 >  8:240995/1‑149 (MQ=255)
                                                                                 cccccTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTggg                                                              <  6:180893/149‑1 (MQ=255)
                                                                                   cccTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCa                                                            <  7:289033/149‑1 (MQ=255)
                                                                                    ccTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTggg                                                              >  1:167628/1‑146 (MQ=255)
                                                                                          tGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCt                                                     <  8:191094/149‑1 (MQ=255)
                                                                                          tGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCt                                                     >  4:21559/1‑149 (MQ=255)
                                                                                          tGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCt                                                     >  3:224349/1‑149 (MQ=255)
                                                                                            cGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTgg                                                   >  7:330628/1‑149 (MQ=255)
                                                                                            cGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTgg                                                   <  4:226399/149‑1 (MQ=255)
                                                                                             gAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGc                                                  <  4:48372/149‑1 (MQ=255)
                                                                                             gAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGc                                                  <  7:228657/149‑1 (MQ=255)
                                                                                               gTTCTCGCGACTGTTCTCCAGCGTCTCCAGCCACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCCCAACACCGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGGTGGAGGGCACAATGCTGGCg                                                 >  8:173755/1‑148 (MQ=255)
                                                                                                ttCTCGCGACTGTTCTCCAGAGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGa                                               <  4:291223/149‑1 (MQ=255)
                                                                                                    cgcgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAAt                                           >  1:52061/1‑149 (MQ=255)
                                                                                                     gcgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATg                                          >  6:163267/1‑149 (MQ=255)
                                                                                                      cgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGc                                         <  2:52061/149‑1 (MQ=255)
                                                                                                      cgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGc                                         <  4:174339/149‑1 (MQ=255)
                                                                                                         cTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                      <  8:330628/149‑1 (MQ=255)
                                                                                                          tGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAc                                     >  8:46625/1‑149 (MQ=255)
                                                                                                            ttCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACca                                                                                 <  8:189423/103‑1 (MQ=255)
                                                                                                            ttCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAc                                     <  2:49745/147‑1 (MQ=255)
                                                                                                              ctcCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcag                                   >  5:78955/1‑147 (MQ=255)
                                                                                                               tcCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcag                                >  2:219010/1‑149 (MQ=255)
                                                                                                                 cAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagca                              <  7:240995/149‑1 (MQ=255)
                                                                                                                  aGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGcc                                        <  3:95336/138‑1 (MQ=255)
                                                                                                                  aGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGcc                                        >  4:95336/1‑138 (MQ=255)
                                                                                                                        tcCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCGGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGaaaa                        >  1:68409/1‑148 (MQ=255)
                                                                                                                         ccAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAAt                       >  5:93389/1‑148 (MQ=255)
                                                                                                                         ccAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCAACCACGGCAATCACCAGCAGCAAATCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGGAGGAAAAt                       >  7:228748/1‑148 (MQ=255)
                                                                                                                          ccccAACGCGTTCTGCCTGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACCCCA‑CCACCACCGCATCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATc                      <  8:228748/145‑1 (MQ=255)
                                                                                                                          cAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                      >  6:181644/1‑132 (MQ=255)
                                                                                                                          cAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                      <  5:181644/132‑1 (MQ=255)
                                                                                                                               aCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtg                                                                <  7:50917/101‑1 (MQ=255)
                                                                                                                               aCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtg                                                                >  8:50917/1‑101 (MQ=255)
                                                                                                                               aCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGa                                               >  6:23399/1‑118 (MQ=255)
                                                                                                                               aCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGa                                               <  5:23399/118‑1 (MQ=255)
                                                                                                                                       tGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACGACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGacaa        >  8:93380/1‑149 (MQ=255)
                                                                                                                                         cATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCg            <  6:93389/143‑1 (MQ=255)
                                                                                                                                           tccgcgTGCACACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGc    >  3:185125/3‑149 (MQ=255)
                                                                                                                                            gcgcGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCg   >  3:75787/1‑149 (MQ=255)
                                                                                                                                             cgcgTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCgg  <  4:224349/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AAATCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCA‑CCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGG  >  NC_002947/6104612‑6104901

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