Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 609,783 G→T 6.1% L356M (CTG→ATG) ‡ PP_0525 ← B12 family TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947609,7830GT 6.1% 869.6 / 30.8 260L356M (CTG→ATG) ‡PP_0525B12 family TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (116/128);  new base T (5/11);  total (121/139)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 273 total aligned reads displayed.

ATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCGG  >  NC_002947/609634‑609932
                                                                                                                                                     |                                                                                                                                                     
aTTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAAcagcag                                                                                                                                                       <  2:2796332/150‑1 (MQ=255)
 ttCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAgg                                                                                                                                                      >  2:578358/1‑150 (MQ=255)
   cTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCt                                                                                                                                                    >  1:3020468/1‑150 (MQ=255)
    ttGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtg                                                                                                                                                   >  2:1803284/1‑150 (MQ=255)
      gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                 <  1:850875/150‑1 (MQ=255)
      gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                 <  1:3627473/150‑1 (MQ=255)
      gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                 <  2:2418132/150‑1 (MQ=255)
      gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                 <  1:2351614/150‑1 (MQ=255)
      gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                                 <  1:2881323/150‑1 (MQ=255)
       tCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgttg                                                                                                                                                >  1:5219542/1‑150 (MQ=255)
       tCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgttg                                                                                                                                                <  1:5554977/150‑1 (MQ=255)
        cGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttgg                                                                                                                                               <  1:3053272/150‑2 (MQ=255)
        cGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttgg                                                                                                                                               <  1:1170560/150‑2 (MQ=255)
        cGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgttgt                                                                                                                                               <  2:4761534/150‑1 (MQ=255)
        cGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgttgt                                                                                                                                               <  2:493528/150‑1 (MQ=255)
         gTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTtgtg                                                                                                                                              <  1:708693/150‑1 (MQ=255)
           ggCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGtt                                                                                                                                            >  2:5611908/1‑150 (MQ=255)
           ggCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACACAACAGCCTGCTGTTgggta                                                                                                                                            <  2:6084505/150‑2 (MQ=255)
             cgcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                          <  2:5481913/150‑1 (MQ=255)
                  gCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAg                                                                                                                                     >  1:4175868/1‑150 (MQ=255)
                   ccAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAgg                                                                                                                                    <  1:3700771/150‑1 (MQ=255)
                   ccAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAgg                                                                                                                                    <  1:6033191/150‑1 (MQ=255)
                    cAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                   >  1:2740634/1‑150 (MQ=255)
                    cAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                   >  1:6788543/1‑150 (MQ=255)
                    cAAGCTCGGTGTAAAAGTATTCGTAATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGATCACCCAggg                                                                                                                                   >  2:785458/1‑150 (MQ=255)
                     aaGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                  >  2:5551723/1‑150 (MQ=255)
                        cTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAg                                                                                                                               <  2:1332083/150‑1 (MQ=255)
                        cTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAg                                                                                                                               <  2:1092226/150‑1 (MQ=255)
                         tCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAgg                                                                                                                              <  1:5273150/150‑1 (MQ=255)
                         tCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAgg                                                                                                                              >  2:5684737/1‑150 (MQ=255)
                         tCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAgg                                                                                                                              <  2:4318440/150‑1 (MQ=255)
                                  aaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGc                                                                                                                     <  1:1794875/150‑1 (MQ=255)
                                   aGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGcc                                                                                                                    >  2:2615515/1‑150 (MQ=255)
                                            cATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCtt                                                                                                           <  1:6518731/150‑1 (MQ=255)
                                             aTGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTg                                                                                                          <  1:164325/150‑1 (MQ=255)
                                             aTGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTg                                                                                                          >  1:3670482/1‑150 (MQ=255)
                                               gAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGt                                                                                                        >  1:2285467/1‑150 (MQ=255)
                                               gAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGt                                                                                                        <  2:3371970/150‑1 (MQ=255)
                                                    cGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCgg                                                                                                   <  2:623048/150‑1 (MQ=255)
                                                    cGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCgg                                                                                                   >  2:1610659/1‑150 (MQ=255)
                                                       tGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                                <  1:6298533/150‑1 (MQ=255)
                                                       tGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                                <  2:6091694/150‑1 (MQ=255)
                                                          tGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtaggta                                                                                             >  1:1524048/1‑150 (MQ=255)
                                                              cGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGa                                                                                         >  1:5913131/1‑150 (MQ=255)
                                                              cGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGa                                                                                         <  2:1141853/150‑1 (MQ=255)
                                                                 aCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAg                                                                                      <  2:2710306/150‑1 (MQ=255)
                                                                       cATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTgcgc                                                                                <  2:5019408/150‑1 (MQ=255)
                                                                       cATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTgcgc                                                                                >  2:3488841/1‑150 (MQ=255)
                                                                           aTTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGa                                                                            >  2:4458395/1‑150 (MQ=255)
                                                                           aTTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGa                                                                            <  1:4360969/150‑1 (MQ=255)
                                                                            ttCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAt                                                                                            >  1:4000210/1‑133 (MQ=255)
                                                                            ttCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAt                                                                                            <  2:4000210/133‑1 (MQ=255)
                                                                             tCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                          >  1:789128/1‑150 (MQ=255)
                                                                             tCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGt                                                                          <  2:904985/150‑1 (MQ=255)
                                                                               cAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATTGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTCGGTTGTTCTGACAAGTTGCTTGGTCGCGTTAGGAATTGAAAGGTGCGCATGAGTgg                                                                        >  1:2040456/1‑150 (MQ=255)
                                                                               cAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTgg                                                                        >  1:489322/1‑150 (MQ=255)
                                                                                   ggCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGa                                                                    <  1:946009/150‑1 (MQ=255)
                                                                                     cTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGc                                                                  <  1:1066718/150‑1 (MQ=255)
                                                                                      tGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCt                                                                 >  2:3465400/1‑150 (MQ=255)
                                                                                       gTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCtt                                                                <  1:4624729/150‑1 (MQ=255)
                                                                                         cTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGt                                                              <  1:3235074/150‑1 (MQ=255)
                                                                                            ccGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtctc                                                           >  2:4918052/1‑150 (MQ=255)
                                                                                                  aGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCttt                                                     >  2:3784173/1‑150 (MQ=255)
                                                                                                      cTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATc                                                 <  2:5273451/150‑1 (MQ=255)
                                                                                                       tCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcg                                                >  2:6742596/1‑150 (MQ=255)
                                                                                                       tCACTCGCATGGACACGGTCTTCGAACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGCCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGGAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcg                                                >  2:732960/1‑150 (MQ=255)
                                                                                                         aCTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgcg                                              <  2:4541309/150‑1 (MQ=255)
                                                                                                         aCTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATcgcg                                              <  2:1548804/150‑1 (MQ=255)
                                                                                                           tCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGt                                            <  2:3288518/150‑1 (MQ=255)
                                                                                                             gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGttt                                          >  2:958260/1‑150 (MQ=255)
                                                                                                               aTGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTc                                        >  2:2350347/1‑150 (MQ=255)
                                                                                                                 ggACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGc                                      >  1:3460586/1‑150 (MQ=255)
                                                                                                                  gACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCt                                     <  1:4735673/150‑1 (MQ=255)
                                                                                                                  gACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCt                                     <  1:4001403/150‑1 (MQ=255)
                                                                                                                    cacGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCtgt                                   <  1:5352775/150‑1 (MQ=255)
                                                                                                                     acGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                   <  2:4497194/53‑1 (MQ=255)
                                                                                                                      cGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                                 <  1:3047730/150‑1 (MQ=255)
                                                                                                                      cGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCGTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                                 <  1:5630319/150‑1 (MQ=255)
                                                                                                                          cTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCa                             <  1:3815296/150‑1 (MQ=255)
                                                                                                                           ttCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAg                            >  1:5521108/1‑150 (MQ=255)
                                                                                                                           ttCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAg                            >  1:6564520/1‑150 (MQ=255)
                                                                                                                               tACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCa                        >  2:175710/1‑150 (MQ=255)
                                                                                                                                 ccAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCaaa                      <  1:2615515/150‑1 (MQ=255)
                                                                                                                                     tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                  <  1:1673107/150‑1 (MQ=255)
                                                                                                                                     tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                  <  2:3057229/150‑1 (MQ=255)
                                                                                                                                     tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                  <  2:4595589/150‑1 (MQ=255)
                                                                                                                                     tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                  <  1:4406617/150‑1 (MQ=255)
                                                                                                                                       gCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGt                >  1:6442686/1‑150 (MQ=255)
                                                                                                                                        cTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTg               <  2:2116314/150‑1 (MQ=255)
                                                                                                                                         tACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGc              >  1:4990471/1‑150 (MQ=255)
                                                                                                                                          acacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCaa                       <  2:4670224/137‑1 (MQ=255)
                                                                                                                                           cccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc            <  2:4173268/150‑1 (MQ=255)
                                                                                                                                           cccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc            <  2:5856011/150‑1 (MQ=255)
                                                                                                                                           cccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc            <  2:6305922/150‑1 (MQ=255)
                                                                                                                                            ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTcc           >  2:51578/1‑150 (MQ=255)
                                                                                                                                            acAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTcc           >  1:3771835/2‑150 (MQ=255)
                                                                                                                                             cAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCTGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCa          <  2:3589759/150‑1 (MQ=255)
                                                                                                                                             cAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCa          >  2:1522248/1‑150 (MQ=255)
                                                                                                                                                cagcagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTg       <  2:6489499/150‑1 (MQ=255)
                                                                                                                                                   cagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTc    <  2:3020468/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTCGG  >  NC_002947/609634‑609932

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: