Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 994,951 T→A 6.6% intergenic (+27/+98) PP_0859 → / ← PP_0860 ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947994,9510TA 6.6% 307.3 / 9.4 107intergenic (+27/+98)PP_0859/PP_0860ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (48/50);  new base A (4/3);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 108 total aligned reads displayed.

GCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994803‑995099
                                                                                                                                                    |                                                                                                                                                     
gCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTg                                                                                                                                                      >  1:1080990/1‑150 (MQ=255)
gCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTg                                                                                                                                                      >  2:1439326/1‑150 (MQ=255)
 ccTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                                     >  2:835548/1‑150 (MQ=255)
  cTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGcc                                                                                                                                                    >  2:2379550/1‑150 (MQ=255)
      tCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAACCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCcgcg                                                                                                                                                <  2:1064015/150‑1 (MQ=255)
       gcaggCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACATCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGa                                                                                                                                               >  2:1556589/5‑150 (MQ=255)
             cGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                         <  1:2217820/150‑1 (MQ=255)
              ggggAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTCCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggcc                                                                                                                                        <  1:513083/150‑1 (MQ=255)
                 gAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggc                                                                                                                                     <  2:1612401/150‑1 (MQ=255)
                 gAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggc                                                                                                                                     <  2:1080990/150‑1 (MQ=255)
                  aGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                    <  1:1974189/150‑1 (MQ=255)
                   ggCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGccc                                                                                                                                   <  2:2234252/150‑1 (MQ=255)
                    gCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCt                                                                                                                                  >  2:2113738/1‑150 (MQ=255)
                          ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                         <  2:2121666/137‑1 (MQ=255)
                          ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                         >  1:2121666/1‑137 (MQ=255)
                          ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTAGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagc                                                                                                                           >  2:192914/1‑150 (MQ=255)
                              tgtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAAt                                                                                                                       >  1:2178482/1‑150 (MQ=255)
                               gtTCACTCCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAAtt                                                                                                                      <  1:49799/150‑1 (MQ=255)
                                   aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                  >  1:1821347/1‑150 (MQ=255)
                                   aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                  >  1:533994/1‑150 (MQ=255)
                                   aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCAATGACAGCCAAtgcgca                                                                                                                  <  2:1900796/150‑6 (MQ=255)
                                    cTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCAc                                                                                                                 <  2:533994/150‑1 (MQ=255)
                                       cAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                                    <  2:1929897/144‑1 (MQ=255)
                                       cAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGACCTGCCAGCCAATTTc                                                                                                                    >  1:1929897/1‑144 (MQ=255)
                                                gcgcGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                                     <  2:2242600/150‑1 (MQ=255)
                                                    ggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATCCCTTCGAGCTGCACTGACAGCAAGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGACCGTTCACCGTTg                                                                                                 >  1:626762/1‑150 (MQ=255)
                                                            tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgcc                                                                                         >  1:425415/1‑150 (MQ=255)
                                                               cggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAcc                                                                                      <  1:835548/150‑1 (MQ=255)
                                                               cggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAcc                                                                                      <  2:2178482/150‑1 (MQ=255)
                                                                ggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCg                                                                                     >  1:582345/1‑150 (MQ=255)
                                                                    tACCGGTCCAAGGTACCGGCCAGTCTGGATGCCGATACACTCGAGCGGCAGTGACAGGAGGCGCGGCTTAGTCGCTGGCTTGCCCTCGAACTGGCCGGCCGTCCCAGCCAAATTCCACCAGCCGTTCACCGTTGCCGCCGCCACCGTAtt                                                                                 >  2:345103/1‑150 (MQ=19)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACCGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  2:1551118/150‑1 (MQ=255)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:1556589/150‑1 (MQ=255)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  2:1336284/150‑1 (MQ=255)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:1127947/150‑1 (MQ=255)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:60605/150‑1 (MQ=255)
                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:561252/150‑1 (MQ=255)
                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                                     <  2:401600/124‑1 (MQ=255)
                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                                     >  1:401600/1‑124 (MQ=255)
                                                                            cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                                    >  1:606732/1‑107 (MQ=255)
                                                                            cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                                    <  2:606732/107‑1 (MQ=255)
                                                                                   ccGGCCCATCTTGATGCCGATACCTTCGCGCTGCACTTACAACAGGGCCGGCCTCTTTCCCGGCTTTCCCGCGAAAAGGCCGGGCCTGCCAGCCCAATTTCCCCGGCCGTTCACCCTTGCCGCCGCCACCGTATTCCCGGGCAAACCCGc                                                                  >  2:2197192/1‑150 (MQ=16)
                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGGAACCGTGTTCGCTCGTAAACCCGc                                                                  >  2:1756989/1‑150 (MQ=37)
                                                                                     ggCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                <  1:2379550/150‑1 (MQ=255)
                                                                                      gCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTTTACCCGcgaa                                                               <  1:2248138/150‑4 (MQ=255)
                                                                                                  gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                   >  2:487360/1‑150 (MQ=255)
                                                                                                   ccGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            >  1:198997/1‑124 (MQ=255)
                                                                                                   ccGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  2:198997/124‑1 (MQ=255)
                                                                                                    cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGCAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 >  2:591932/1‑150 (MQ=255)
                                                                                                         ccTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                            <  1:4919/150‑1 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgc                                                                                                                                >  2:1738927/1‑65 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgc                                                                                                                                <  1:1738927/65‑1 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATa                                               <  1:26695/145‑1 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATa                                               >  2:26695/1‑145 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  1:358660/1‑150 (MQ=255)
                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTGGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  1:466331/1‑150 (MQ=255)
                                                                                                               aGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCg                                      <  1:192914/150‑1 (MQ=255)
                                                                                                                    cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGat                                 >  2:1185246/1‑150 (MQ=255)
                                                                                                                      cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            >  2:2387153/1‑105 (MQ=255)
                                                                                                                      cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:2387153/105‑1 (MQ=255)
                                                                                                                      cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGatat                               >  2:1516218/1‑150 (MQ=255)
                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCATGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATc                              >  2:1790629/1‑150 (MQ=255)
                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCACTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATc                              <  1:2330436/150‑1 (MQ=255)
                                                                                                                         aCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                            <  1:1516218/150‑1 (MQ=255)
                                                                                                                         aCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                            >  2:553634/1‑150 (MQ=255)
                                                                                                                           agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgc                                                                              <  2:1808628/98‑1 (MQ=255)
                                                                                                                           agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgc                                                                              >  1:1808628/1‑98 (MQ=255)
                                                                                                                           agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  1:1318206/1‑115 (MQ=255)
                                                                                                                           agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             <  2:1318206/115‑1 (MQ=255)
                                                                                                                              agGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCAGCTCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATc                       <  1:1229865/150‑1 (MQ=255)
                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          <  2:1793816/129‑1 (MQ=255)
                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  1:1793816/1‑129 (MQ=255)
                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt                     >  1:115489/1‑150 (MQ=255)
                                                                                                                                 cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc                    <  1:1859508/150‑1 (MQ=255)
                                                                                                                                 cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc                    <  2:2304380/150‑1 (MQ=255)
                                                                                                                                 cccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc                    <  1:487360/150‑1 (MQ=255)
                                                                                                                                   cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  2:1394268/1‑107 (MQ=255)
                                                                                                                                   cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             <  1:1394268/107‑1 (MQ=255)
                                                                                                                                   cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCATATCCTCGa                  >  1:1910688/1‑150 (MQ=255)
                                                                                                                                    ggCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                 <  2:425415/150‑1 (MQ=255)
                                                                                                                                    ggCTTCTTCGCGGGCTGGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                 <  2:466331/150‑1 (MQ=255)
                                                                                                                                     gccctttcGCGGGCAAGCCCGCGAAGAAGCCGGGCCTGCTAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATc                >  1:1087557/6‑150 (MQ=21)
                                                                                                                                       ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                 >  2:1668611/1‑147 (MQ=255)
                                                                                                                                       ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                 <  1:1668611/147‑1 (MQ=255)
                                                                                                                                       ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa              >  2:74171/1‑150 (MQ=255)
                                                                                                                                         tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:584892/85‑1 (MQ=18)
                                                                                                                                         tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            >  2:584892/2‑86 (MQ=18)
                                                                                                                                         tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCGACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAca            >  1:1132696/2‑149 (MQ=255)
                                                                                                                                          ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCt                                                              >  1:1600583/1‑99 (MQ=255)
                                                                                                                                          ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg           >  2:2050651/1‑150 (MQ=255)
                                                                                                                                            cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct         <  2:115489/150‑1 (MQ=255)
                                                                                                                                             gcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCGGCCAATTTCCACCGGCCGTTCACCGTTGCCGACGCCACCGTATTCGCGGGGAAACCCGCTCCTCCCCACCGATCAATACACACCTCGCCGACATCGCCCATCCTCCATCAActgctt        >  2:1628611/1‑149 (MQ=21)
                                                                                                                                              cgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGt       >  1:1405433/1‑150 (MQ=255)
                                                                                                                                               gggCTTGCCGGCAAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTa      <  2:626762/150‑1 (MQ=255)
                                                                                                                                               gggCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTa      <  1:2050651/150‑1 (MQ=255)
                                                                                                                                                ggCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAc     <  2:1405433/150‑1 (MQ=255)
                                                                                                                                                 gCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacacac                                                        >  2:2108474/1‑98 (MQ=255)
                                                                                                                                                 gCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacacac                                                        <  1:2108474/98‑1 (MQ=255)
                                                                                                                                                 gCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAca    >  1:290804/1‑150 (MQ=255)
                                                                                                                                                  cTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTCGTGCTGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAcag   <  1:635973/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994803‑995099

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