Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,296,911 C→T 8.9% intergenic (+34/‑64) gloA → / → PP_3767 lactoylglutathione lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,296,9110CT 8.9% 182.9 / 6.1 67intergenic (+34/‑64)gloA/PP_3767lactoylglutathione lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base T (4/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATGCCGC  >  NC_002947/4296767‑4297060
                                                                                                                                                 |                                                                                                                                                     
gTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCg                                                                                                                                                  <  1:1203416/150‑1 (MQ=255)
  gCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGc                                                                                                                                                <  2:1063714/150‑1 (MQ=255)
   ccGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGcc                                                                                                                                               >  1:692116/1‑150 (MQ=255)
    cGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGc                                                                                                                                                >  2:224137/1‑148 (MQ=255)
    cGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGc                                                                                                                                                <  1:224137/148‑1 (MQ=255)
       tCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCtt                                                                                                                                           <  2:978859/150‑1 (MQ=255)
            aaaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgt                                                                                                                                      <  1:505086/150‑1 (MQ=255)
             aaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                     <  2:527977/150‑1 (MQ=255)
              aCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTgg                                                                                                                                    <  1:1197601/150‑1 (MQ=255)
               cGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGt                                                                                                                                   <  2:527458/150‑1 (MQ=255)
                  cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                                >  2:185356/1‑150 (MQ=255)
                  cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                                <  1:789531/150‑1 (MQ=255)
                   tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                               >  1:1343542/1‑150 (MQ=255)
                       ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                           >  2:677197/1‑150 (MQ=255)
                       ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCCGTGGTATATGAAg                                                                                                                           <  2:950442/150‑1 (MQ=255)
                        gACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAgg                                                                                                                          >  1:370133/1‑150 (MQ=255)
                            gCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAAc                                                                                                                      >  1:166516/1‑150 (MQ=255)
                             ccGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACt                                                                                                                     >  2:849949/1‑150 (MQ=255)
                                   gAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGTtgtg                                                                                                                                     <  1:286676/128‑1 (MQ=255)
                                   gAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                     >  2:286676/1‑128 (MQ=255)
                                   gAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTca                                                                                                               <  2:678751/150‑1 (MQ=255)
                                     accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCAGCTTTCAGGCCGGGAAATCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                             >  1:728060/1‑150 (MQ=255)
                                     accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                             <  2:1343542/150‑1 (MQ=255)
                                               cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                          >  2:706410/1‑143 (MQ=255)
                                               cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                          <  1:706410/143‑1 (MQ=255)
                                               cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCg                                                                                                   >  1:274503/1‑150 (MQ=255)
                                                 tCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAg                                                                                                 <  1:849949/150‑1 (MQ=255)
                                                                   gTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATg                                                                                     <  2:1310972/144‑1 (MQ=255)
                                                                   gTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATg                                                                                     >  1:1310972/1‑144 (MQ=255)
                                                                       ttGGGTCGAAGTGATCCAGCCAACCGANCTCGAAGGCTAANCGTTCGNCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATAGGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTc                                                                           >  1:1166441/1‑150 (MQ=255)
                                                                       ttGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTc                                                                           >  1:1081428/1‑150 (MQ=255)
                                                                           gTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAAtc                                                                                >  2:5604/1‑141 (MQ=255)
                                                                           gTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtg                                                                       >  1:787630/1‑150 (MQ=255)
                                                                           gTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtg                                                                       >  1:1173655/1‑150 (MQ=255)
                                                                           gTCGAAGTGATCCAGCCAACCGAGCGCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAAtc                                                                                <  1:5604/141‑1 (MQ=255)
                                                                            tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                      <  2:274503/150‑1 (MQ=255)
                                                                            tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCGAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                      <  1:631321/150‑1 (MQ=255)
                                                                              gAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                     <  1:1211527/85‑1 (MQ=255)
                                                                              gAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                     >  2:1211527/1‑85 (MQ=255)
                                                                              gAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGa                                                                    >  2:247525/1‑150 (MQ=255)
                                                                               aaGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAgg                                                                                                                          <  2:978058/95‑1 (MQ=255)
                                                                               aaGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAgg                                                                                                                          >  1:978058/1‑95 (MQ=255)
                                                                                 gTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATa                                                                 <  2:1173655/150‑1 (MQ=255)
                                                                                            aaCCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGcaca                                                      >  2:1048249/1‑150 (MQ=255)
                                                                                            aaCCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCacac                                                     <  2:868873/150‑1 (MQ=255)
                                                                                             aCCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACcaca                                                                                                      >  1:920769/1‑101 (MQ=255)
                                                                                             aCCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACcaca                                                                                                      <  2:920769/101‑1 (MQ=255)
                                                                                                     cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCATGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                          >  1:1376561/1‑105 (MQ=255)
                                                                                                     cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCATGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                          <  2:1376561/105‑1 (MQ=255)
                                                                                                     cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATa                                                                                                           <  1:300757/88‑1 (MQ=255)
                                                                                                     cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATa                                                                                                           >  2:300757/1‑88 (MQ=255)
                                                                                                        aGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                          <  1:1266691/102‑1 (MQ=255)
                                                                                                        aGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                          >  2:1266691/1‑102 (MQ=255)
                                                                                                        aGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCa                                          <  2:1081428/150‑1 (MQ=255)
                                                                                                        aGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCa                                          <  1:185356/150‑1 (MQ=255)
                                                                                                            tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                      >  1:834898/1‑150 (MQ=255)
                                                                                                            tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                      >  2:1096666/1‑150 (MQ=255)
                                                                                                            tAAGCGTTCGTCGCGCCGCGCACGTCAGGCCGGGAACTCCTGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGACTACCGCGAGCACACCCTCCCTGCCAACGt                                      >  2:424152/1‑150 (MQ=255)
                                                                                                                       cgcgGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTTTACAAAGGTAACTTCTTCACATACCACAGAGCAGTAAGGACGATGCAATCACTTCACTTTGGATACCTCGAGCACACCATCACTGCCAGCGTGATGCCCAATc                           <  1:1096666/150‑1 (MQ=255)
                                                                                                                       cgcgGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATc                           <  2:834898/150‑1 (MQ=255)
                                                                                                                                  tCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcc                 >  1:1239738/1‑150 (MQ=255)
                                                                                                                                    aagaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcccc               >  2:62617/4‑150 (MQ=255)
                                                                                                                                      gaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGc                                                                                                    >  1:271360/2‑63 (MQ=38)
                                                                                                                                      gaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGc                                                                                                    <  2:271360/62‑1 (MQ=38)
                                                                                                                                        ggCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCAGACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGcc           >  1:792018/1‑150 (MQ=255)
                                                                                                                                            gggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGccc                               >  1:1348818/1‑126 (MQ=255)
                                                                                                                                            gggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGccc                               <  2:1348818/126‑1 (MQ=255)
                                                                                                                                                 cTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCa                                                                                   <  2:721717/69‑1 (MQ=255)
                                                                                                                                                 cTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCa                                                                                   >  1:721717/1‑69 (MQ=255)
                                                                                                                                                 cTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATgccgc  <  1:677197/150‑1 (MQ=255)
                                                                                                                                                 cTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATgccgc  <  2:787630/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATGCCGC  >  NC_002947/4296767‑4297060

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