Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,205 A→G 7.1% intergenic (+77/‑110) 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,2050AG 7.1% 200.3 / 3.9 70intergenic (+77/‑110)hemN/fnrAoxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base G (2/3);  total (34/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 = 7.66e-01

TCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGCC  >  NC_002947/4850060‑4850353
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                        >  2:95295/1‑150 (MQ=255)
          tGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                        >  1:139813/1‑150 (MQ=255)
           ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGttaga                                                                                                                                       >  1:874871/1‑147 (MQ=255)
           ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                           <  2:874955/146‑1 (MQ=255)
           ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGt                                                                                                                                           >  1:874955/1‑146 (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:594694/150‑1 (MQ=255)
                                                                               gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                                   >  2:1127753/1‑150 (MQ=255)
                                                                                     ccGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            >  1:1133578/1‑70 (MQ=38)
                                                                                     ccGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  2:1133578/71‑2 (MQ=38)
                                                                                                aGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGc                                                                                                         <  1:831580/95‑1 (MQ=255)
                                                                                                aGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGc                                                                                                         >  2:831580/1‑95 (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                                                                                                                   >  1:1840457/1‑78 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgc                                                                                                                   <  2:1840457/78‑1 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                           <  2:1576273/150‑1 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                           >  1:885189/1‑150 (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)
                                                                                                                                                cACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGAATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGcc  >  1:315259/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGCC  >  NC_002947/4850060‑4850353

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: