Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,203 C→T 6.9% 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 6.9% 327.9 / 9.4 114intergenic (+75/‑112)hemN/fnrAoxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (52/54);  new base T (3/5);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 117 total aligned reads displayed.

TCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGC  >  NC_002947/4850060‑4850352
                                                                                                                                               |                                                                                                                                                       
tCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTg                                                                                                                                                   >  1:661264/1‑150 (MQ=255)
      tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                             <  2:1874750/150‑2 (MQ=255)
          tGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                         >  2:351964/1‑150 (MQ=255)
          tGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                         >  1:58407/1‑150 (MQ=255)
              gtTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCa                                                                                                                                     <  2:1920496/150‑1 (MQ=255)
              gtTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCa                                                                                                                                     <  1:439540/150‑1 (MQ=255)
                     gcgTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCa                                                                                                                            >  2:629677/1‑150 (MQ=255)
                              gCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATcg                                                                                                                   >  2:1411692/1‑150 (MQ=255)
                               cGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATcgc                                                                                                                  <  1:266885/150‑1 (MQ=255)
                                    cacaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGcaca                                                                                                             <  2:1546607/150‑1 (MQ=255)
                                      caACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACt                                                                                                           >  2:2233712/1‑150 (MQ=255)
                                         ccGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        <  1:351964/150‑1 (MQ=255)
                                         ccGCCAGCGTTATTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        <  1:2176365/150‑1 (MQ=255)
                                           gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAGGGCGTTCAACCGCATGGACAGGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCt                                                                                                      >  2:1952222/1‑150 (MQ=255)
                                            ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATAAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgc                                                                                                                                   >  1:411500/1‑122 (MQ=255)
                                                          gAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTg                                                                                       >  2:1265396/1‑150 (MQ=255)
                                                           aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                             >  1:2160888/1‑96 (MQ=255)
                                                           aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                             <  2:2160888/97‑2 (MQ=255)
                                                            gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGGGAg                                                                                     >  1:1108718/1‑150 (MQ=255)
                                                            gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTCGCCTACACCCTGTGTCGTGAg                                                                                     <  2:58407/150‑1 (MQ=255)
                                                             tGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGt                                                                                    >  2:660856/1‑150 (MQ=255)
                                                              gATCTGAAAAACAATGCGTTCAACCGCATGTACAAGCTTCCTTGTCAGGCACACTATGACCACACAATGCCTGACGGCAGGCCACCTGCCGGGGCTTCGCCCGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCTTGAGtt                                                                                   >  2:982070/1‑150 (MQ=255)
                                                               aTCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        >  1:1454022/1‑128 (MQ=255)
                                                               aTCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        <  2:1454022/128‑1 (MQ=255)
                                                                tCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGcc                                                                                                       <  2:617534/128‑1 (MQ=255)
                                                                tCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGcc                                                                                                       >  1:617534/1‑128 (MQ=255)
                                                                 cTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAcc                                                                                <  1:1545987/150‑1 (MQ=255)
                                                                  tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAccc                                                                               >  1:300450/1‑150 (MQ=255)
                                                                  tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAccc                                                                               <  2:661264/150‑1 (MQ=255)
                                                                  tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTAGCAACCAGGCGCGTGAGCAT‑‑CAAACACGGTCCGACACCCTGTGTCGTGAGTTAccc                                                                               >  2:1834225/1‑150 (MQ=255)
                                                                  tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAACACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAccc                                                                               >  2:1629353/1‑150 (MQ=255)
                                                                    aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGtt                                                                                                                                           >  2:1624884/1‑90 (MQ=255)
                                                                    aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGtt                                                                                                                                           <  1:1624884/90‑1 (MQ=255)
                                                                    aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        >  1:2091785/1‑123 (MQ=255)
                                                                      aaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTAc                                                                           <  2:1556928/150‑1 (MQ=255)
                                                                        aCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                         <  1:1469375/150‑1 (MQ=255)
                                                                          aaGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACt                                                                       >  2:2340627/1‑150 (MQ=255)
                                                                           aGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACtt                                                                      <  1:290529/150‑1 (MQ=255)
                                                                               gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCTTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                                  <  1:1952222/150‑1 (MQ=255)
                                                                               gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                                  >  1:1750940/1‑150 (MQ=255)
                                                                                ttCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtg                                                                 <  2:300450/150‑1 (MQ=255)
                                                                                 tCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgt                                                                <  1:1265396/150‑1 (MQ=255)
                                                                                    aCCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAAC‑‑CCGCACCCACATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCt                                                             >  1:1250469/1‑150 (MQ=255)
                                                                                         aTGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      >  1:1380102/1‑120 (MQ=255)
                                                                                         aTGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      <  2:1380102/120‑1 (MQ=255)
                                                                                            gACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcaca                                                     >  2:172562/1‑150 (MQ=255)
                                                                                             aCAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACaa                                                    >  2:753311/1‑150 (MQ=255)
                                                                                              cAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAg                                                   <  2:1750940/150‑1 (MQ=255)
                                                                                               aaGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      >  2:2911/1‑114 (MQ=255)
                                                                                                aGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTCCCCCTACGCCTTATGTGTGCTTTCCCACAAGGa                                                 <  1:2233712/150‑1 (MQ=255)
                                                                                                 gCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGc                                                                                                        >  2:1292514/1‑94 (MQ=255)
                                                                                                  cTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATc                                               >  1:2140445/1‑150 (MQ=255)
                                                                                                    tCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGa                                             >  1:314727/1‑150 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGNCACCTGNCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTTTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                          >  1:1766392/1‑150 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                          <  2:2140445/150‑1 (MQ=255)
                                                                                                       ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGGCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                          >  1:752009/1‑150 (MQ=255)
                                                                                                         gccAGCCACACTATGAGCACACAAGGACTCCCGGCAAGCCAACTGCCGGGGTTTCGCCAGC‑‑ACGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTCCGACTTATGTGAGCTTTCACACAAGGATCGATTCaa                                        >  2:129452/3‑150 (MQ=255)
                                                                                                         gTCAGGCACACTATGAGCACACAAGGCCTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCaa                                        <  2:1098501/150‑1 (MQ=255)
                                                                                                          tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCtt                                                            >  2:1971366/1‑129 (MQ=255)
                                                                                                          tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCtt                                                            <  1:1971366/129‑1 (MQ=255)
                                                                                                              gCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      <  1:1797153/99‑1 (MQ=255)
                                                                                                              gCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      >  2:1797153/1‑99 (MQ=255)
                                                                                                               cacaCTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgtg                                                               <  2:395200/121‑1 (MQ=255)
                                                                                                               cacaCTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgtg                                                               >  1:395200/1‑121 (MQ=255)
                                                                                                                acacTATGAGCACACAAGGACTCCCGGCAGGCCAACTACCGCGGCTTCGTTAGC‑‑TCGCACGAGCATCGCTCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCg                                 >  1:657964/1‑150 (MQ=255)
                                                                                                                 cacTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                                >  1:777249/1‑150 (MQ=255)
                                                                                                                  acTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAg                               >  2:1029762/1‑150 (MQ=255)
                                                                                                                  acTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTCTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGt                                                                                        <  2:541185/93‑1 (MQ=255)
                                                                                                                    tATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCaaa                                       <  1:128839/140‑1 (MQ=255)
                                                                                                                    tATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCaaa                                       >  2:128839/1‑140 (MQ=255)
                                                                                                                     aTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACCCTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCa                            <  2:1704928/150‑1 (MQ=255)
                                                                                                                       gAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGcc                                                                                                       >  1:1133026/1‑73 (MQ=255)
                                                                                                                       gAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGcc                                                                                                       <  2:1133026/73‑1 (MQ=255)
                                                                                                                       gAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTc                                                                                          <  1:218703/86‑1 (MQ=255)
                                                                                                                       gAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTc                                                                                          >  2:218703/1‑86 (MQ=255)
                                                                                                                         gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTc                         >  1:2257480/1‑149 (MQ=255)
                                                                                                                         gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTc                         <  2:2257480/149‑1 (MQ=255)
                                                                                                                           acacaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCaaa                      >  1:289920/1‑150 (MQ=255)
                                                                                                                              caAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                  >  1:2059275/1‑135 (MQ=255)
                                                                                                                              caAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                  <  2:2059275/135‑1 (MQ=255)
                                                                                                                               ccgatcTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCTAAAAAATTTCGTACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgc                  <  2:1269679/145‑1 (MQ=255)
                                                                                                                               aaGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgc                  <  1:1411692/150‑1 (MQ=255)
                                                                                                                                    accccGGCAGGTGGCCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                         >  1:247708/3‑90 (MQ=255)
                                                                                                                                    accccGGCAGGTGGCCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                         <  2:247708/88‑1 (MQ=255)
                                                                                                                                     tCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcac            >  1:1333995/1‑150 (MQ=255)
                                                                                                                                     ccccGGCAGGTGGCCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCttt                                                           >  1:2352768/2‑103 (MQ=255)
                                                                                                                                      cccGGCAGGTGGCCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCa                            <  1:187255/133‑1 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgtg                                                               >  2:2107437/1‑97 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCt                                                             <  1:905999/99‑1 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCt                                                             >  2:905999/1‑99 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                                >  2:1331392/1‑128 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGCCGGGGTTTAGGCTTC‑‑TCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgcgagatcgg          >  1:1269679/1‑143 (MQ=255)
                                                                                                                                       ccGGCAGGCCACCTGACGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                                <  1:1331392/128‑1 (MQ=255)
                                                                                                                                          gCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAAc       <  1:2048339/150‑1 (MQ=255)
                                                                                                                                              gCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATc                                               <  1:1557688/106‑1 (MQ=255)
                                                                                                                                              gCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATc                                               >  2:1557688/1‑106 (MQ=255)
                                                                                                                                              gCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACaa        <  2:995042/145‑1 (MQ=255)
                                                                                                                                              gCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACaa        >  1:995042/1‑145 (MQ=255)
                                                                                                                                              cccACCTGCCGGCGTTTCGCCACC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAgg   <  2:752009/149‑1 (MQ=255)
                                                                                                                                               ccACCTGCCGGGGTTTCGCCAGC‑‑GCGCCCCATCACCGCGCCCTCGGTGCCACACCCTGTGTCGTGAGTTAACCTCACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGc  <  1:2340627/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
TCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGC‑‑GCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGC  >  NC_002947/4850060‑4850352

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: