Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,430,085 G→T 8.3% 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.3% 218.9 / 9.1 72intergenic (+38/‑146)PP_3924/PP_3925hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base T (2/4);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

CGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCA  >  NC_002947/4429948‑4430231
                                                                                                                                         |                                                                                                                                                  
cggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGt                                                                                                                                         >  8:18167/1‑149 (MQ=255)
cggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGt                                                                                                                                         <  7:31711/149‑1 (MQ=255)
 ggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGtt                                                                                                                                        <  1:68886/149‑1 (MQ=255)
    gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttt                                                                                                                                     >  4:338811/1‑149 (MQ=255)
    gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGt                                                                                                                                         >  1:13898/1‑145 (MQ=255)
    gAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGt                                                                                                                                         <  2:13898/145‑1 (MQ=255)
       aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                  >  8:53687/1‑149 (MQ=255)
       aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                  <  4:199593/149‑1 (MQ=255)
       aCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtttttttt                                                                                                                                  >  2:331312/1‑149 (MQ=255)
        ccGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTg                                                                                                                                 >  6:198712/1‑149 (MQ=255)
         cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGc                                                                                                                                >  6:230793/1‑149 (MQ=255)
          gACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGcc                                                                                                                               <  7:18167/149‑1 (MQ=255)
              gCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGc                                                                                                                           <  5:198712/149‑1 (MQ=255)
                     ggCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGt                                                                                                                    <  3:338811/149‑1 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                   >  6:335612/1‑149 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                   >  2:222228/1‑149 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                   <  7:224957/149‑1 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                   >  8:224957/1‑149 (MQ=255)
                            aTGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGa                                                                                                             >  5:165633/1‑149 (MQ=255)
                             tGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGTTCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGa                                                                                                             <  1:225621/148‑1 (MQ=255)
                             tGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGAt                                                                                                            <  3:22557/149‑1 (MQ=255)
                                 tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCATGCATTCGGTGCAGTGATccac                                                                                                        <  7:53687/149‑1 (MQ=255)
                                   gCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                       <  1:222228/148‑1 (MQ=255)
                                      aGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGc                                                                                                    >  5:81066/1‑148 (MQ=255)
                                           cGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAacac                                                                                               <  1:331312/148‑1 (MQ=255)
                                              gTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTc                                                                                           <  6:165633/149‑1 (MQ=255)
                                                 gAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAcc                                                                                        >  5:17715/1‑149 (MQ=255)
                                                   aGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtt                                                                                                                                        >  6:315900/1‑99 (MQ=255)
                                                   aGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGtt                                                                                                                                        <  5:315900/99‑1 (MQ=255)
                                                       gACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTg                                                                                  >  6:196569/1‑149 (MQ=255)
                                                            gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGcaca                                                                             >  7:212830/1‑149 (MQ=255)
                                                                     gATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCTCCAGCAACACCCTCACCTATTTGGCACATGGGACGAt                                                                    >  5:276052/1‑149 (MQ=255)
                                                                     gATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGa                                                                     <  5:335612/148‑1 (MQ=255)
                                                                        gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                         <  6:80471/125‑1 (MQ=255)
                                                                        gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                         >  5:80471/1‑125 (MQ=255)
                                                                         aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                >  7:289656/1‑149 (MQ=255)
                                                                          aaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCTCCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                               <  6:276052/149‑1 (MQ=255)
                                                                          aaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                               <  8:212830/149‑1 (MQ=255)
                                                                           aCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                              <  8:289656/149‑1 (MQ=255)
                                                                             tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCaa                                                            <  5:230793/149‑1 (MQ=255)
                                                                                 cGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTAGGTTCAGTGATACAACAGCAACACCCTCACCTATTTGGCCCATGGGACGATTGTTCCAACCCt                                                        >  3:102170/1‑149 (MQ=255)
                                                                                  gAGCTGTACCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTATGCACTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGCTACCATTGTTCCAACCCTg                                                       >  6:80987/1‑149 (MQ=255)
                                                                                       gTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGa                                                  >  1:152777/1‑149 (MQ=255)
                                                                                           ccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAt                                              <  6:17715/149‑1 (MQ=255)
                                                                                            cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                             <  6:81066/149‑1 (MQ=255)
                                                                                                cTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGAt                                         <  5:196569/149‑1 (MQ=255)
                                                                                                   aGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTAc                                      <  2:152777/149‑1 (MQ=255)
                                                                                                       cTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGt                                                                 <  8:227716/118‑1 (MQ=255)
                                                                                                       cTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                                  >  8:257335/1‑149 (MQ=255)
                                                                                                       cTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCAACAGCAACACCCTCACCTATTTGGCACATGGGACGATTGt                                                                 >  7:227716/1‑118 (MQ=255)
                                                                                                        tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGc                                                     >  1:253291/1‑129 (MQ=255)
                                                                                                        tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGc                                                     <  2:253291/129‑1 (MQ=255)
                                                                                                        tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                                  >  2:20270/1‑148 (MQ=255)
                                                                                                           gAATACCGGATGATCGAAAAACCGAAGCGTGAGGGCTTGGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAg                               <  4:102170/148‑1 (MQ=255)
                                                                                                              tACCGGATGATCTAAATACCGAAGCCTGACGGCTTCTGTTTTTTTTGGCCTGCATTCGGTTCAGTGATCCACCAGCAAGACGCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGGGGACCATTAGATTACCCAATAGTAc                            <  5:80987/148‑1 (MQ=255)
                                                                                                                ccGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTg                         >  2:315190/1‑149 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCGCATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAg                               <  7:92903/142‑1 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCGCATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAg                               >  8:92903/1‑142 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCGACCCTGGCGGACCATTAGATTACCCAATAGTAc                            <  1:29412/145‑1 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCGACCCTGGCGGACCATTAGATTACCCAATAGTAc                            >  2:29412/1‑145 (MQ=255)
                                                                                                                       caaaaaaaaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCa                                               >  2:316724/5‑120 (MQ=38)
                                                                                                                       caaaaaaaaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCa                                               <  1:316724/116‑1 (MQ=38)
                                                                                                                       atCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                  <  1:20270/149‑1 (MQ=255)
                                                                                                                       atCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACATATTTTGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATACTACTTTCAATc                    >  4:75635/1‑147 (MQ=255)
                                                                                                                        tCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                            <  2:122542/122‑1 (MQ=255)
                                                                                                                        tCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                            >  1:122542/1‑122 (MQ=255)
                                                                                                                         cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                            <  2:269222/121‑1 (MQ=255)
                                                                                                                         cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                            >  1:269222/1‑121 (MQ=255)
                                                                                                                                 cGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGaa        <  1:315190/149‑1 (MQ=255)
                                                                                                                                       cTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCa  <  7:257335/149‑1 (MQ=255)
                                                                                                                                        tGACGGCTTCGGTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTg                                                       >  5:107871/1‑94 (MQ=255)
                                                                                                                                        tGACGGCTTCGGTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTg                                                       <  6:107871/94‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCA  >  NC_002947/4429948‑4430231

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: