Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,125,766 T→G 7.9% S171A (TCC→GCC)  PP_4512 → DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,125,7660TG 7.9% 206.3 / 8.6 76S171A (TCC→GCC) PP_4512DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (42/28);  new base G (2/4);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

AAGCCGTCTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACGGGTGTG  >  NC_002947/5125620‑5125915
                                                                                                                                                   |                                                                                                                                                     
aaGCCGTCTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCg                                                                                                                                                    >  1:216097/1‑150 (MQ=255)
   ccGTCTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCAATCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcc                                                                                                                                                 >  2:810780/1‑150 (MQ=255)
       cTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcccgat                                                                                                                                             <  2:1114509/150‑5 (MQ=255)
       cTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcccgat                                                                                                                                             <  2:558819/150‑5 (MQ=255)
          ggTCGTGCTGCTCCCCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTTTCTGCACTGCTCTCCGGCCTGCTCTTCACGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcccgattgg                                                                                                                                          <  1:810780/150‑8 (MQ=255)
              gTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCTGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCa                                                                                                                                      <  1:359523/150‑1 (MQ=255)
              gTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACAAATCGGGCCGGCTCCGGCCGTTTGCCTGCa                                                                                                                                      <  2:1180547/150‑1 (MQ=255)
               tGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGCGCCACACTCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAg                                                                                                                                     >  2:748652/1‑150 (MQ=255)
               tGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGCCCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAg                                                                                                                                     <  2:1157481/150‑1 (MQ=255)
                   ggTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTAGCctacactaca                                                                                                                                 <  2:1064294/150‑1 (MQ=255)
                    gTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAg                                                                                                                                >  2:844698/1‑150 (MQ=255)
                    gTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAg                                                                                                                                >  1:452474/1‑150 (MQ=255)
                     tCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAgg                                                                                                                               <  2:1249738/150‑1 (MQ=255)
                          aTCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCa                                                                                                                          >  1:1073303/1‑150 (MQ=255)
                           tCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCAc                                                                                                                         >  1:1115637/1‑150 (MQ=255)
                              gCCTGATACTGGTCACACCAGCATTCCACCTGGCCAGCGAAGCCACCGGCGCCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGCTTGCCTGCAGTACAGGCCGCACtgt                                                                                                                      >  1:628269/1‑150 (MQ=255)
                               ccTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACtgtg                                                                                                                     <  2:452474/150‑1 (MQ=255)
                                         gTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           >  2:783066/1‑150 (MQ=255)
                                         gTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           <  1:753582/150‑1 (MQ=255)
                                          acacaCCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCCGCTATGGAAAACGTTCTCGGGCCTGCCGCACGCGCTGCTGTCACTGACAAATCGTGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtgg                                                                                                          >  1:739007/2‑150 (MQ=255)
                                              acCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtggtggt                                                                                                      <  2:819050/150‑1 (MQ=255)
                                              acCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGCAGCGCAGCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAAGCCGGAATGTGGCAAGACCtgttggt                                                                                                      >  1:798449/1‑150 (MQ=255)
                                              acCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCCGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtggtggt                                                                                                      >  2:971927/1‑150 (MQ=255)
                                                cAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCCGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtggtggtgg                                                                                                    <  1:971927/150‑1 (MQ=255)
                                                cAGCATTCGACCTGGCCAGCGAAGCAACCGGCGGCGTACTAAGGGGCCAGGTCTCGGGCGGACTGTTTGCACTTCTGGGACTTACGAGTCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtggtggtgg                                                                                                    >  2:1231227/1‑150 (MQ=255)
                                                           cTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           >  2:1081358/1‑132 (MQ=255)
                                                           cTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           <  1:1081358/132‑1 (MQ=255)
                                                             ggCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCctgct                                                                                       >  2:738292/1‑150 (MQ=255)
                                                                cAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACtgtg                                                                                                                     <  1:547583/117‑1 (MQ=255)
                                                                cAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACtgtg                                                                                                                     >  2:547583/1‑117 (MQ=255)
                                                                  gCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcccgat                                                                                                                                             >  1:793736/1‑87 (MQ=255)
                                                                  gCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGAGCCGGcccgat                                                                                                                                             <  2:793736/91‑5 (MQ=255)
                                                                    gAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGt                                                                                  >  1:1047903/1‑148 (MQ=255)
                                                                    gAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTc                                                                                >  1:347864/1‑150 (MQ=255)
                                                                    gAAGCCACCGGCGGCCTGCGCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGt                                                                                  <  2:1047903/148‑1 (MQ=255)
                                                                         cACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAgc                                                                           >  2:1209125/1‑150 (MQ=255)
                                                                         cACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAgc                                                                           >  1:1068288/1‑150 (MQ=255)
                                                                          aCCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAgcg                                                                          >  1:1403011/1‑150 (MQ=255)
                                                                                        ctGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGa                                                            <  1:197332/150‑1 (MQ=255)
                                                                                        cggggCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGGGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGa                                                            <  2:628269/148‑1 (MQ=255)
                                                                                         tGGGCACTGCCACCGGGCCTGCGGTTCTCGCTGCTGTCACTGACCACTCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGAc                                                           <  2:442773/150‑1 (MQ=255)
                                                                                                tGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCgg                                                    >  2:65808/1‑150 (MQ=255)
                                                                                                 gCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCtgt                                                                                              <  1:273492/107‑1 (MQ=255)
                                                                                                 gCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCtgt                                                                                              >  2:273492/1‑107 (MQ=255)
                                                                                                   tctcGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGccc                                                 >  2:375791/1‑150 (MQ=255)
                                                                                                    ctcGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCt                                                <  2:347864/150‑1 (MQ=255)
                                                                                                     tcGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGccc                                                 <  2:801174/148‑1 (MQ=255)
                                                                                                     tcGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGccc                                                 >  1:801174/1‑148 (MQ=255)
                                                                                                     tcGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTg                                               >  1:252602/1‑150 (MQ=255)
                                                                                                             gctgTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCt                                       >  2:31944/1‑150 (MQ=255)
                                                                                                              ctgTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCtg                                      >  2:268918/1‑150 (MQ=255)
                                                                                                               tgTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCtgt                                     <  2:1073303/150‑1 (MQ=255)
                                                                                                               tgTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCtgt                                     <  1:424009/150‑1 (MQ=255)
                                                                                                                gTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTTGTTGGCCTGTGCCTGCTGCCGTTCGCAGCGCCAGGCAAAGTCCAAGTTGCGGATCTGGGCTGGGtgtg                                    >  1:1022550/1‑150 (MQ=255)
                                                                                                                     tcGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           <  2:305962/75‑1 (MQ=255)
                                                                                                                     tcGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCtg                                                                                                           >  1:305962/1‑75 (MQ=255)
                                                                                                                     tcGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATc                                >  2:778615/1‑150 (MQ=255)
                                                                                                                     tcGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATc                                >  2:1061948/1‑150 (MQ=255)
                                                                                                                     tcGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATc                                >  2:1258297/1‑150 (MQ=255)
                                                                                                                      cGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCg                               <  1:31944/150‑1 (MQ=255)
                                                                                                                       gctgctGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGc                              >  2:600359/1‑150 (MQ=255)
                                                                                                                        ctgctgTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGcc                             <  1:600359/150‑1 (MQ=255)
                                                                                                                         tgctgTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGa                                             <  2:1068157/133‑1 (MQ=255)
                                                                                                                         tgctgTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGa                                             >  1:1068157/1‑133 (MQ=255)
                                                                                                                           ctgTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCt                                                >  2:424466/1‑128 (MQ=255)
                                                                                                                           ctgTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCt                                                <  1:424466/128‑1 (MQ=255)
                                                                                                                                aCTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTAGTGGGCCTGTGCCTGCTGACGTTCGCAGCGCCCGGCATGGGGGCAGTTGCGGCCCCGGCCTGAGGGTAGATCGCCCGGCTCgg                     >  2:85723/1‑150 (MQ=255)
                                                                                                                                  tGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCCGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCa                   >  2:1054193/1‑150 (MQ=255)
                                                                                                                                    aCCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCAtc                 >  1:976635/1‑150 (MQ=255)
                                                                                                                                      cAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCAtctt               >  1:1209068/1‑150 (MQ=255)
                                                                                                                                       aaTCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCAtcttc              <  2:252602/150‑1 (MQ=255)
                                                                                                                                       aaTCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTCGCTGTGGATCGCCCTGCTCGGCAtcttc              >  2:177506/1‑150 (MQ=255)
                                                                                                                                           gggCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACCGCATACTGATCGCCCTGCTCGGCATCTTCTGCa          <  1:748652/150‑1 (MQ=255)
                                                                                                                                              ccGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACgg       <  2:1320961/150‑1 (MQ=255)
                                                                                                                                               cGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACggg      >  1:879137/1‑150 (MQ=255)
                                                                                                                                                 gCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACGGgtg    >  1:840202/1‑150 (MQ=255)
                                                                                                                                                   tCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACGGGtgtg  >  1:142333/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGCCGTCTTGGTCGTGCTGGTCAGCATCGGCCTGATACTGGTCACACCAGCATTCGACCTGGCCAGCGAAGCCACCGGCGGCCTGCTCTGGGCACTGCTCTCGGGCCTGCTGTTCTCGCTGCTGTCACTGACCAATCGGGCCGGCTCCGGCCGTTTGCCTGCAGTACAGGCCGCACTGTGGCAAAACCTGGTGGTGGGCCTGTGCCTGCTGCCGTTCGCAGCGCCCGGCCTGGCCGACGTTGCGGCCCTGGACTGGCTGTGGATCGCCCTGCTCGGCATCTTCTGCACGGGTGTG  >  NC_002947/5125620‑5125915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: