Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,768,986 C→A 7.1% R629L (CGC→CTC)  PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,768,9860CA 7.1% 204.3 / 4.3 65R629L (CGC→CTC) PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base A (3/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GCTTCGAACATCTGCCGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTACAACAGATGAACACCCAACTCGAT  >  NC_002947/5768842‑5769134
                                                                                                                                                |                                                                                                                                                    
gCTTCGAACATCTGCCGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCCGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCt                                                                                                                                                   <  8:83556/148‑1 (MQ=255)
  ttCGAACATCTGCCGATCCACGTCCCGCGCCCACAATTGCAGAGCCAACTGACCAGGCTGCGGTTTCTCGGGGAGGGTGTAAGGGTCGACCAGAGATTGGGGCAGGTCGGCCTCACGGACACCACCATTGTCGGAGACCGCGAGCTGAg                                                                                                                                                >  6:198659/1‑149 (MQ=255)
    cGAACATCTGCCGATCCACGTACCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCt                                                                                                                                              >  7:245338/1‑149 (MQ=255)
        cATCTGCCGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCt                                                                                                                                          <  4:180528/149‑1 (MQ=255)
              ccGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCt                                                                                                                                          <  4:84849/143‑1 (MQ=255)
              ccGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCt                                                                                                                                          >  3:84849/1‑143 (MQ=255)
              ccGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATc                                                                                                                                     <  5:7919/148‑1 (MQ=255)
                 aTCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCca                                                                                                                                   >  2:209460/1‑147 (MQ=255)
                   ccACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCAcag                                                                                                                                >  5:172273/1‑148 (MQ=255)
                   ccACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCAcag                                                                                                                                >  2:204850/1‑148 (MQ=255)
                           cgcgCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGtt                                                                                                                       <  7:83579/149‑1 (MQ=255)
                            gcgcCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTg                                                                                                                      >  6:166468/1‑149 (MQ=255)
                             cgcCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGt                                                                                                                     <  6:218599/149‑1 (MQ=255)
                              gcCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGCAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTc                                                                                                                    <  8:46373/149‑1 (MQ=255)
                               cccAGAATTGCAGCGACAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCCCCATTGTCGTAGACCGCGCCCTGAGCTTGCTCGATCCACAGCAGCACGTTGTcc                                                                                                                   >  6:158212/1‑149 (MQ=255)
                                 cAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTg                                                                                                                 <  4:135786/149‑1 (MQ=255)
                                   gAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGc                                                                                                                >  3:264894/1‑148 (MQ=255)
                                        gCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGcc                                                                                                          >  1:153826/1‑149 (MQ=255)
                                             gcCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCgg                                                                                                     <  6:172273/149‑1 (MQ=255)
                                             gcCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCgg                                                                                                     >  4:137387/1‑149 (MQ=255)
                                                aaCTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGggcg                                                                                                  >  8:166648/1‑149 (MQ=255)
                                                  cTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGggcggc                                                                                                <  6:121711/149‑1 (MQ=255)
                                                  cTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAACCCCCCGggctgc                                                                                                >  2:196472/1‑149 (MQ=255)
                                                         ggCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCg                                                                                          >  8:235133/1‑148 (MQ=255)
                                                            tGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGAGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATc                                                                                      >  3:315232/1‑149 (MQ=255)
                                                            tGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGAGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATc                                                                                      >  3:320234/1‑149 (MQ=255)
                                                               ggTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCaa                                                                                   <  3:103346/149‑1 (MQ=255)
                                                                   tctcCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAg                                                                                <  5:166468/148‑1 (MQ=255)
                                                                    ctcCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTg                                                                              >  6:89668/1‑149 (MQ=255)
                                                                     tcCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGc                                                                             >  7:33410/1‑149 (MQ=255)
                                                                       cGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGcc                                                                            <  7:38299/148‑1 (MQ=255)
                                                                       cGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCt                                                                           <  7:235133/149‑1 (MQ=255)
                                                                       cGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGAGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGc                                                                             <  4:320234/147‑1 (MQ=255)
                                                                       cGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGAGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGc                                                                             <  4:315232/147‑1 (MQ=255)
                                                                         ggAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGa                                                                         <  5:134335/149‑1 (MQ=255)
                                                                                gtAAGGGTCGAGCAGCCATTGGGGCAGGTGTGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGt                                                                  >  2:289882/1‑149 (MQ=255)
                                                                                gtAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGt                                                                  >  5:187189/1‑149 (MQ=255)
                                                                                 tAAGGGTCGAGCAGCCATTGGGGCAGGTGTGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTc                                                                 <  1:289882/149‑1 (MQ=255)
                                                                                 tAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTc                                                                 <  3:56976/149‑1 (MQ=255)
                                                                                   aGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCg                                                                <  1:209460/148‑1 (MQ=255)
                                                                                     ggTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTAc                                                             >  5:100675/1‑149 (MQ=255)
                                                                                          agcagcCCTTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTCCAAAGCCCCCGGGCGGCTGTTCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAAt                                                        <  5:198659/149‑1 (MQ=255)
                                                                                              gcCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAgc                                                    <  8:61329/149‑1 (MQ=255)
                                                                                               ccATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATc                                                                                                                                     <  5:223228/67‑1 (MQ=255)
                                                                                               ccATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATc                                                                                                                                     >  6:223228/1‑67 (MQ=255)
                                                                                                cATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAgcgc                                                  <  5:89668/149‑1 (MQ=255)
                                                                                                 aTTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGcc                                                 >  4:32793/1‑149 (MQ=255)
                                                                                                       gCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCg                                           <  8:33410/149‑1 (MQ=255)
                                                                                                       gCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCg                                           <  4:264894/149‑1 (MQ=255)
                                                                                                           gTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGgtg                                        >  5:350989/1‑148 (MQ=255)
                                                                                                                    ccGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGc                              <  6:350989/149‑1 (MQ=255)
                                                                                                                    ccGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGGCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCggg                               <  8:245338/148‑1 (MQ=255)
                                                                                                                     cGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACt                                                            <  7:172064/118‑1 (MQ=255)
                                                                                                                     cGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACt                                                            >  8:172064/1‑118 (MQ=255)
                                                                                                                     cGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCa                             <  5:158212/149‑1 (MQ=255)
                                                                                                                      ggAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAg                            <  7:166648/149‑1 (MQ=255)
                                                                                                                            caccaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGTTCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTa                          <  3:273236/145‑1 (MQ=255)
                                                                                                                            caccaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGTTCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTa                          >  4:273236/1‑145 (MQ=255)
                                                                                                                            caccaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCggg                               <  1:187181/140‑1 (MQ=255)
                                                                                                                            caccaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTacaac                      <  1:204850/149‑1 (MQ=255)
                                                                                                                            caccaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATAGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCggg                               >  2:187181/1‑140 (MQ=255)
                                                                                                                             accaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGgtgt                                       <  5:77836/131‑1 (MQ=255)
                                                                                                                             accaTTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGgtgt                                       >  6:77836/1‑131 (MQ=255)
                                                                                                                                       tAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAG‑CCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCg                                           >  1:96687/1‑116 (MQ=255)
                                                                                                                                       tAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAG‑CCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCg                                           <  2:96687/116‑1 (MQ=255)
                                                                                                                                             gcgcgcTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAgcgc                                                  <  4:250112/104‑1 (MQ=255)
                                                                                                                                             gcgcgcTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAgcgc                                                  >  3:250112/1‑104 (MQ=255)
                                                                                                                                             gcgcgcTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTACAACAGATGAACACCCAACTc     <  6:187189/149‑1 (MQ=255)
                                                                                                                                                cgcTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTACAACAGATGAACACCCAACTCGAt  <  5:132185/149‑1 (MQ=255)
                                                                                                                                                cgcTGAGCTTGCTCGATCCACAGCAGCACGGTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTCGCGCCTTTTCGGTGTACAGCGGGCAGTACAACAGATGAACACCCAACTCGAt  <  2:153826/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTTCGAACATCTGCCGATCCACGTCCCGCGCCCAGAATTGCAGCGCCAACTGAACAGGCTGCGGTTTCTCCGGGAGGGTGTAAGGGTCGAGCAGCCATTGGGGCAGGTGCGGCTCCCGGAAACCACCATTGTCGTAGACCGCGCGCTGAGCTTGCTCGATCCACAGCAGCACGTTGTCCTGCAAAGCCCCCGGGCGGCTGATCGCATCCAACAGTGCCTGAGCGCTGTCGTACTGAATTAGCGCCTTTTCGGTGTACAGCGGGCAGTACAACAGATGAACACCCAACTCGAT  >  NC_002947/5768842‑5769134

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