Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,430,085 G→T 8.2% intergenic (+38/‑146) PP_3924 → / → PP_3925 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,430,0850GT 8.2% 216.2 / 8.4 73intergenic (+38/‑146)PP_3924/PP_3925hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (3/3);  total (34/39)
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.18e-01

TAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGAC  >  NC_002947/4429944‑4430228
                                                                                                                                             |                                                                                                                                                
tAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCtt                                                                                                                                           >  4:172131/1‑149 (MQ=255)
   gcggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCgg                                                                                                                                        <  6:315927/149‑1 (MQ=255)
    cggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGt                                                                                                                                       <  3:172131/149‑1 (MQ=255)
     ggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtt                                                                                                                                      >  3:16203/1‑149 (MQ=255)
        gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttt                                                                                                                                   <  2:106912/149‑1 (MQ=255)
           aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGTGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                <  3:211283/149‑1 (MQ=255)
           aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                <  2:253337/149‑1 (MQ=255)
            ccGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTg                                                                                                                               >  3:169893/1‑149 (MQ=255)
             cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCgg                                                                                                                                        >  3:59313/1‑139 (MQ=255)
             cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCgg                                                                                                                                        <  4:59313/139‑1 (MQ=255)
               aCGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATTCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttt                                                                                                                                  <  7:263718/143‑1 (MQ=255)
                  gCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTg                                                                                                                          >  1:113529/1‑148 (MQ=255)
                          gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                 <  1:196849/149‑1 (MQ=255)
                           ccTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                 >  7:142770/1‑148 (MQ=255)
                            cTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                 >  1:205100/1‑147 (MQ=255)
                               cATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGt                                                                                                             <  1:137997/148‑1 (MQ=255)
                                 tGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTTTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGAt                                                                                                          <  3:232711/149‑1 (MQ=255)
                                 tGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGAt                                                                                                          >  7:234194/1‑149 (MQ=255)
                                   ggTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATc                                                                                                         <  8:116181/148‑1 (MQ=255)
                                     tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccac                                                                                                      <  1:52169/149‑1 (MQ=255)
                                     tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccac                                                                                                      <  4:169893/149‑1 (MQ=255)
                                      ggCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccac                                                                                                      >  1:168814/1‑148 (MQ=255)
                                       gCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                     <  2:168814/148‑1 (MQ=255)
                                               cGGGGTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACAcc                                                                                            <  4:97300/149‑1 (MQ=255)
                                                        ggCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCAt                                                                                                                       >  4:273963/1‑113 (MQ=255)
                                                        ggCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCAt                                                                                                                       <  3:273963/113‑1 (MQ=255)
                                                         gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGc                                                                                                  <  6:56157/133‑1 (MQ=255)
                                                         gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGc                                                                                                  >  5:56157/1‑133 (MQ=255)
                                                           atcATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCTGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                <  1:299357/146‑1 (MQ=255)
                                                              atgatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGca                                                                            <  8:234194/149‑1 (MQ=255)
                                                                gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGcaca                                                                          >  7:241704/1‑149 (MQ=255)
                                                                  tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCt                                                                                                                           <  7:29153/99‑1 (MQ=255)
                                                                  tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCt                                                                                                                           >  8:29153/1‑99 (MQ=255)
                                                                   gCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATgg                                                                       >  7:124555/1‑149 (MQ=255)
                                                                    cTGTAGATGAAAATCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATggg                                                                      <  8:11575/149‑1 (MQ=255)
                                                                      gTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTTCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTAACCTATTTGGCACATGGGa                                                                     >  2:183210/1‑148 (MQ=255)
                                                                       tCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGAc                                                                    <  8:142770/148‑1 (MQ=255)
                                                                              aaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                            >  2:84120/1‑149 (MQ=255)
                                                                               aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                       <  2:227374/122‑1 (MQ=255)
                                                                               aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                            >  8:73593/1‑148 (MQ=255)
                                                                               aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                            >  2:108487/1‑148 (MQ=255)
                                                                               aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGAACGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                       >  1:227374/1‑122 (MQ=255)
                                                                                 tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                           <  4:223988/147‑1 (MQ=255)
                                                                                 tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                           <  8:275019/147‑1 (MQ=255)
                                                                                 tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                           >  3:223988/1‑147 (MQ=255)
                                                                                 tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                           >  7:275019/1‑147 (MQ=255)
                                                                                   aGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCCGCCACACCCTCACCTATGTGGCACATGGGACGATTGTTCCAAcc                                                       >  4:196506/1‑149 (MQ=255)
                                                                                    gctatcTGTTCCCGTCCTGTGGGCTGGGAATACCGGATGATCGAAAAACCGAACCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTACTTTGGCACATGGGACGATTGTTCCAAcc                                                       <  1:183210/144‑1 (MQ=255)
                                                                                     cGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCt                                                     >  2:188968/1‑149 (MQ=255)
                                                                                         cTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCg                                                 >  5:28313/1‑149 (MQ=255)
                                                                                           gTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCCAAGCCGGACGGCTTCGGTTTTTTTTTCCCTGCATTCGGTGCAGTGATCCCCCCGCAACACCCACCACAATTTTGCACAAGGGACGATTGTTACAACCCTGGCggc                                               >  2:320899/1‑148 (MQ=255)
                                                                                            ttCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAc                                              <  4:16203/149‑1 (MQ=255)
                                                                                            ttCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATNCGGTGCAGTGATGCACCAACAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAc                                              >  1:101143/1‑149 (MQ=255)
                                                                                                 gTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                         <  2:205100/149‑1 (MQ=255)
                                                                                                    cTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAc                                              <  2:93078/141‑1 (MQ=255)
                                                                                                    cTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAc                                              >  1:93078/1‑141 (MQ=255)
                                                                                                    cTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                       <  1:320899/148‑1 (MQ=255)
                                                                                                      gAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATg                                                                        >  4:268144/1‑113 (MQ=255)
                                                                                                      gAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATg                                                                        <  3:268144/113‑1 (MQ=255)
                                                                                                        gggCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTAcc                                  <  8:241704/149‑1 (MQ=255)
                                                                                                         ggCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTAccc                                 <  1:108487/149‑1 (MQ=255)
                                                                                                         ggCTGGGAAACCCGGATGATCGAAACACCGAAGCCTGCCGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTAccc                                 <  8:238148/149‑1 (MQ=255)
                                                                                                            tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                               >  4:305425/1‑148 (MQ=255)
                                                                                                                           atCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAGCATTAGATTACCCAATAGTACTTGCTTTc                 >  2:18255/1‑147 (MQ=255)
                                                                                                                            tCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga               <  2:113529/148‑1 (MQ=255)
                                                                                                                              aaaaaaCCGAAGCCGTCAGGCTTCTGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                       <  7:31982/121‑1 (MQ=255)
                                                                                                                              aaaaaaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                       >  8:31982/2‑122 (MQ=39)
                                                                                                                                  aaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAg        >  8:20252/1‑149 (MQ=255)
                                                                                                                                   aCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAg        >  1:40754/1‑148 (MQ=255)
                                                                                                                                   aCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACtt                       <  4:205550/133‑1 (MQ=255)
                                                                                                                                   aCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACtt                       >  3:205550/1‑133 (MQ=255)
                                                                                                                                     cGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGaa     <  8:124555/149‑1 (MQ=255)
                                                                                                                                        aGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGAc  <  7:73593/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
TAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGAC  >  NC_002947/4429944‑4430228

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: