Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,586,136 G→C 5.5% intergenic (+42/‑500) PP_t39 → / → PP_2266 tRNA‑Leu/DNA‑directed RNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,586,1360GC 5.5% 235.7 / 3.6 74intergenic (+42/‑500)PP_t39/PP_2266tRNA‑Leu/DNA‑directed RNA polymerase
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base C (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GCGCCTTTGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCATCCACT  >  NC_002947/2585988‑2586281
                                                                                                                                                     |                                                                                                                                                  
gcgcCTTTGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACgg                                                                                                                                                    >  5:88764/1‑149 (MQ=255)
    cTTTGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTAGAGTCTCTTCGTCAGCACAACGAATTAGAAAAAATCTCCCCGGATTCATGAAAGCCACACGGCTc                                                                                                                                                 >  2:186322/1‑148 (MQ=255)
      ttGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTACAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAg                                                                                                                                             >  3:236383/1‑149 (MQ=255)
       tGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCaa                                                                                                                                              <  8:147766/147‑1 (MQ=255)
             acgCGTTACGCTGACGTGGGGAAATTGCTAGTCACTCTCGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCg                                                                                                                                       <  4:113015/146‑1 (MQ=37)
              aGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCgg                                                                                                                                      <  6:88764/148‑1 (MQ=255)
                agTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGa                                                                                                                                   <  1:131982/148‑1 (MQ=255)
                  ttACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGAtt                                                                                                                                 <  3:237771/149‑1 (MQ=255)
                        ggACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTc                                                                                                                           <  5:98013/149‑1 (MQ=255)
                            gtggtgAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGcc                                                                                                                        >  2:9769/1‑148 (MQ=255)
                            gtggtgAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCt                                                                                                                       <  8:274632/149‑1 (MQ=255)
                                 gAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCAc                                                                                                                  >  7:186271/1‑149 (MQ=255)
                                   aaTTTGTAGACACACTGTATTTAGGTTCCAGGGGCGCAAGCTGTAAGTGTTCGAGTCTCTTCGTCCGCTCGACGAATTAGAATAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACtt                                                                                                                <  3:253714/149‑1 (MQ=255)
                                   aaTTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCACAGGGTTCGGGGATTTTTTTCGCCTATCAc                                                                                                                  <  1:214714/148‑1 (MQ=255)
                                           gACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAAt                                                                                                        <  1:8479/149‑1 (MQ=255)
                                              acacTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCgg                                                                                                                                      >  1:325631/1‑116 (MQ=255)
                                              acacTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCgg                                                                                                                                      <  2:325631/116‑1 (MQ=255)
                                                acTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATa                                                                                                    >  5:45994/1‑148 (MQ=255)
                                                acTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATAc                                                                                                   >  6:96037/1‑149 (MQ=255)
                                                acTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATAc                                                                                                   >  2:305999/1‑149 (MQ=255)
                                                acTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCAGCACCACGAATTAGAAAAAATCTCCACGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGCGATTTTTTTCGACTATCACTTGAACGAATGATAc                                                                                                   >  6:314411/1‑149 (MQ=38)
                                                    gATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGaa                                                                                                         <  3:17231/139‑1 (MQ=255)
                                                    gATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGaa                                                                                                         >  4:17231/1‑139 (MQ=255)
                                                          ggTTCCAGCTGCGCAAGCTGTAAGAGTTCTAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGTATTCCTTAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCATATCACTTGCACGAATGATACCTTTCCACat                                                                                         <  1:186322/149‑1 (MQ=38)
                                                          ggTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATATCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTATTTCGCCTATCACTTGCACGAATGATACCTTTCcaca                                                                                          >  2:352402/1‑148 (MQ=255)
                                                           gtaccAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGCCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACAAATGATAACTTTCCACat                                                                                         >  1:54907/4‑148 (MQ=255)
                                                           gTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTACTTCGTCCGCACCACGAATTAGAAAAAATATCCCCGGATTCCTGAAAGCCACACGCCTCAAGGGTTCGGGGATTTTTTTCGCCTATCCCTTGCACAAATGATACCTTTCcaca                                                                                          >  1:235586/1‑148 (MQ=38)
                                                               cAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATaa                                                                                     <  5:257843/148‑1 (MQ=255)
                                                                aGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATc                                                                                   <  5:96037/149‑1 (MQ=255)
                                                                   ggcgCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGcc                                                                                <  1:325515/149‑1 (MQ=255)
                                                                      gcAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCTCGAATGATACCTTTCCACATATAATCGCCTtc                                                                             >  4:297020/1‑149 (MQ=255)
                                                                       cAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCCCACGGCTCAAGGGTTCGGGGAGTTTTTTCGCCTAACACTTTCACGAATGAAAACTTTCCACATAAAATCGCATtct                                                                            >  8:145092/1‑149 (MQ=38)
                                                                             gTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATc                                                                                   <  2:192295/136‑1 (MQ=255)
                                                                             gTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATc                                                                                   >  1:192295/1‑136 (MQ=255)
                                                                              tAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCaaa                                                                     <  3:160961/149‑1 (MQ=255)
                                                                              tAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCa                                                                       >  6:136030/1‑147 (MQ=255)
                                                                                agagTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAaga                                                                   <  5:136030/149‑1 (MQ=255)
                                                                                agagTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAaga                                                                   <  1:153752/149‑1 (MQ=255)
                                                                                   gTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAagagag                                                                >  4:93950/1‑149 (MQ=255)
                                                                                    ttCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGc                                                               <  1:9769/149‑1 (MQ=255)
                                                                                        aGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCaaaa                                                           <  5:76158/149‑1 (MQ=255)
                                                                                        aGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCaa                                                             >  3:109126/1‑147 (MQ=255)
                                                                                         gTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAg                                                          <  3:93950/149‑1 (MQ=255)
                                                                                          tctcTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGCCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACAAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAgg                                                         <  2:235586/149‑1 (MQ=255)
                                                                                          tctcTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGCCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACAAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAgg                                                         <  2:54907/149‑1 (MQ=255)
                                                                                            tctTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAggg                                                        >  5:369307/1‑148 (MQ=255)
                                                                                                  gTCCTCACCACGAATTAGAAAAAATCTCCCAGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGcc                                                                                <  2:181969/119‑1 (MQ=255)
                                                                                                  gTCCGCACCACGAATTAGAAAAAATCTACCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACTAATGATACCTTTACAAATATAATCGca                                                                                >  1:181969/1‑118 (MQ=255)
                                                                                                     cGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTg                                               <  6:74164/149‑1 (MQ=255)
                                                                                                      gttccacGAATTAGAAAAAATGTCCCCGTATTCCTGTAAGCCACACTGCTCAATGTTTCGTTGATTTTTTTCGCCTATCACTTGCACGATTGATTCCTTTCCTCATATAATCGCCTTCTCACCAAAGAGAGCAAAAGGGACACTTTTgg                                              <  5:314411/146‑1 (MQ=255)
                                                                                                      gCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAATCCACACGTCTCAATGGTTCGGTGATTTTTTTCGACTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAACGAGAGCAAAAGGGACACTTTTgg                                              <  1:222052/149‑1 (MQ=255)
                                                                                                      gCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCAACAAATAATCGCCTACTCAGCAAAGAGAGCAAAAGGGACACTTTTgg                                              >  6:92751/1‑149 (MQ=255)
                                                                                                         ccacGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGCGATTTTTTACCCCTATCACTGGCACCAATGATACCTTACCACATATAAACGCCTTCTCAGCAAAGACACAAAAAGGGACACTTTTGGAac                                           >  6:17563/1‑149 (MQ=255)
                                                                                                          cacGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGaca                                          <  8:186271/149‑1 (MQ=255)
                                                                                                          cacGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGaca                                          >  6:230220/1‑149 (MQ=255)
                                                                                                          cacGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGaca                                          >  4:131028/1‑149 (MQ=255)
                                                                                                           acGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTGCGGGGATTTTTTTCGCCTATCACTTGACCGAATGATACCTTTCCACATAGAATCGCCTTGTCAGCAAAGAGAGCATAAGGGACACTTGTGGGacc                                          >  4:216740/1‑147 (MQ=255)
                                                                                                                   gAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCAc                                                                                                                  >  3:361915/1‑68 (MQ=255)
                                                                                                                   gAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCAc                                                                                                                  <  4:361915/68‑1 (MQ=255)
                                                                                                                     aaaaaTCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCtt                               >  3:160066/1‑149 (MQ=255)
                                                                                                                        aaTCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAACGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCttttt                            <  3:216740/149‑1 (MQ=255)
                                                                                                                              cccGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGTGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGa                          <  6:369307/145‑1 (MQ=255)
                                                                                                                                cGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCtctt                    <  1:352402/149‑1 (MQ=255)
                                                                                                                                 ggATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCtcttc                   <  1:305999/149‑1 (MQ=255)
                                                                                                                                  gATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCtcttct                  <  6:45994/149‑1 (MQ=255)
                                                                                                                                  gATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCtcttct                  >  2:100152/1‑149 (MQ=255)
                                                                                                                                       cTGAAAGCCACACGGCTCAAGGGTTCGTGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGcc             <  7:71697/149‑1 (MQ=255)
                                                                                                                                       cTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGAc                         >  5:257501/1‑137 (MQ=255)
                                                                                                                                       cTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCGCTTCTGAGcc             >  8:171528/1‑149 (MQ=255)
                                                                                                                                        tGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGccc            >  7:84710/1‑149 (MQ=255)
                                                                                                                                          aaaGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGcccc           >  4:167367/1‑148 (MQ=255)
                                                                                                                                            aGCCACACGGCTCAAGGGTTCGGTGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTc         <  5:92751/148‑1 (MQ=255)
                                                                                                                                             gCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCAt       >  4:158830/1‑149 (MQ=255)
                                                                                                                                              ccACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCATc      <  5:230220/149‑1 (MQ=255)
                                                                                                                                                  accgCTCAATGGTTAGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCATCCACt  <  5:17563/146‑1 (MQ=255)
                                                                                                                                                  acGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCATCCACt  <  3:131028/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGCCTTTGTTTTAAGCGTTACGCGGACGTGGTGAAATTGGTAGACACACTGGATTTAGGTTCCAGCGGCGCAAGCTGTAAGAGTTCGAGTCTCTTCGTCCGCACCACGAATTAGAAAAAATCTCCCCGGATTCCTGAAAGCCACACGGCTCAAGGGTTCGGGGATTTTTTTCGCCTATCACTTGCACGAATGATACCTTTCCACATATAATCGCCTTCTCAGCAAAGAGAGCAAAAGGGACACTTTTGGGACACAAAGGCCCTTTTTGACCTCTTCTGAGCCCCTCATCCACT  >  NC_002947/2585988‑2586281

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