Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,603,669 T→C 6.3% intergenic (+1000/‑502) PP_4074 → / → PP_4076 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,603,6690TC 6.3% 186.5 / 2.1 63intergenic (+1000/‑502)PP_4074/PP_4076hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/2);  total (30/33)
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.00e-01

TCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAAGA  >  NC_002947/4603524‑4603812
                                                                                                                                                 |                                                                                                                                               
tCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATccc                                                                                                                                              >  2:67212/1‑149 (MQ=255)
    ccTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTAcc                                                                                                                                          <  4:203430/149‑1 (MQ=255)
        ccTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCTCTACCGCTGGTGCCCCAAAATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAg                                                                                                                                       <  5:344296/148‑1 (MQ=255)
        ccTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGa                                                                                                                                        <  1:67212/147‑1 (MQ=255)
          tGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTa                                                                                                                                    >  7:77967/1‑149 (MQ=255)
            ttGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACt                                                                                                                                  >  4:49465/1‑149 (MQ=255)
               gCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCt                                                                                                                                             >  4:123932/1‑135 (MQ=255)
               gCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCt                                                                                                                                             <  3:123932/135‑1 (MQ=255)
                  aGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTc                                                                                                                            >  1:48425/1‑149 (MQ=255)
                         ttAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCg                                                                                                                      >  5:323219/1‑148 (MQ=255)
                         ttAACGACAGCTTTGCCCGCTACCGATGGTGACCCAACATCATCGGCCCGAAATCGGGAGGTGCAGTAGAAGGCCTGCCGGTGCATCTGTAAGAATCGCTGGGCCAGCTGCATGCAAAGATCCCTACCGACGTACTGATTTCCGACCg                                                                                                                      >  5:163787/1‑148 (MQ=255)
                          tAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGc                                                                                                                     >  1:360544/1‑148 (MQ=255)
                              gACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                  >  1:193099/1‑147 (MQ=255)
                              gACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAg                                                                                                                  >  3:169712/1‑147 (MQ=255)
                               aCAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAg                                                                                                               <  2:48425/149‑1 (MQ=255)
                                cAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGt                                                                                                              <  8:96302/149‑1 (MQ=255)
                                  gCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTc                                                                                                            <  6:323219/149‑1 (MQ=255)
                                   cTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTc                                                                                                            >  6:365824/1‑148 (MQ=255)
                                    tttGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTc                                                                                                            >  7:13165/1‑147 (MQ=255)
                                       gCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGc                                                                                                        <  8:221520/148‑1 (MQ=255)
                                          cGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGc                                                                                                    >  5:372454/1‑149 (MQ=255)
                                                 gCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAagag                                                                                             >  4:289784/1‑149 (MQ=255)
                                                  cTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAgg                                                                                            <  2:360544/149‑1 (MQ=255)
                                                   tGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAggg                                                                                           >  2:67869/1‑149 (MQ=255)
                                                    ggTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGACGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCAGACAGCAAGGACTACCAGCTGGCCGAAGAggtc                                                                                          >  8:260428/1‑147 (MQ=255)
                                                    ggTGCCCCAAAAGCATCGGCCCGCATGCGGGAGGTGCCGTCGATGGCCTGGCGGTGCATCTGTTCGTCTCGGTGGGCCAGCTGGCTGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGGCGAAGATGGc                                                                                          <  8:155588/149‑1 (MQ=255)
                                                      tGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCtt                                                                                        <  5:365824/149‑1 (MQ=255)
                                                        ccccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTa                                                                                      <  2:193099/149‑1 (MQ=255)
                                                          ccAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATc                                                                                    <  1:227397/149‑1 (MQ=255)
                                                           cAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCATACCGAGGTACTGATTTCCGACCGCAAGGACTTCGAGCTGGCCGAAGAGGGCTTTATCg                                                                                   >  6:285953/1‑149 (MQ=255)
                                                              catcatCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGccc                                                                                >  4:337763/1‑149 (MQ=255)
                                                                tcatcGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCt                                                                               >  1:7654/1‑148 (MQ=255)
                                                                  atcGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGa                                                                             <  4:169712/148‑1 (MQ=255)
                                                                      gCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATg                                                                        <  8:77967/149‑1 (MQ=255)
                                                                      gCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATg                                                                        <  2:7654/149‑1 (MQ=255)
                                                                      gCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGACCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAag                                                                                               >  4:277785/1‑126 (MQ=255)
                                                                      gCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGACCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAag                                                                                               <  3:277785/126‑1 (MQ=255)
                                                                                          gTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGACCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAAcgc                                                    >  2:204250/1‑149 (MQ=255)
                                                                                          gTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGACCCCTACCGAGGTACTGATTTCAGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAAcgc                                                    >  6:178178/1‑149 (MQ=255)
                                                                                           tCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAAcgcc                                                   >  1:165309/1‑149 (MQ=255)
                                                                                              aaGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCcgcg                                                <  2:165309/149‑1 (MQ=255)
                                                                                                      ccGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTAGTGATTTCCGACCGCTAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTtctc                                        <  2:289958/149‑1 (MQ=255)
                                                                                                      ccGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTtcttct                                         <  3:49465/148‑1 (MQ=255)
                                                                                                                 gTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCTGCGTTCTTCTCAGCCAACTCgg                             >  2:122334/1‑149 (MQ=255)
                                                                                                                  tACGAATCGCTGGTCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCAGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACt                                <  7:260428/145‑1 (MQ=255)
                                                                                                                  tACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAATGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGt                            <  3:289784/149‑1 (MQ=255)
                                                                                                                  tACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGt                            >  7:160987/1‑149 (MQ=255)
                                                                                                                     gAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAgg                                                                                            <  4:346957/82‑1 (MQ=255)
                                                                                                                     gAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAgg                                                                                            >  3:346957/1‑82 (MQ=255)
                                                                                                                                ccAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGc                    >  6:175548/1‑143 (MQ=255)
                                                                                                                                ccAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGc                    <  5:175548/143‑1 (MQ=255)
                                                                                                                                  aGCTGCAGGCCAAGACCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTGAGCCAACTCGGTGCAGAAGCCGAAAAAc            <  1:204250/149‑1 (MQ=255)
                                                                                                                                  aGCTGCAGGCCAAGACCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAAc            <  5:178178/149‑1 (MQ=255)
                                                                                                                                    cTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGttc                                             <  8:59073/114‑1 (MQ=255)
                                                                                                                                    cTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGttc                                             >  7:59073/1‑114 (MQ=255)
                                                                                                                                       cAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTc         <  3:337763/147‑1 (MQ=255)
                                                                                                                                          gCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCaa    <  8:160987/149‑1 (MQ=255)
                                                                                                                                          gCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCaa    <  6:372454/149‑1 (MQ=255)
                                                                                                                                            cAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAAGa  <  1:67869/149‑1 (MQ=255)
                                                                                                                                            cAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAAGa  <  8:13165/149‑1 (MQ=255)
                                                                                                                                            cAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAAGa  >  3:308116/1‑149 (MQ=255)
                                                                                                                                                aTCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTa                                                                                      <  8:151379/61‑1 (MQ=255)
                                                                                                                                                aTCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTa                                                                                      >  7:151379/1‑61 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCAGCCTTCCTGTTGGCCAGCTGCATTAACGACAGCTTTGCCCGCTACCGCTGGTGCCCCAACATCATCGGCCCGCAATCGGGAGGTGCCGTCGAAGGCCTGCCGGTGCATCTGTACGAATCGCTGGGCCAGCTGCAGGCCAAGATCCCTACCGAGGTACTGATTTCCGACCGCAAGGAGTTCGAGCTGGCCGAAGAGGGCTTTATCGCCCTGACCATGCGCAAGGACAGCGACAACGCCGCGTTCTTCTCAGCCAACTCGGTGCAGAAGCCGAAAAACTTCCCCAAGA  >  NC_002947/4603524‑4603812

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