Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,294,122 A→G 7.4% intergenic (+17/+23) pdxH → / ← PP_1130 pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,294,1220AG 7.4% 262.0 / 6.5 91intergenic (+17/+23)pdxH/PP_1130pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/46);  new base G (2/5);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCTG  >  NC_002947/1293973‑1294269
                                                                                                                                                     |                                                                                                                                                   
tGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGCCCTTTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCaa                                                                                                                                                     >  1:1779546/1‑150 (MQ=255)
   gCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCg                                                                                                                                                  <  1:1299592/150‑1 (MQ=255)
    cATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGCCGGGCAGTGGCTGCGAGAGCGGCTGGCGCGCTGAGCAACTGAGGCCGCCAAGCgg                                                                                                                                                 >  2:275159/1‑150 (MQ=255)
     aTTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGc                                                                                                                                                >  1:517184/1‑150 (MQ=255)
           gatcttACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgt                                                                                                                                          <  1:2250156/145‑2 (MQ=255)
             ggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCa                                                                                                                                        <  2:2155999/150‑1 (MQ=255)
               aTACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAtt                                                                                                                                      <  1:1476905/150‑1 (MQ=255)
               aTACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAtt                                                                                                                                      <  1:837338/150‑1 (MQ=255)
                tACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCCCCCTGGGCCGCCAAGCGGCCCCAGCCATTg                                                                                                                                     <  1:62163/150‑1 (MQ=255)
                 aCCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGGAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGcccc                                                                                                                                             <  1:1023244/141‑1 (MQ=255)
                 aCCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCACCTGGGGCCGAATAGCGGATCCAGCGAtata                                                                                                                                    >  2:2250156/1‑147 (MQ=255)
                 aCCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGcccc                                                                                                                                             >  2:1023244/1‑141 (MQ=255)
                 aCCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGa                                                                                                                                    >  1:2154980/1‑150 (MQ=255)
                   cGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATc                                                                                                                                  <  1:979691/150‑1 (MQ=255)
                   cGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATc                                                                                                                                  <  2:2154980/150‑1 (MQ=255)
                     tttGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                                <  1:708062/150‑1 (MQ=255)
                               gAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCgg                                                                                                                      >  2:180615/1‑150 (MQ=255)
                                 gcgcATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGGTCAAGCTACCCCGGAt                                                                                                                    >  2:691019/1‑150 (MQ=255)
                                     attgAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                                >  2:973150/4‑150 (MQ=255)
                                     aTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                                <  2:1671940/150‑1 (MQ=255)
                                     aTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCCGAGTGTCTTGGTCCGCCAAGCGGCACCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                                <  1:1227852/150‑1 (MQ=255)
                                      tCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCt                                                                                                               >  2:197689/1‑150 (MQ=255)
                                              tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGAACTGGCTGCGAGCGAGAGTGGCGCGCTGAGCAGCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTc                                                                                                       >  1:1451487/1‑150 (MQ=255)
                                              tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGCGGAACGGCAGTGGCTGCGGGTGCGGCTGGCGCCCTGGGCAACGGGCGCGGCGAAGCGGCCCCAGCCATTGATCAAGCTATCCCGGATAAGGTgccgccgc                                                                                                       <  2:325133/150‑3 (MQ=255)
                                                  aGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCgcgc                                                                                                   >  1:549027/1‑150 (MQ=255)
                                                    gggcgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgttg                                                                                                                                        <  1:373244/111‑4 (MQ=255)
                                                        gggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAac                                                                                             <  2:549027/150‑1 (MQ=255)
                                                         ggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAacc                                                                                            >  1:2002223/1‑150 (MQ=255)
                                                             aGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAGCTGGGGCCGCTTGGCGGCCCCAgttgct                                                                                                                                      >  1:1601742/1‑99 (MQ=38)
                                                             aGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAGCTGGGGCCGCTTGGCGGCCCCAgttgct                                                                                                                                      <  2:1601742/104‑6 (MQ=37)
                                                             aGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAccc                                                                                        <  1:1643210/150‑1 (MQ=255)
                                                              gtcgtcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCCGCCATTGATCAAGCTACCCCTGATAACCTGCCGCCTCGCGCTCCAACCAcccc                                                                                       <  2:1518796/150‑1 (MQ=255)
                                                                cgtcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGc                                                                                     <  2:2169093/150‑1 (MQ=255)
                                                                  tcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCa                                                                                   >  1:1563722/1‑150 (MQ=255)
                                                                    tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGt                                                                                 >  2:1184964/1‑150 (MQ=255)
                                                                     gCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                   >  2:1598184/1‑83 (MQ=255)
                                                                     gCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                   <  1:1598184/83‑1 (MQ=255)
                                                                     gCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                                <  2:2152778/150‑1 (MQ=255)
                                                                      cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCa                                                                                                                                 >  2:729449/1‑100 (MQ=255)
                                                                      cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCa                                                                                                                                 <  1:729449/100‑1 (MQ=255)
                                                                       aTGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGc                                                                              >  2:610182/1‑150 (MQ=255)
                                                                       aTGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCCCCAACCACCCCGCCAGTTGc                                                                              >  1:444841/1‑150 (MQ=255)
                                                                       aTGATCGTTTGAACTATCGGCTGGTGGACGGACAGTGTCTGCGAGAGCGGCTGGCTCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGATACCCCGGATAACCTGCCTGCTCTCGCTGCAACCACCCCGCCAGTTGc                                                                              <  2:1779546/150‑1 (MQ=255)
                                                                             gTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGc                                                                        <  2:1451487/150‑1 (MQ=255)
                                                                              tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCaa                                                                                                                                                     <  2:1990212/72‑1 (MQ=255)
                                                                              tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCaa                                                                                                                                                     >  1:1990212/1‑72 (MQ=255)
                                                                              tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAt                                                                                                                                       <  1:393420/86‑1 (MQ=255)
                                                                              tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAt                                                                                                                                       >  2:393420/1‑86 (MQ=255)
                                                                               ttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGccc                                                                      <  1:180615/150‑1 (MQ=255)
                                                                                  aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGTTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAc                                                                   <  1:2239535/150‑1 (MQ=255)
                                                                                  aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAc                                                                   >  2:1409119/1‑150 (MQ=255)
                                                                                   aCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAcc                                                                  <  1:841544/150‑1 (MQ=255)
                                                                                      aTCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtc                                                               >  1:421890/1‑150 (MQ=255)
                                                                                         ggCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                            <  2:517184/150‑1 (MQ=255)
                                                                                               tggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCcgc                                                      <  1:197689/150‑1 (MQ=255)
                                                                                                    gggCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgtt                                                                                                                                         >  2:496726/1‑60 (MQ=25)
                                                                                                    gggCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgtt                                                                                                                                         <  1:496726/62‑3 (MQ=25)
                                                                                                     gccAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAGTCTTCCTTAGCGTTTGATCAAGCACCACAGTATAACCTGCCGCCTCCCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgc                                                <  1:524243/148‑1 (MQ=255)
                                                                                                      gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                                <  2:1239625/117‑1 (MQ=255)
                                                                                                      gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                                >  1:1239625/1‑117 (MQ=255)
                                                                                                          gggCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCaa                                           <  2:1990491/149‑1 (MQ=255)
                                                                                                             cTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCt                                        >  2:717642/1‑150 (MQ=255)
                                                                                                             cTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCt                                        >  2:1794042/1‑150 (MQ=255)
                                                                                                              tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                             <  2:1174165/80‑1 (MQ=255)
                                                                                                              tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                             >  1:1174165/1‑80 (MQ=255)
                                                                                                              tGCGAGAGCGGCTGGCGCCCTAAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGGTCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCtt                                       >  2:1389142/1‑150 (MQ=255)
                                                                                                               gCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCttt                                      <  1:378303/150‑1 (MQ=255)
                                                                                                                 gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTcc                                    <  1:1147090/150‑1 (MQ=255)
                                                                                                                 gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTGTcc                                    <  1:1184964/150‑1 (MQ=255)
                                                                                                                   gagCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCaa                                           <  1:1126313/141‑1 (MQ=255)
                                                                                                                   gagCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCaa                                           >  2:1126313/1‑141 (MQ=255)
                                                                                                                      cGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCAt                               <  1:1794042/150‑1 (MQ=255)
                                                                                                                         cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGaa                            >  2:1924116/1‑150 (MQ=255)
                                                                                                                         cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGaa                            <  2:421890/150‑1 (MQ=255)
                                                                                                                         cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGaa                            <  1:410258/150‑1 (MQ=255)
                                                                                                                         cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGaa                            <  1:973150/150‑1 (MQ=255)
                                                                                                                            gcgcCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgc                         >  2:1825151/1‑150 (MQ=255)
                                                                                                                              gcCCTGAGAAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgc                       >  1:1898879/1‑150 (MQ=255)
                                                                                                                               cccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgcg                      >  1:949795/1‑150 (MQ=255)
                                                                                                                                ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGt                     >  2:2009936/1‑150 (MQ=255)
                                                                                                                                      cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa               <  1:610182/150‑1 (MQ=255)
                                                                                                                                       aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgc                         <  2:1339297/139‑1 (MQ=255)
                                                                                                                                       aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg              <  1:717642/150‑1 (MQ=255)
                                                                                                                                       aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg              >  1:1358623/1‑150 (MQ=255)
                                                                                                                                       aaCGGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgc                         >  1:1339297/1‑139 (MQ=255)
                                                                                                                                             ggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCga        <  2:1563722/150‑2 (MQ=255)
                                                                                                                                              gccgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtc                                                               >  1:1396259/1‑94 (MQ=255)
                                                                                                                                              gccgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtc                                                               <  2:1396259/94‑1 (MQ=255)
                                                                                                                                              gccgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGc       >  2:556297/1‑150 (MQ=255)
                                                                                                                                                 gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGcccc    >  2:647759/1‑150 (MQ=255)
                                                                                                                                                  ccAAGCGTCCTCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   <  2:492568/150‑1 (MQ=255)
                                                                                                                                                  ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCATGATCGGCCCCt   <  1:1825151/150‑1 (MQ=255)
                                                                                                                                                  ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   >  1:1395872/1‑150 (MQ=255)
                                                                                                                                                   cAAGCGGCCCAAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCACGATCGGCCCCTg  <  1:1389142/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCTG  >  NC_002947/1293973‑1294269

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