Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,138,413 T→C 6.5% intergenic (‑228/‑89) PP_3641 ← / → PP_3643 cytosine/purine/uracil/thiamine/allantoin permease family protein/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,138,4130TC 6.5% 180.3 / 2.6 61intergenic (‑228/‑89)PP_3641/PP_3643cytosine/purine/uracil/thiamine/allantoin permease family protein/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base C (2/2);  total (27/34)
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

GGTATTCCAAACTTTTGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGT  >  NC_002947/4138265‑4138560
                                                                                                                                                      |                                                                                                                                                   
ggTATTCCAAACTTTTGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTAGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAAtt                                                                                                                                                     >  3:201817/1‑149 (MQ=255)
 gTATTCCAAACTTTTGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTc                                                                                                                                                    <  3:217022/149‑1 (MQ=255)
       cAAACTTTTGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCtt                                                                                                                                               <  6:328536/148‑1 (MQ=255)
            ttttGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCaaaa                                                                                                                                         >  4:253394/1‑149 (MQ=255)
            ttttGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATCCCGCTTGCaaaa                                                                                                                                         >  3:250651/1‑149 (MQ=255)
            ttttGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTAAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCACGGTTGAAAATCCCTCTTGCaaaa                                                                                                                                         >  3:246391/1‑149 (MQ=255)
                gACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAgg                                                                                                                                     >  4:134235/1‑149 (MQ=255)
                       cTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACg                                                                                                                              <  2:123821/149‑1 (MQ=255)
                       cTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACg                                                                                                                              <  3:6209/149‑1 (MQ=255)
                          aTCCGCATGAAAAAACGGCCACAGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACgg                                                                                                                             <  2:136803/147‑1 (MQ=255)
                              gCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATAc                                                                                                                        >  7:172127/1‑148 (MQ=255)
                                aTGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCAt                                                                                                                     >  4:136884/1‑149 (MQ=255)
                                 ttaaaaaaCGGCAACTTTCCGGCTTGAAAAACCCTTCCTAAAAACC‑CCTTAACGGTCAATGGCAAGGTTTTTGGGCTTATTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCAt                                                                                                                     <  1:261334/147‑1 (MQ=255)
                                 tGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCAt                                                                                                                     >  6:214991/1‑148 (MQ=255)
                                 tGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATc                                                                                                                    <  7:303691/149‑1 (MQ=255)
                                 tGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATc                                                                                                                    <  5:214991/149‑1 (MQ=255)
                                  gAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCAt                                                                                                                     >  4:150188/1‑147 (MQ=255)
                                   aaaaaaCGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCAt                                                                                                                     >  1:151790/1‑146 (MQ=255)
                                   aaaaaaCGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAg                                                                                                                  <  2:80099/149‑1 (MQ=255)
                                             cACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAAc                                                                                                        <  3:134235/149‑1 (MQ=255)
                                                     cGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATaaaaa                                                                                                >  3:25455/1‑149 (MQ=255)
                                                      gCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAAc                                                                                               >  2:6719/1‑149 (MQ=255)
                                                           cAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCTCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACcgcg                                                                                           <  3:258363/148‑1 (MQ=255)
                                                           cAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACcgcg                                                                                           <  3:136884/148‑1 (MQ=255)
                                                            aaaaCCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGa                                                                                         >  2:160800/1‑149 (MQ=255)
                                                            aaaaCCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGa                                                                                         >  5:147439/1‑149 (MQ=255)
                                                             aaaCCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGac                                                                                        <  2:151790/149‑1 (MQ=255)
                                                                       gACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGcc                                                                              >  1:607/1‑149 (MQ=255)
                                                                        aCAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGt                                                                                                                                    >  8:150064/1‑94 (MQ=255)
                                                                        aCAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGt                                                                                                                                    <  7:150064/94‑1 (MQ=255)
                                                                         cAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCt                                                                            <  4:25455/149‑1 (MQ=255)
                                                                              cGCCTCACCTCTCAA‑TGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATCCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCaga                                                                       <  4:246391/149‑1 (MQ=255)
                                                                              cGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATCCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCaga                                                                       <  4:250651/149‑1 (MQ=255)
                                                                                   cACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGa                                                                                                                 <  6:301060/102‑1 (MQ=255)
                                                                                   cACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGa                                                                                                                 >  5:301060/1‑102 (MQ=255)
                                                                                    aCCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTACAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACt                                                                 >  6:216770/1‑149 (MQ=255)
                                                                                       gCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTccc                                                              <  8:172127/149‑1 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCaga                                                                       <  3:240636/130‑1 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCaga                                                                       >  4:240636/1‑130 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGAc                                                    <  3:253394/149‑1 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGAc                                                    >  8:127736/1‑149 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCCCAGCCCATCAGAGGTACTCCCGAAGATCgc                                                     >  6:66180/1‑147 (MQ=255)
                                                                                                  aCGTTTTT‑CGCGGCTTTCGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCTGACCACCGCAGCCCTTCAGAGGTACTCCCg                                                             <  5:66180/140‑1 (MQ=255)
                                                                                                         t‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTg                                             >  8:174692/1‑149 (MQ=255)
                                                                                                              ggCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTc                                                                          >  8:208849/1‑116 (MQ=255)
                                                                                                              ggCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTc                                                                          <  7:208849/116‑1 (MQ=255)
                                                                                                               gCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTACAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTAc                                        <  5:216770/149‑1 (MQ=255)
                                                                                                                    gggCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGa                                    <  7:174692/148‑1 (MQ=255)
                                                                                                                            cATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTg                            <  6:147439/148‑1 (MQ=255)
                                                                                                                            cATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTg                            <  7:127736/148‑1 (MQ=255)
                                                                                                                             aTCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTggg                          <  1:6719/149‑1 (MQ=255)
                                                                                                                                gcgAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACa                                       <  1:104846/133‑1 (MQ=255)
                                                                                                                                gcgAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACa                                       >  2:104846/1‑133 (MQ=255)
                                                                                                                                    aaaCGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTc                   >  5:168876/1‑149 (MQ=255)
                                                                                                                                      aCGCAGGGTTGAAAATTCCGCTTGCAAAATATGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTccc                 <  4:201817/149‑1 (MQ=255)
                                                                                                                                         cAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTc              <  6:17796/149‑1 (MQ=255)
                                                                                                                                             gTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTccc                 >  7:156625/1‑142 (MQ=255)
                                                                                                                                             gTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTccc                 <  8:156625/142‑1 (MQ=255)
                                                                                                                                             gTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGc          <  2:607/149‑1 (MQ=255)
                                                                                                                                              ttGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCg         >  8:85589/1‑149 (MQ=255)
                                                                                                                                                 aaaaTTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTa      >  2:264285/1‑149 (MQ=255)
                                                                                                                                                    aTTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTTTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCg   <  7:85589/149‑1 (MQ=255)
                                                                                                                                                     ttCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGt  <  3:150188/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGTATTCCAAACTTTTGACAACACTGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAA‑GGCCACGTTTTT‑CGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGT  >  NC_002947/4138265‑4138560

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