Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,131,111 T→A 8.0% intergenic (+23/+89) PP_3634 → / ← PP_3635 GNAT family acetyltransferase/sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,131,1110TA 8.0% 221.2 / 6.2 80intergenic (+23/+89)PP_3634/PP_3635GNAT family acetyltransferase/sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (5/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCC  >  NC_002947/4130967‑4131256
                                                                                                                                                |                                                                                                                                                 
aTCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGccc                                                                                                                                               <  5:3485/149‑1 (MQ=255)
 tCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcg                                                                                                                                              >  5:114341/1‑149 (MQ=255)
   tGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcgcg                                                                                                                                            >  4:289338/1‑149 (MQ=255)
   tGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcgcg                                                                                                                                            >  2:176999/1‑149 (MQ=255)
    gCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATGCGATGGCCTGTTCGCGGGCTTGCCCgcgc                                                                                                                                           >  5:363357/1‑149 (MQ=255)
    gCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcg                                                                                                                                            >  3:352116/1‑148 (MQ=255)
    gCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcg                                                                                                                                            >  3:346299/1‑148 (MQ=255)
     cTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                          <  4:352116/149‑3 (MQ=255)
      ttAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGcc                                                                                                                                          >  6:419582/1‑148 (MQ=255)
        aaCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                          <  4:346299/146‑3 (MQ=255)
         aCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCcaca                                                                                                                                      <  3:269350/149‑1 (MQ=255)
          ccAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAg                                                                                                                                     >  2:101121/1‑149 (MQ=255)
          ccAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAg                                                                                                                                     <  1:101121/149‑1 (MQ=255)
           cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAgg                                                                                                                                    <  3:27819/149‑1 (MQ=255)
           cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAg                                                                                                                                     >  1:318138/1‑148 (MQ=255)
                      ggATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAg                                                                                                                         <  3:289338/149‑1 (MQ=255)
                                      gCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                         <  8:42429/149‑1 (MQ=255)
                                         gAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTg                                                                                                           <  8:367581/144‑1 (MQ=255)
                                         gAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTg                                                                                                           >  7:367581/1‑144 (MQ=255)
                                              gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                 <  8:50490/149‑1 (MQ=255)
                                                tGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTgg                                                                                               <  6:114341/149‑1 (MQ=255)
                                                 gcgctttGAGCGGCTGTATGAGCTGCCCTCTCAAGCCCAGGATTGTGTGCGCCAGCGAATGAACGTGCATTGATCGGATGGCCTGTTGGGGGGCTTCCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTgg                                                                                               <  2:318138/142‑1 (MQ=255)
                                                 gcgcGATGAGCGGCTGTATGTGCTGCCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAAGGTGCATTGATCCGATGGCGTGTTCGCGGTCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTgg                                                                                               <  7:141560/148‑1 (MQ=255)
                                                 gcgcGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTgg                                                                                               >  8:103216/1‑148 (MQ=255)
                                                  cgcgATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGa                                                                                             <  3:309302/149‑1 (MQ=255)
                                                     gATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCg                                                                                          >  3:35030/1‑149 (MQ=255)
                                                         aGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGc                                                                                       >  4:305261/1‑148 (MQ=255)
                                                           cGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAg                                                                                    <  5:14809/149‑1 (MQ=255)
                                                             gCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                            <  3:288394/123‑1 (MQ=255)
                                                             gCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                            >  4:288394/1‑123 (MQ=255)
                                                                 gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCcgc                                                                               >  1:338405/1‑148 (MQ=255)
                                                                    gAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATAc                                                                                                                               >  2:135537/1‑97 (MQ=255)
                                                                       cTGGCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGaagaag                                                                        >  3:53992/1‑149 (MQ=255)
                                                                         gACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                      >  7:59833/1‑149 (MQ=255)
                                                                            cGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCa                                                                   >  7:134013/1‑149 (MQ=255)
                                                                              cTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                 >  2:398152/1‑149 (MQ=255)
                                                                                 aaGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTgg                                                               >  7:278960/1‑148 (MQ=255)
                                                                                 aaGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGc                                                              >  8:25142/1‑149 (MQ=255)
                                                                                 aaGCCCAGCAATGTGGGCGCCAGCCCAACAACGTGCATTGAACCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCATAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGc                                                              <  1:176999/149‑1 (MQ=38)
                                                                                    cccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCtt                                                           >  6:206620/1‑149 (MQ=255)
                                                                                       aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                        <  6:363357/149‑1 (MQ=255)
                                                                                        ggAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                       >  8:101079/1‑149 (MQ=255)
                                                                                        ggAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                       <  6:174315/149‑1 (MQ=255)
                                                                                              tgtgCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                  >  2:58824/1‑148 (MQ=255)
                                                                                                       gcgcAAGAACGTGCATTGATCGGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAgcgc                                        <  2:103062/149‑1 (MQ=255)
                                                                                                       gcgcAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAgcgc                                        >  6:288980/1‑149 (MQ=255)
                                                                                                         gcAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCg                                      <  7:101079/149‑1 (MQ=255)
                                                                                                           atcaaCGTGCATTTATCCTAGGGCGGTTTCGCGGGCTTGCCGGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGAt                                    <  2:338405/146‑1 (MQ=255)
                                                                                                           aaGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCGGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGAt                                    >  7:263796/1‑149 (MQ=255)
                                                                                                            aGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                   <  7:25142/149‑1 (MQ=255)
                                                                                                            aGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                   >  7:221725/1‑149 (MQ=255)
                                                                                                            aGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                   <  8:263796/149‑1 (MQ=255)
                                                                                                                 gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGcc                              <  8:221725/149‑1 (MQ=255)
                                                                                                                   gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                            >  5:389727/1‑69 (MQ=255)
                                                                                                                   gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                            <  6:389727/69‑1 (MQ=255)
                                                                                                                   gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                  <  4:139724/127‑1 (MQ=255)
                                                                                                                   gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                  >  3:139724/1‑127 (MQ=255)
                                                                                                                     aTTGATCCGATGGCCTGTTCGCGCGCTAGCCCGCCCCCACAGGAATACCTTAAGGCCCAAGGCCGGTGATGTACCTCTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCACCCATGCACCCAGCGCCGATCAAGCCCaag                          >  1:349738/1‑147 (MQ=255)
                                                                                                                         aTCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGCGCCCGCGAAGAAGCCCCATTGGCCTTTCTGCCATGCACCCAGCGCCGAGGCAGCCCACgcgc                       >  3:238451/1‑148 (MQ=255)
                                                                                                                          tCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACgcgc                       >  3:377785/1‑147 (MQ=255)
                                                                                                                            cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCtga                   <  8:278960/149‑1 (MQ=255)
                                                                                                                             gATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                  >  5:122142/1‑149 (MQ=255)
                                                                                                                                 gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTTATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCAt               <  8:59833/148‑1 (MQ=255)
                                                                                                                                 gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATc              <  7:103216/149‑1 (MQ=255)
                                                                                                                                     gTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGc          >  4:304203/1‑149 (MQ=255)
                                                                                                                                     gTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGc          <  4:35030/149‑1 (MQ=255)
                                                                                                                                      ttCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAg           <  4:377785/147‑1 (MQ=255)
                                                                                                                                      ttCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca         <  1:58824/149‑1 (MQ=255)
                                                                                                                                      ttCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca         <  5:206620/149‑1 (MQ=255)
                                                                                                                                       tCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                  >  1:303382/1‑107 (MQ=255)
                                                                                                                                       tCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                  <  2:303382/107‑1 (MQ=255)
                                                                                                                                       gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca         >  5:166926/2‑148 (MQ=255)
                                                                                                                                       gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca         >  3:182811/2‑148 (MQ=255)
                                                                                                                                         gcgGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGcaca       >  1:402278/1‑148 (MQ=255)
                                                                                                                                         gcgGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAACGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGcaca       >  5:4972/1‑148 (MQ=255)
                                                                                                                                           gggCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCGCCCAGCGCCGATGCAGCCTACGCGCCAGAGCATCCAGCACAAAc    >  7:151029/1‑149 (MQ=255)
                                                                                                                                            ggCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCGTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAAcc   >  7:221694/1‑149 (MQ=255)
                                                                                                                                             gCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAAccc  <  4:53992/149‑1 (MQ=255)
                                                                                                                                              cTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                       <  4:163497/95‑1 (MQ=255)
                                                                                                                                              cTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                       >  3:163497/1‑95 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCC  >  NC_002947/4130967‑4131256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: