Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,795,625 A→T 5.0% intergenic (‑24/+33) PP_5073 ← / ← hemE hypothetical protein/uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,795,6250AT 5.0% 1311.1 / 32.7 390intergenic (‑24/+33)PP_5073/hemEhypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (190/179);  new base T (5/16);  total (195/195)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 397 total aligned reads displayed.

GCCAGCCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG  >  NC_002947/5795476‑5795774
                                                                                                                                                     |                                                                                                                                                     
gccagccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCa                                                                                                                                                       >  1:3710502/1‑150 (MQ=255)
 ccagccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAt                                                                                                                                                      >  2:1138429/1‑150 (MQ=255)
  cagccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATc                                                                                                                                                     >  2:8781409/1‑150 (MQ=255)
  cagccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATc                                                                                                                                                     >  1:6630240/1‑150 (MQ=255)
     ccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACTGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGTGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAATACCCGGCATCGaa                                                                                                                                                  >  1:7665438/1‑150 (MQ=255)
     ccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGaa                                                                                                                                                  <  2:8685591/150‑1 (MQ=255)
     ccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGaa                                                                                                                                                  <  2:7546017/150‑1 (MQ=255)
       gCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGc                                                                                                                                                >  2:5710944/1‑150 (MQ=255)
          aaTGCGCCAGTGCTGAACTAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                             <  2:431418/150‑1 (MQ=255)
          aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                             <  2:3225086/150‑1 (MQ=255)
          aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                             <  2:3737084/150‑1 (MQ=255)
           aTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCggg                                                                                                                                            <  2:4492586/150‑1 (MQ=255)
            tGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGt                                                                                                                                           >  2:6324437/1‑150 (MQ=255)
            tGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGt                                                                                                                                           <  2:2962399/150‑1 (MQ=255)
             gcgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtt                                                                                                                                          <  1:7560923/150‑1 (MQ=255)
             gcgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtt                                                                                                                                          <  1:1950291/150‑1 (MQ=255)
             gcgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGtt                                                                                                                                          <  1:4778989/150‑1 (MQ=255)
             gcgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGtt                                                                                                                                          <  2:3528693/150‑1 (MQ=255)
               gcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttt                                                                                                                                        <  1:313139/150‑2 (MQ=255)
               gcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTa                                                                                                                                        <  1:7828394/150‑1 (MQ=255)
               gcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTa                                                                                                                                        <  2:7414898/150‑1 (MQ=255)
                ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGttttt                                                                                                                                       <  1:2664166/150‑3 (MQ=255)
                ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTAc                                                                                                                                       <  2:4073788/150‑1 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttttt                                                                                                                                      <  2:2770645/150‑4 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttttt                                                                                                                                      <  1:6957399/150‑4 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttttt                                                                                                                                      <  1:8299003/150‑4 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                      <  1:7208296/150‑1 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                      >  2:4146488/1‑150 (MQ=255)
                 cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTACGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGattttt                                                                                                                                      <  2:5408321/150‑4 (MQ=255)
                    tGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                   >  1:7428433/1‑150 (MQ=255)
                    tGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                   <  1:1044825/150‑1 (MQ=255)
                      cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgca                                                                                                                                 <  2:6915088/150‑1 (MQ=255)
                      cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgca                                                                                                                                 <  1:8544615/150‑1 (MQ=255)
                       tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcag                                                                                                                                >  2:1342904/1‑150 (MQ=255)
                       tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcag                                                                                                                                >  1:5253256/1‑150 (MQ=255)
                        gAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATTGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagt                                                                                                                               <  2:3624235/150‑1 (MQ=255)
                          aCAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgc                                                                                                                             >  2:5089233/1‑150 (MQ=255)
                           cAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgca                                                                                                                            <  1:565688/150‑1 (MQ=255)
                           cAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgca                                                                                                                            <  2:1868286/150‑1 (MQ=255)
                            aaGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                           >  1:2358190/1‑150 (MQ=255)
                                    gAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCAc                                                                                                                   >  2:853867/1‑150 (MQ=255)
                                            gggCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTa                                                                                                           >  2:3121962/1‑150 (MQ=255)
                                              gCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACt                                                                                                         <  1:306349/150‑1 (MQ=255)
                                                  gaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAACAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcgc                                                                                                     <  2:1485188/150‑1 (MQ=255)
                                                  gaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcgc                                                                                                     <  1:2300245/150‑1 (MQ=255)
                                                   aTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGcc                                                                                                    <  1:5864829/150‑1 (MQ=255)
                                                     aaCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGa                                                                                                  >  1:4562930/1‑150 (MQ=255)
                                                     aaCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGa                                                                                                  <  1:5434616/150‑1 (MQ=255)
                                                      aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAATACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                 >  2:6707271/1‑150 (MQ=255)
                                                      aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                 >  1:2946281/1‑150 (MQ=255)
                                                      aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                 >  1:1402229/1‑150 (MQ=255)
                                                      aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                 >  1:8067009/1‑150 (MQ=255)
                                                      aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                 >  1:2760632/1‑150 (MQ=255)
                                                        ggATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAg                                                                                               >  2:7749688/1‑150 (MQ=255)
                                                        ggATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTGAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAg                                                                                               >  1:604494/1‑150 (MQ=255)
                                                         gATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGt                                                                                              <  2:2143690/150‑1 (MQ=255)
                                                         gATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGt                                                                                              >  2:2620615/1‑150 (MQ=255)
                                                                aTGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAg                                                                                                                       >  1:6682728/1‑118 (MQ=255)
                                                                       tCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt                                                                                >  1:3289511/1‑150 (MQ=255)
                                                                          ggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGTTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgat                                                                             >  2:8206065/1‑150 (MQ=255)
                                                                          ggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgat                                                                             <  1:3922068/150‑1 (MQ=255)
                                                                          ggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgat                                                                             <  1:5533306/150‑1 (MQ=255)
                                                                              cTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                         <  2:8322579/150‑1 (MQ=255)
                                                                               tGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGTGCCGACAGTTCGTGCACGGCGTTGATGAATg                                                                        >  2:2319711/1‑150 (MQ=255)
                                                                               tGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATg                                                                        >  1:3310624/1‑150 (MQ=255)
                                                                                 tttcggtttcggTGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCa                                                                      >  1:1245744/1‑150 (MQ=255)
                                                                                    cggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                                   >  1:5032263/1‑150 (MQ=255)
                                                                                    cggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                                   >  1:3553906/1‑150 (MQ=255)
                                                                                       tttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCa                                                                >  2:394311/1‑150 (MQ=255)
                                                                                          cgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGt                                                             <  2:8075576/150‑1 (MQ=255)
                                                                                            aTGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTc                                                           >  1:4089397/1‑150 (MQ=255)
                                                                                            aTGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTc                                                           >  1:2602964/1‑150 (MQ=255)
                                                                                               ggACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCCAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCgg                                                        <  2:6637014/150‑1 (MQ=255)
                                                                                               ggACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCgg                                                        <  1:8228259/150‑1 (MQ=255)
                                                                                               ggACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCgg                                                        <  2:1520774/150‑1 (MQ=255)
                                                                                                gACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCggg                                                       <  1:4360622/150‑1 (MQ=255)
                                                                                                 aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGt                                                      >  2:1581894/1‑150 (MQ=255)
                                                                                                 aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGt                                                      >  2:532049/1‑150 (MQ=255)
                                                                                                 aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGt                                                      <  2:6263432/150‑1 (MQ=255)
                                                                                                  ccAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTc                                                     <  2:8486806/150‑1 (MQ=255)
                                                                                                          gCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCg                                             >  1:6241163/1‑150 (MQ=255)
                                                                                                            gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGc                                           >  1:4431907/1‑150 (MQ=255)
                                                                                                            gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAACCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGc                                           <  1:1615145/150‑1 (MQ=255)
                                                                                                               gaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTg                                        >  1:4462162/1‑150 (MQ=255)
                                                                                                               gaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTg                                        >  2:3433511/1‑150 (MQ=255)
                                                                                                                    cTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                   >  2:4561281/1‑150 (MQ=255)
                                                                                                                     tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCg                                  <  1:7153942/150‑1 (MQ=255)
                                                                                                                     tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCg                                  <  2:4392427/150‑1 (MQ=255)
                                                                                                                      gAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGt                                 <  1:2192014/150‑1 (MQ=255)
                                                                                                                       aaTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                                <  1:6381478/150‑1 (MQ=255)
                                                                                                                            aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa                           <  2:4254970/150‑1 (MQ=255)
                                                                                                                               cTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGt                        >  1:5768592/1‑150 (MQ=255)
                                                                                                                                      tACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAacac                 <  1:6293957/150‑1 (MQ=255)
                                                                                                                                       aCAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAAcaca                >  2:7828064/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATg              >  1:1843547/1‑150 (MQ=255)
                                                                                                                                              aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAg         <  2:1776453/150‑1 (MQ=255)
                                                                                                                                              aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAg         <  1:4908724/150‑1 (MQ=255)
                                                                                                                                                 cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc      >  2:7178268/1‑150 (MQ=255)
                                                                                                                                                 cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc      >  1:1451493/1‑150 (MQ=255)
                                                                                                                                                 cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc      >  1:8366325/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGCCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG  >  NC_002947/5795476‑5795774

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: