Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 994,950 T→A 5.5% intergenic (+26/+99) 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,9500TA 5.5% 607.8 / 16.2 182intergenic (+26/+99)PP_0859/PP_0860ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (82/90);  new base A (7/3);  total (89/93)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 186 total aligned reads displayed.

GAGCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994801‑995099
                                                                                                                                                     |                                                                                                                                                     
gagCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCt                                                                                                                                                       <  1:1397484/150‑1 (MQ=255)
     tGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGccc                                                                                                                                                  >  1:2619522/1‑150 (MQ=255)
     tGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGccc                                                                                                                                                  <  2:2664809/150‑1 (MQ=255)
        tCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCcgcg                                                                                                                                               >  1:1187180/1‑150 (MQ=255)
           gTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTAGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:2475094/150‑1 (MQ=255)
            tGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           >  2:3626039/1‑150 (MQ=255)
            tGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           <  2:2618876/150‑1 (MQ=255)
            tGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           >  1:4736042/1‑150 (MQ=255)
             gCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTTTTCGCGGGCAAGCCCGCGaagaa                                                                                                                                          <  2:2805678/150‑4 (MQ=255)
             gCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGaaaaa                                                                                                                                          <  2:3004009/150‑3 (MQ=255)
             gCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                          <  2:334419/150‑1 (MQ=255)
               cGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                        <  2:3686147/150‑1 (MQ=255)
                    aGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                   >  2:4213435/1‑150 (MQ=255)
                     ggCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGccc                                                                                                                                  <  2:2758992/150‑1 (MQ=255)
                      gCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCt                                                                                                                                 <  2:2593150/150‑1 (MQ=255)
                         gACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgcc                                                                                                                              >  1:2292272/1‑150 (MQ=255)
                          aCGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgcca                                                                                                                             <  2:544252/150‑1 (MQ=255)
                                   tCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                                    >  2:2974138/1‑150 (MQ=255)
                                     aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                  >  2:3813076/1‑150 (MQ=255)
                                     aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCa                                                                                                                  >  2:422042/1‑150 (MQ=255)
                                       tGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCAcc                                                                                                                >  2:4560468/1‑150 (MQ=255)
                                       tGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCAcc                                                                                                                <  2:5064942/150‑1 (MQ=255)
                                             aCTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGt                                                                                                          <  1:4981193/150‑1 (MQ=255)
                                                     cggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGtt                                                                                                  <  1:2545420/150‑1 (MQ=255)
                                                     cggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGtt                                                                                                  <  1:364085/150‑1 (MQ=255)
                                                       gcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgc                                                                                                <  1:4995709/150‑1 (MQ=255)
                                                          tgtgtattgatcggtgtgtaggagcgggtttacccgcgaatacggtggcggcggcaacggtgaacggccggtggaaattggctggcaGGGCTTGCCCGCGAAAGGGCCGGCCCTgcca                                                                                                                             <  2:1151817‑M1/31‑1 (MQ=255)
                                                              tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggcc                                                                                                                                       >  1:3829807/1‑104 (MQ=255)
                                                              tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgcc                                                                                         >  1:891266/1‑150 (MQ=255)
                                                                gcggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAc                                                                                       >  2:504638/1‑150 (MQ=255)
                                                                gcggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAc                                                                                       >  1:4568874/1‑150 (MQ=255)
                                                                 cggcgTCCCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCCCTGACAGCAGGCCCGGCTTCTTCGCGGGCTGGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCAcc                                                                                      <  1:1040511/150‑1 (MQ=255)
                                                                         cGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgc                                                                              <  2:4700714/150‑1 (MQ=255)
                                                                         cGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgc                                                                              <  1:782427/150‑1 (MQ=255)
                                                                          gggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgcg                                                                             >  2:578228/1‑150 (MQ=255)
                                                                           ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  2:4175211/150‑1 (MQ=255)
                                                                           ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGGTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                            <  1:4590790/150‑1 (MQ=255)
                                                                                aGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTaaa                                                                       <  1:3389421/150‑1 (MQ=255)
                                                                                aGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTaaa                                                                       >  2:2964880/1‑150 (MQ=255)
                                                                                 gTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAc                                                                      <  1:578228/150‑1 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                                  <  1:4213435/150‑1 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                                  <  1:4213710/150‑1 (MQ=255)
                                                                                       ggCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAcc                                                                     >  1:3531851/1‑145 (MQ=255)
                                                                                       ggCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAcc                                                                     <  2:3531851/145‑1 (MQ=255)
                                                                                         ccAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCGGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCt                                                              <  2:2940425/150‑1 (MQ=255)
                                                                                          cAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGTAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTCCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  2:2384952/1‑150 (MQ=255)
                                                                                          cAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             <  2:2302074/150‑1 (MQ=255)
                                                                                                 gATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCg                                                      >  2:952752/1‑150 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                   <  1:2974138/150‑1 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCGCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 >  1:1338864/1‑150 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 <  2:4212713/150‑1 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 >  1:2418066/1‑150 (MQ=255)
                                                                                                       gATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAAt                                                >  1:5070234/1‑150 (MQ=255)
                                                                                                         tACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATac                                              <  2:1898069/150‑1 (MQ=255)
                                                                                                          aCCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                             <  2:4055633/150‑1 (MQ=255)
                                                                                                           ccTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                            <  2:1050303/150‑1 (MQ=255)
                                                                                                           ccTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                            <  2:1049962/150‑1 (MQ=255)
                                                                                                            cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCCATAcaca                                           >  2:3905729/1‑150 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          <  2:4525099/150‑1 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  2:4033042/1‑150 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  1:2384045/1‑150 (MQ=255)
                                                                                                                gAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAATCCCGCTCCTACACACCGATCAATACACAtctc                                       <  2:4248689/150‑1 (MQ=255)
                                                                                                                   cTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                                    <  1:3592061/150‑1 (MQ=255)
                                                                                                                     gCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGa                                  >  2:2289763/1‑150 (MQ=255)
                                                                                                                      cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGat                                 >  1:2172716/1‑150 (MQ=255)
                                                                                                                        cTGGCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGatag                               >  2:4649807/1‑149 (MQ=255)
                                                                                                                          gACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCg                             <  2:891266/150‑1 (MQ=255)
                                                                                                                           aCAGCAGGCCCGTCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                            <  1:4870805/150‑1 (MQ=255)
                                                                                                                           aCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                            <  2:4770739/150‑1 (MQ=255)
                                                                                                                           aCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                            <  1:1122034/150‑1 (MQ=255)
                                                                                                                             ggcagGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCCCTCCTACACACCCAACAATACACATCTCGCCGATATCGccc                          >  1:114601/2‑150 (MQ=21)
                                                                                                                             agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgc                                                                                          <  1:1601278/86‑1 (MQ=255)
                                                                                                                             agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgc                                                                                          >  2:1601278/1‑86 (MQ=255)
                                                                                                                             agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacacac                                                        <  2:1412149/120‑1 (MQ=255)
                                                                                                                                agGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATc                       >  2:803930/1‑150 (MQ=255)
                                                                                                                                  gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt                     >  1:3060206/1‑150 (MQ=255)
                                                                                                                                    ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCg                   >  2:2128681/1‑150 (MQ=255)
                                                                                                                                    ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCg                   <  1:1856373/150‑1 (MQ=255)
                                                                                                                                    ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCg                   >  2:155403/1‑150 (MQ=255)
                                                                                                                                     cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCACCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa                  >  1:4405564/1‑150 (MQ=255)
                                                                                                                                     cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGGTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCCATCAATACACATCTCGCCGATATCGCCCATCCTCGa                  >  2:4963519/1‑150 (MQ=255)
                                                                                                                                      ggccctttcGCGGGCAAGCCCGCGAAGAAGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt                 >  2:5247204/7‑150 (MQ=25)
                                                                                                                                       gCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATc                <  1:2834863/150‑1 (MQ=255)
                                                                                                                                         cctttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa              >  1:746946/4‑150 (MQ=255)
                                                                                                                                          ctttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAc             >  1:69432/3‑150 (MQ=255)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAGGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaaat            >  1:825240/2‑148 (MQ=255)
                                                                                                                                            ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCACTACACCTCTCGCCGATATCGCCCATCCTCGCTCAAcgg           >  2:4142961/1‑148 (MQ=255)
                                                                                                                                            ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg           >  2:4113640/1‑150 (MQ=255)
                                                                                                                                            ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg           >  1:594960/1‑150 (MQ=255)
                                                                                                                                             tcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgc          <  1:1835234/150‑1 (MQ=255)
                                                                                                                                             tcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgc          >  1:1959738/1‑150 (MQ=255)
                                                                                                                                              cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct         <  1:3083773/150‑1 (MQ=255)
                                                                                                                                              cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct         <  1:3221940/150‑1 (MQ=255)
                                                                                                                                              cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct         <  1:4510380/150‑1 (MQ=255)
                                                                                                                                               gcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  1:3022988/1‑97 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGt       >  2:1623921/1‑150 (MQ=255)
                                                                                                                                                 GGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             <  2:3027325‑M1/95‑1 (MQ=255)
                                                                                                                                                 GGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  1:3027325‑M1/1‑95 (MQ=255)
                                                                                                                                                   GCTTGCCCGCGCAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACa    >  2:1804978‑M1/1‑150 (MQ=255)
                                                                                                                                                   GCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACa    <  1:5265318‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994801‑995099

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: