Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,430,086 A→C 8.3% intergenic (+39/‑145) 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,0860AC 8.3% 213.1 / 10.0 72intergenic (+39/‑145)PP_3924/PP_3925hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base C (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.95e-01

CGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATGA  >  NC_002947/4429948‑4430234
                                                                                                                                          |                                                                                                                                                    
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                                                                                                                      >  8:224957/1‑149 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                      >  2:222228/1‑149 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                      <  7:224957/149‑1 (MQ=255)
                      gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTg                                                                                                                      >  6:335612/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                                                                                            >  5:80471/1‑125 (MQ=255)
                                                                        gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                            <  6:80471/125‑1 (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                                  >  8:92903/1‑142 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCGCATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAg                                  <  7:92903/142‑1 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCGACCCTGGCGGACCATTAGATTACCCAATAGTAc                               >  2:29412/1‑145 (MQ=255)
                                                                                                                 cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCGACCCTGGCGGACCATTAGATTACCCAATAGTAc                               <  1:29412/145‑1 (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                                               >  1:122542/1‑122 (MQ=255)
                                                                                                                        tCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                               <  2:122542/122‑1 (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                                                          <  6:107871/94‑1 (MQ=255)
                                                                                                                                        tGACGGCTTCGGTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTg                                                          >  5:107871/1‑94 (MQ=255)
                                                                                                                                          aCGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATGa  <  4:81099/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATGA  >  NC_002947/4429948‑4430234

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: