Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 609,782 A→C 7.2% L356R (CTG→CGG) ‡ 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,7820AC 7.2% 460.6 / 16.1 142L356R (CTG→CGG) ‡PP_0525B12 family TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (73/58);  new base C (7/4);  total (80/62)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTC  >  NC_002947/609633‑609930
                                                                                                                                                     |                                                                                                                                                    
gATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAAcagca                                                                                                                                                      <  2:1744800/150‑1 (MQ=255)
 aTTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAAcagcag                                                                                                                                                     <  2:1934214/150‑1 (MQ=255)
  ttCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAgg                                                                                                                                                    >  1:2456500/1‑150 (MQ=255)
   tCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGc                                                                                                                                                   >  2:4864351/1‑150 (MQ=255)
     ttGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtg                                                                                                                                                 >  1:722346/1‑150 (MQ=255)
       gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                               <  1:4512464/150‑1 (MQ=255)
       gTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                               <  2:2353701/150‑1 (MQ=255)
        tCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttg                                                                                                                                              <  1:3053134/150‑1 (MQ=255)
         cGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGTTGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTGCTCACTCGCATGGACAGGGTCTTCGTACCAGTCGCTACCCAACAGCCTGCtgttgg                                                                                                                                             <  1:4533530/150‑2 (MQ=255)
          gTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTtgtg                                                                                                                                            <  1:2068863/150‑1 (MQ=255)
           tGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTgtgt                                                                                                                                           <  1:1351339/150‑1 (MQ=255)
              cgcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACCCGGTCTTCGTACCCGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                        >  2:801687/1‑150 (MQ=255)
              cgcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                        <  2:1030172/150‑1 (MQ=255)
               gcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCAc                                                                                                                                       <  1:4864351/150‑1 (MQ=255)
               gcAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCAc                                                                                                                                       >  1:286704/1‑150 (MQ=255)
                   gCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAg                                                                                                                                   >  1:2762753/1‑150 (MQ=255)
                   gCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAg                                                                                                                                   >  1:2763316/1‑150 (MQ=255)
                   gCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACACAACAGCCTGCTGTTgggta                                                                                                                                          <  1:2894671/143‑2 (MQ=255)
                    ccAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAgg                                                                                                                                  <  2:3309056/150‑1 (MQ=255)
                    ccAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAgg                                                                                                                                  <  2:3309093/150‑1 (MQ=255)
                     cAAGCTCGGTGTAATAGTATTCGTCATGAAAGCGATGATGCACCAACAGCGCATGATAACAGCTGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTAACAGTCGCTACCAAACAGAAGGATGATGTGTTCACCCAggc                                                                                                                                 >  2:3156499/1‑149 (MQ=255)
                     cAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAggg                                                                                                                                 >  2:2185559/1‑150 (MQ=255)
                      aaGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                >  2:1500045/1‑150 (MQ=255)
                      aaGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGc                                                                                                                                >  1:689863/1‑150 (MQ=255)
                                tAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTc                                                                                                                      <  1:1195420/150‑1 (MQ=255)
                                 aaaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCt                                                                                                                     <  2:773525/150‑1 (MQ=255)
                                  aaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTg                                                                                                                    >  2:313161/1‑150 (MQ=255)
                                  aaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTg                                                                                                                    <  1:590892/150‑1 (MQ=255)
                                  aaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTg                                                                                                                    >  2:809778/1‑150 (MQ=255)
                                  aaaGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTg                                                                                                                    >  1:3842966/1‑150 (MQ=255)
                                    aGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGcc                                                                                                                  >  1:5076514/1‑150 (MQ=255)
                                        ttCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGt                                                                                                              >  1:2574919/1‑150 (MQ=255)
                                        ttCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGt                                                                                                              >  1:5002445/1‑150 (MQ=255)
                                          cGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTg                                                                                                            <  2:286704/150‑1 (MQ=255)
                                           gTCATGAAAGCGATGCTGGACGAACAGCGCATGAGTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGc                                                                                                           <  2:3482184/150‑1 (MQ=255)
                                            tCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCt                                                                                                          <  1:4082290/150‑1 (MQ=255)
                                             cATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCtt                                                                                                         <  1:350159/150‑1 (MQ=255)
                                               tGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTgg                                                                                                       >  1:786352/1‑150 (MQ=255)
                                                    gCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCtgtt                                                                                                                                               <  2:4683559/105‑1 (MQ=255)
                                                           tGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtaggta                                                                                           >  1:428431/1‑150 (MQ=255)
                                                             gACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTAtt                                                                                         >  2:4674040/1‑150 (MQ=255)
                                                                gAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGaa                                                                                      >  2:2477883/1‑150 (MQ=255)
                                                                gAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGaa                                                                                      <  1:275320/150‑1 (MQ=255)
                                                                 aaCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGaaa                                                                                     <  1:2329070/150‑1 (MQ=255)
                                                                 aaCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGaaa                                                                                     <  1:4674040/150‑1 (MQ=255)
                                                                  aCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                        >  1:4981296/1‑98 (MQ=255)
                                                                  aCAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCa                                                                                                                                        <  2:4981296/98‑1 (MQ=255)
                                                                         aTGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGGTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCt                                                                             >  1:5189142/1‑150 (MQ=255)
                                                                         aTGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCt                                                                             <  2:3363740/150‑1 (MQ=255)
                                                                         aTGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCt                                                                             >  1:4142419/1‑150 (MQ=255)
                                                                         aTGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCt                                                                             >  1:1164150/1‑150 (MQ=255)
                                                                          tGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCtt                                                                            <  2:4527782/150‑1 (MQ=255)
                                                                           gATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTg                                                                           >  1:3565978/1‑150 (MQ=255)
                                                                            aTTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGa                                                                          <  2:3224239/150‑1 (MQ=255)
                                                                                cAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTgg                                                                      >  2:1483466/1‑150 (MQ=255)
                                                                                cAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTgg                                                                      >  1:1865950/1‑150 (MQ=255)
                                                                                    ggCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGa                                                                  <  2:148168/150‑1 (MQ=255)
                                                                                       tGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCt                                                               >  2:431657/1‑150 (MQ=255)
                                                                                       tGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCt                                                               >  2:1321971/1‑150 (MQ=255)
                                                                                          cTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGt                                                            >  2:3839417/1‑150 (MQ=255)
                                                                                           ttCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtc                                                           >  1:1558817/1‑150 (MQ=255)
                                                                                            tCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtct                                                          >  1:3704759/1‑150 (MQ=255)
                                                                                             ccGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCACCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGGTCAGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtctc                                                         <  2:5433111/150‑1 (MQ=255)
                                                                                             ccGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtctc                                                         <  2:2396168/150‑1 (MQ=255)
                                                                                             ccGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGtctc                                                         <  2:2949397/150‑1 (MQ=255)
                                                                                                     tCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCa                                                 <  2:2332264/150‑1 (MQ=255)
                                                                                                      ccTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGc                                                                                                                   >  2:2939104/1‑83 (MQ=255)
                                                                                                       cTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATc                                               <  2:5173723/150‑1 (MQ=255)
                                                                                                              gCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGttt                                        <  1:5257437/150‑1 (MQ=255)
                                                                                                               cATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTCGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGtttt                                       <  1:3393375/150‑1 (MQ=255)
                                                                                                                aTGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTc                                      <  1:4420709/150‑1 (MQ=255)
                                                                                                                aTGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTc                                      >  1:2627640/1‑150 (MQ=255)
                                                                                                                 tGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCg                                     <  1:3374044/150‑1 (MQ=255)
                                                                                                                 tGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCg                                     >  1:2742726/1‑150 (MQ=255)
                                                                                                                  ggACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGc                                    <  1:5058845/150‑1 (MQ=255)
                                                                                                                     cacGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAg                                                                 >  1:2325167/1‑118 (MQ=255)
                                                                                                                     cacGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAg                                                                 <  2:2325167/118‑1 (MQ=255)
                                                                                                                       cGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                               <  2:5213012/150‑1 (MQ=255)
                                                                                                                       cGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTgg                               <  2:2704143/150‑1 (MQ=255)
                                                                                                                        ggTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTggg                              <  2:2456500/150‑1 (MQ=255)
                                                                                                                               gTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTc                       >  1:4498215/1‑150 (MQ=255)
                                                                                                                                   cAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAAc                   >  2:753423/1‑150 (MQ=255)
                                                                                                                                      tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                >  2:1005147/1‑150 (MQ=255)
                                                                                                                                      tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                <  1:1500045/150‑1 (MQ=255)
                                                                                                                                      tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                <  2:2762753/150‑1 (MQ=255)
                                                                                                                                      tCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGa                <  2:2763316/150‑1 (MQ=255)
                                                                                                                                           acacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCt           >  2:2833749/4‑150 (MQ=255)
                                                                                                                                            cccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc          <  2:18770/150‑1 (MQ=255)
                                                                                                                                            cacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc          >  2:3857622/3‑150 (MQ=255)
                                                                                                                                            cacAACAGCCTGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTc          >  2:2133183/3‑150 (MQ=255)
                                                                                                                                             ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCggtag                                                                                              >  2:1762873/1‑65 (MQ=255)
                                                                                                                                             ccAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTcc         >  1:913743/1‑150 (MQ=255)
                                                                                                                                              cAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGc                                                      <  2:2539816/104‑1 (MQ=255)
                                                                                                                                              cAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGc                                                      >  1:2539816/1‑104 (MQ=255)
                                                                                                                                              cAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCa        >  2:1496742/1‑150 (MQ=255)
                                                                                                                                               aaCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCAc       >  1:315067/1‑150 (MQ=255)
                                                                                                                                               aaCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCAc       >  1:3694078/1‑150 (MQ=255)
                                                                                                                                                aCAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACt      <  2:689863/150‑1 (MQ=255)
                                                                                                                                                    cagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTc  <  2:457044/150‑1 (MQ=255)
                                                                                                                                                    cagGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTc  <  2:4356258/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATTCTTGTCGTGGCGCAGGCCAAGCTCGGTGTAAAAGTATTCGTCATGAAAGCGATGCTGGACGAACAGCGCATGATTCCAGCGGCTGTCTTCCGTGAAGTCCTCACTCGCATGGACACGGTCTTCGTACCAGTCGCTACCCAACAGCAGGCTGTTGTGTTCACCCAGGGCTAGGTTGTTCTGCCAAGTTGCTTGGTCGCGGTAGGTATTGAAAGGTGCGCTTGAGTGGCTGAGCTTGTCTCGGCTTTCATCGCGGTTTTCGCTGTGGGCCAGCTCAAACCGAGTGCTCCACTGCTC  >  NC_002947/609633‑609930

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: