Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,203 C→T 7.2% intergenic (+75/‑112) hemN → / → fnrA oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,850,2030CT 7.2% 200.4 / 4.4 69intergenic (+75/‑112)hemN/fnrAoxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGCGCTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGG  >  NC_002947/4850054‑4850351
                                                                                                                                                     |                                                                                                                                                    
gcgcgcTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGc                                                                                                                                                      <  1:897621/150‑1 (MQ=255)
gcgcgcTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGc                                                                                                                                                      <  2:226169/150‑1 (MQ=255)
      tCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTg                                                                                                                                                >  2:327606/1‑150 (MQ=255)
       cGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGc                                                                                                                                               <  2:29854/150‑1 (MQ=255)
           aTGCATGGTGTTCGCCGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCggg                                                                                                                                           <  2:199907/150‑1 (MQ=255)
            tGCGTGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTAGGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                          <  2:76465/150‑2 (MQ=255)
             tcttGGTGTTAAACGCATATCTGGAGTTGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                         <  2:874871/147‑1 (MQ=255)
             gCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACCCTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                         <  1:1448981/150‑1 (MQ=255)
              cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGtt                                                                                                                                        >  1:1904917/1‑150 (MQ=255)
                tGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                      >  1:139813/1‑150 (MQ=255)
                tGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                      >  2:95295/1‑150 (MQ=255)
                 ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGttaga                                                                                                                                     >  1:874871/1‑147 (MQ=255)
                 ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                         >  1:874955/1‑146 (MQ=255)
                 ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                         <  2:874955/146‑1 (MQ=255)
                   tgtTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGcc                                                                                                                                   <  1:201234/150‑1 (MQ=255)
                    gtTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCa                                                                                                                                  <  2:69101/150‑1 (MQ=255)
                         aCGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgcgcg                                                                                                                             >  1:14705/1‑150 (MQ=255)
                                  tGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcatcat                                                                                                                    >  1:402774/1‑150 (MQ=255)
                                  tGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCATCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcatcat                                                                                                                    >  2:1895709/1‑150 (MQ=255)
                                   ggCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCatcatc                                                                                                                   <  1:525613/150‑1 (MQ=255)
                                        tGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGc                                                                                                                                    <  1:969478/128‑1 (MQ=255)
                                        tGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGc                                                                                                                                    >  2:969478/1‑128 (MQ=255)
                                                 gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCt                                                                                                     >  1:1576273/1‑150 (MQ=255)
                                                  ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcatc                                                                                                                      >  2:1222687/1‑132 (MQ=255)
                                                  ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcatc                                                                                                                      <  1:1222687/132‑1 (MQ=255)
                                                  ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTa                                                                                                    >  1:533073/1‑150 (MQ=255)
                                                    aGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTaca                                                                                                  <  2:755887/150‑1 (MQ=255)
                                                     gCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGc                                                                                                                                    >  1:752653/1‑115 (MQ=255)
                                                     gCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGc                                                                                                                                    <  2:752653/115‑1 (MQ=255)
                                                       gTTTCTAGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCTGCGCGCACAATCATCGCGCACACTGGCCTACAccc                                                                                               >  2:1340536/1‑150 (MQ=255)
                                                                 aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                     <  1:1044905/150‑1 (MQ=255)
                                                                 aGTGAGCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                     <  2:453900/150‑1 (MQ=255)
                                                                       cTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAcc                                                                               <  1:186266/150‑1 (MQ=255)
                                                                        tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCATCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAccc                                                                              <  1:95295/150‑1 (MQ=255)
                                                                            aaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTAGGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTAc                                                                          <  2:162707/150‑1 (MQ=255)
                                                                               cAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGAc                                                                       <  1:265015/150‑1 (MQ=255)
                                                                                  ggCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTa                                                                    <  2:402774/150‑1 (MQ=255)
                                                                                   gCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAt                                                                   <  2:533073/150‑1 (MQ=255)
                                                                                     gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                                 >  2:1127753/1‑150 (MQ=255)
                                                                                     gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                                 <  2:594694/150‑1 (MQ=255)
                                                                                           ccGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                          >  1:1133578/1‑70 (MQ=38)
                                                                                           ccGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                          <  2:1133578/71‑2 (MQ=38)
                                                                                                      aGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGc                                                                                                       >  2:831580/1‑95 (MQ=255)
                                                                                                      aGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGc                                                                                                       <  1:831580/95‑1 (MQ=255)
                                                                                                       gCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTcc                                                        <  1:1743195/141‑1 (MQ=255)
                                                                                                       gCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTcc                                                        >  2:1743195/1‑141 (MQ=255)
                                                                                                       gCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGAt                                               <  2:14705/150‑1 (MQ=255)
                                                                                                          tCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGa                                            <  1:1895709/150‑1 (MQ=255)
                                                                                                            cTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCACGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGAtt                                          <  1:327606/150‑1 (MQ=255)
                                                                                                             ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgc                                                                                                                 <  2:1840457/78‑1 (MQ=255)
                                                                                                             ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgc                                                                                                                 >  1:1840457/1‑78 (MQ=255)
                                                                                                             ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                         >  1:885189/1‑150 (MQ=255)
                                                                                                             ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                         <  2:1576273/150‑1 (MQ=255)
                                                                                                               gTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTa                                                                                 >  1:177329/1‑108 (MQ=255)
                                                                                                               gTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTa                                                                                 <  2:177329/108‑1 (MQ=255)
                                                                                                                tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCaaa                                      >  1:1363635/1‑150 (MQ=255)
                                                                                                                 cAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAAt                                     <  2:1632858/150‑1 (MQ=255)
                                                                                                                    gCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTc                                  >  2:184165/1‑150 (MQ=255)
                                                                                                                    gCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTc                                  >  1:1885366/1‑150 (MQ=255)
                                                                                                                     cacaCTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                 >  1:1283899/1‑150 (MQ=255)
                                                                                                                       cacTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                               <  2:1363635/150‑1 (MQ=255)
                                                                                                                           aTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCa                           <  2:1904917/150‑1 (MQ=255)
                                                                                                                                 acacaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACACGGCTCGCTTCCAATGTCTGACCCAGTCaat                     >  2:1543022/1‑149 (MQ=255)
                                                                                                                                   acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCTACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACt                   >  1:1843209/1‑150 (MQ=255)
                                                                                                                                     aaGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTAAAACTgc                 >  2:1046901/1‑150 (MQ=255)
                                                                                                                                          cTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                 >  2:573213/1‑129 (MQ=255)
                                                                                                                                          cTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                 <  1:573213/129‑1 (MQ=255)
                                                                                                                                           ccccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgt                                                               <  2:1857808/97‑1 (MQ=255)
                                                                                                                                             ccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcacac         >  1:1436771/1‑150 (MQ=255)
                                                                                                                                               ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGAACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACaa       >  2:1071193/1‑150 (MQ=255)
                                                                                                                                                    gCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAgg  >  2:1365045/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCGCTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGG  >  NC_002947/4850054‑4850351

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