Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,787,137 G→T 5.6% intergenic (+104/‑78) tsf → / → pyrH elongation factor Ts/uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,787,1370GT 5.6% 229.7 / 3.1 71intergenic (+104/‑78)tsf/pyrHelongation factor Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAGCGGCGAGGCC  >  NC_002947/1786992‑1787283
                                                                                                                                                 |                                                                                                                                                  
gctgctGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCgg                                                                                                                                                 >  2:119968/1‑149 (MQ=255)
gctgctGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCg                                                                                                                                                  <  5:187835/148‑1 (MQ=255)
gctgctGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCACAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCg                                                                                                                                                  <  4:143132/148‑1 (MQ=255)
gctgctGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCACAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCg                                                                                                                                                  <  4:138325/148‑1 (MQ=255)
    ctgAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGc                                                                                                                                              <  1:119968/148‑1 (MQ=255)
      gAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCAc                                                                                                                                           >  3:245419/1‑149 (MQ=255)
       aaGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGATGAGTTTATACGACTCCCTCTTCGCTGGCACCGCTTCGGTGCCag                                                                                                                                           <  6:220388/148‑2 (MQ=255)
            gCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGc                                                                                                                                         >  7:106902/1‑145 (MQ=255)
            gCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGc                                                                                                                                         <  8:106902/145‑1 (MQ=255)
             ccGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACa                                                                                                                                     <  6:99230/148‑1 (MQ=255)
             ccGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACa                                                                                                                                     <  8:168989/148‑1 (MQ=255)
             ccGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAg                                                                                                                                    >  2:233899/1‑149 (MQ=255)
                   tgggcGCTGCTGCCAAGCAGAAAGACAGCCCCGTCTGACGCCCCAACGAGGCTGCCCGCTCTCGCGCGCAGCGTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCTCCGAAGCGGTGGCACGCTACAGTTAGCa                                                                                                                              <  4:236555/145‑1 (MQ=255)
                   aGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAg                                                                                                                                <  5:26252/147‑1 (MQ=255)
                        gctgctGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCt                                                                                                                          <  1:233899/148‑1 (MQ=255)
                             tgcCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCaaa                                                                                                                     >  1:309411/1‑148 (MQ=255)
                                cAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCTCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTTGCACCGAAGCGGTGCCACGCTACATTTAGCAGGCTGAAAACAGc                                                                                                                 >  4:121214/1‑149 (MQ=255)
                                 aaGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGTGGCACCGCTTCGGTGCCagcgaa                                                                                                                                       >  4:68400/1‑120 (MQ=255)
                                 aaGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGTGGCACCGCTTCGGTGCCagcgaa                                                                                                                                       <  3:68400/125‑6 (MQ=255)
                                  aGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGccc                                                                                                               <  2:50657/149‑1 (MQ=255)
                                             aGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCt                                                                                                                          <  8:114259/127‑1 (MQ=255)
                                               ccccGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTc                                                                                                   >  5:273066/1‑148 (MQ=255)
                                                        tCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGt                                                                                            <  4:245419/146‑1 (MQ=255)
                                                          gCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgc                                                                                        <  4:214156/148‑1 (MQ=255)
                                                           ccccAAAGAGGCTGCCCGCTCACGCGTGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                                                                                      >  4:235739/1‑149 (MQ=255)
                                                           ccccAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                                                                                      >  1:125110/1‑149 (MQ=255)
                                                            cccAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                                                                                      >  1:174404/1‑148 (MQ=255)
                                                             ccAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCa                                                                                    <  1:208794/149‑1 (MQ=255)
                                                               aaaGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAg                                                                                                                  <  7:36411/117‑1 (MQ=38)
                                                               aaaGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAg                                                                                                                  >  8:36411/1‑117 (MQ=38)
                                                               aaaGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAg                                                                                   >  7:182778/1‑148 (MQ=255)
                                                               aaaGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAg                                                                                   >  7:50989/1‑148 (MQ=255)
                                                                   agGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgagag                                                                              <  2:125110/149‑1 (MQ=255)
                                                                    ggCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTg                                                                                                                         >  7:193916/1‑105 (MQ=255)
                                                                    ggCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTg                                                                                                                         <  8:193916/105‑1 (MQ=255)
                                                                         cccGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCa                                                                                                                                            <  6:165580/81‑1 (MQ=255)
                                                                         cccGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCa                                                                                                                                            >  5:165580/1‑81 (MQ=255)
                                                                            gcgtgagcgggcAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAAt                                                                     <  4:7999/143‑1 (MQ=255)
                                                                                        aGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGt                                                         <  8:76556/149‑1 (MQ=255)
                                                                                        aGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGt                                                         >  7:22782/1‑149 (MQ=255)
                                                                                         gCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGt                                                         >  5:248992/1‑148 (MQ=255)
                                                                                          ccTCTTTGTCAAAACGGGGAGTAGTTTATACGACTCCCTCTTCGCTTTCACCTAAGCGGTGCCACGATACAGTTAGCAGGCTGCAAACAGCCCGCACTAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGa                                                       <  3:121214/149‑1 (MQ=255)
                                                                                          ccTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAACCAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGt                                                         >  6:101739/1‑147 (MQ=255)
                                                                                              tttGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAAt                                                                     <  1:39641/131‑1 (MQ=255)
                                                                                              tttGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAAt                                                                     >  2:39641/1‑131 (MQ=255)
                                                                                                tGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcg                                                                                         <  5:53009/109‑1 (MQ=255)
                                                                                                tGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcg                                                                                         >  6:53009/1‑109 (MQ=255)
                                                                                                  tCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGc                                               <  8:69266/149‑1 (MQ=255)
                                                                                                    aaaaCGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGcc                                              <  8:228739/148‑1 (MQ=255)
                                                                                                    aaaaCGGGGAGGAGTTTATACGACTCCCTCTTCGCTCGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGcc                                              >  6:127092/1‑148 (MQ=255)
                                                                                                      aaCGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg                                                                         >  1:169082/1‑119 (MQ=255)
                                                                                                      aaCGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg                                                                         <  2:169082/119‑1 (MQ=255)
                                                                                                           ggaggaGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGc                                      <  2:309411/149‑1 (MQ=255)
                                                                                                                    tataCGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACg                              >  2:21338/1‑148 (MQ=255)
                                                                                                                        cGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCAtt                          <  2:174404/148‑1 (MQ=255)
                                                                                                                             cccTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCt                     >  5:137340/1‑148 (MQ=255)
                                                                                                                                    cGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACt               >  2:89309/1‑147 (MQ=255)
                                                                                                                                     gCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTa             <  4:174289/148‑1 (MQ=255)
                                                                                                                                     gCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTa             <  6:248992/148‑1 (MQ=255)
                                                                                                                                      gtGGCACCGCTTCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAgc           >  5:120917/2‑149 (MQ=255)
                                                                                                                                         gCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAgcggc        >  2:81282/1‑149 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCtat                                   <  2:310822/120‑1 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCtat                                   >  1:310822/1‑120 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACtt              <  7:33120/141‑1 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACtt              >  8:33120/1‑141 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAgcggc        >  5:100569/1‑147 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAgcggc        >  5:323670/1‑147 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAgcggc        >  5:105076/1‑147 (MQ=255)
                                                                                                                                           aCCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAGCGGCGa      <  6:273066/149‑1 (MQ=255)
                                                                                                                                             cGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAGCGGCGAgg    <  8:50989/149‑1 (MQ=255)
                                                                                                                                               aaGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAGCGGCGAGGcc  <  1:81282/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGGCTCAGCAGGTGAGTGGTCGCCAACCTCGCTATAAACGCATTTTGCTCAAACTTAGCGGCGAGGCC  >  NC_002947/1786992‑1787283

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: