Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,860 C→T 6.4% 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,301,8600CT 6.4% 732.5 / 23.9 233intergenic (+34/‑64)gloA/PP_3767lactoylglutathione lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (69/149);  new base T (5/10);  total (74/159)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAT  >  NC_002947/4301712‑4302004
                                                                                                                                                    |                                                                                                                                                
ggAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACt                                                                                                                                                 <  1:4146677/150‑1 (MQ=255)
   aGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTccc                                                                                                                                              >  1:1341927/1‑150 (MQ=255)
     tGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCgg                                                                                                                                            <  1:2965600/150‑1 (MQ=255)
      gCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGc                                                                                                                                           <  1:2015784/150‑1 (MQ=255)
       ccGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGcc                                                                                                                                          >  1:2387983/1‑150 (MQ=255)
          ttCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCctga                                                                                                                                       <  1:140351/150‑3 (MQ=255)
            ccAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTg                                                                                                                                     <  2:4922170/150‑1 (MQ=255)
            ccAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTg                                                                                                                                     <  1:4708571/150‑1 (MQ=255)
                aaaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgt                                                                                                                                 <  1:3053797/150‑1 (MQ=255)
                 aaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                <  1:1698068/150‑1 (MQ=255)
                  aCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTgg                                                                                                                               <  1:5209311/150‑1 (MQ=255)
                    gCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGta                                                                                                                             <  2:3872570/150‑1 (MQ=255)
                     ccTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtat                                                                                                                            <  1:4483469/150‑1 (MQ=255)
                      cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                           <  2:3937209/150‑1 (MQ=255)
                      cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                           <  1:809239/150‑1 (MQ=255)
                      cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                           <  2:4374430/150‑1 (MQ=255)
                       tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                          >  1:3620011/1‑150 (MQ=255)
                       tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                          <  2:140274/150‑1 (MQ=255)
                       ggCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                          <  2:1551934/149‑1 (MQ=255)
                        gCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATg                                                                                                                         >  1:1044420/1‑150 (MQ=255)
                        gCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATg                                                                                                                         <  2:4262239/150‑1 (MQ=255)
                           ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  2:2994670/150‑1 (MQ=255)
                           ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  2:3501411/150‑1 (MQ=255)
                           ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  2:620544/150‑1 (MQ=255)
                           ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  1:3926995/150‑1 (MQ=255)
                                 ccGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACt                                                                                                                <  1:2313201/150‑1 (MQ=255)
                                 ccGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACt                                                                                                                <  1:1618355/150‑1 (MQ=255)
                                   gCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGc                                                                                                              <  1:4226/150‑1 (MQ=255)
                                    cATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCt                                                                                                             <  1:5044867/150‑1 (MQ=255)
                                     aTGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCtt                                                                                                            <  2:3662758/150‑1 (MQ=255)
                                      tccgggtccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTc                                                                                                           >  2:1138228/8‑150 (MQ=255)
                                      tGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTc                                                                                                           <  2:13591/150‑1 (MQ=255)
                                       gAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTca                                                                                                          <  2:4305749/150‑1 (MQ=255)
                                        aaCCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcac                                                                                                         <  1:3208051/150‑1 (MQ=255)
                                         accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                        <  1:139240/150‑1 (MQ=255)
                                         accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                        <  2:5191575/150‑1 (MQ=255)
                                         accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                        <  1:4813967/150‑1 (MQ=255)
                                          ccaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGTCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACAt                                                                                                       >  2:4730536/1‑150 (MQ=255)
                                            accGGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                     <  2:5244912/150‑1 (MQ=255)
                                             ccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAcc                                                                                                    <  1:2097369/150‑1 (MQ=255)
                                                    ttCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGACGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGa                                                                                             <  1:575075/150‑1 (MQ=255)
                                                     tCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAg                                                                                            <  1:3104310/150‑1 (MQ=255)
                                                       gTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGt                                                                                          >  1:216104/1‑150 (MQ=255)
                                                       gTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGt                                                                                          >  2:5264368/1‑150 (MQ=255)
                                                           aGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAgg                                                                                      >  2:694301/1‑150 (MQ=255)
                                                           aGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAgg                                                                                      <  1:1132194/150‑1 (MQ=255)
                                                           aGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAgg                                                                                      <  1:988893/150‑1 (MQ=255)
                                                           aGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAgg                                                                                      >  2:1549297/1‑150 (MQ=255)
                                                            ggACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                     <  1:3394/150‑1 (MQ=255)
                                                            ggACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                     <  1:3772855/150‑1 (MQ=255)
                                                                ccGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 <  2:3137683/150‑1 (MQ=255)
                                                                ccGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 <  2:1986091/150‑1 (MQ=255)
                                                                ccGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 <  2:2616099/150‑1 (MQ=255)
                                                                ccGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 >  2:5564043/1‑150 (MQ=255)
                                                                     cGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTTGTATATGCAGGTAACTGCTTCTCATACCACAGCGAGGTAAGGACGATGCAAt                                                                            <  1:3102413/150‑1 (MQ=255)
                                                                     cGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAAt                                                                            <  2:792264/150‑1 (MQ=255)
                                                                        ttATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAAtctc                                                                         >  1:4464042/1‑150 (MQ=255)
                                                                         tATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATctct                                                                        <  2:3102140/150‑1 (MQ=255)
                                                                          aTTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCtt                                                                       <  2:527504/150‑1 (MQ=255)
                                                                           ttGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTc                                                                      <  2:5026634/150‑1 (MQ=255)
                                                                              ggTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACt                                                                   >  2:3608875/1‑150 (MQ=255)
                                                                                tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                 <  1:2219258/150‑1 (MQ=255)
                                                                                tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                 <  1:4795163/150‑1 (MQ=255)
                                                                                  gAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGa                                                               >  2:2321765/1‑150 (MQ=255)
                                                                                  gAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGa                                                               <  1:2372733/150‑1 (MQ=255)
                                                                                   aaGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGaa                                                              >  1:4561444/1‑150 (MQ=255)
                                                                                   aaGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGaa                                                              <  1:5505285/150‑1 (MQ=255)
                                                                                     gTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATa                                                            <  2:3003068/150‑1 (MQ=255)
                                                                                      tGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATAc                                                           <  2:3679859/150‑1 (MQ=255)
                                                                                        aTCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcg                                                         >  2:4624775/1‑150 (MQ=255)
                                                                                          ccAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACCTACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcgcg                                                       >  1:3916496/1‑150 (MQ=255)
                                                                                          ccAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACAAACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcgcg                                                       >  1:243983/1‑150 (MQ=255)
                                                                                           cAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGa                                                      <  2:822703/150‑1 (MQ=255)
                                                                                            aGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAg                                                     <  2:2062704/150‑1 (MQ=255)
                                                                                              ccAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGca                                                   <  2:2776895/150‑1 (MQ=255)
                                                                                                 aCCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCacac                                                <  2:4256094/150‑1 (MQ=255)
                                                                                                        tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                         <  1:1000451/150‑1 (MQ=255)
                                                                                                         cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTg                                        >  1:4302273/1‑150 (MQ=255)
                                                                                                            aGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCa                                     <  2:3330848/150‑1 (MQ=255)
                                                                                                             ggCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAg                                    >  2:5323521/1‑150 (MQ=255)
                                                                                                                tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                 <  2:940326/150‑1 (MQ=255)
                                                                                                                tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                 >  2:506870/1‑150 (MQ=255)
                                                                                                                      ttCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATggt                           <  1:4955542/150‑3 (MQ=255)
                                                                                                                        cgtcgCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGcccc                         <  2:2387983/150‑1 (MQ=255)
                                                                                                                          tcgCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCAt                       >  2:2230170/1‑150 (MQ=255)
                                                                                                                            gcggcgCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATcc                     <  1:1997900/150‑1 (MQ=255)
                                                                                                                             cggcgCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATccc                    <  2:1511323/150‑1 (MQ=255)
                                                                                                                             cggcgCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATccc                    <  1:2871307/150‑1 (MQ=255)
                                                                                                                               gcgcgcTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGa                  <  2:4330529/150‑1 (MQ=255)
                                                                                                                               gcgcgcTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGa                  <  2:4330324/150‑1 (MQ=255)
                                                                                                                                    cTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAc             <  1:1858059/150‑1 (MQ=255)
                                                                                                                                     tCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcc            >  1:670934/1‑150 (MQ=255)
                                                                                                                                      cGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAccc           <  2:399867/150‑1 (MQ=255)
                                                                                                                                       aagaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcccc          >  1:787701/4‑150 (MQ=255)
                                                                                                                                         gaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTg        >  2:2850346/2‑150 (MQ=255)
                                                                                                                                          aGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGc       <  1:122713/150‑1 (MQ=255)
                                                                                                                                           ggCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGcc      <  1:4042735/150‑1 (MQ=255)
                                                                                                                                            gCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGccc     >  2:4654722/1‑150 (MQ=255)
                                                                                                                                            gCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGccc     >  2:5336974/1‑150 (MQ=255)
                                                                                                                                               gggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCCCACCAGCACTGCCAGCGTGCTGCCCCATCCCGACTCACCCCTGCCCTag  >  2:3424310/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAT  >  NC_002947/4301712‑4302004

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