Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,496 G→C 9.3% A116A (GCG→GCC)  rep → single‑stranded DNA dependent ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,011,4960GC 9.3% 188.8 / 7.8 64A116A (GCG→GCC) repsingle‑stranded DNA dependent ATPase
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (4/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TCGCCACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCG  >  NC_002947/6011351‑6011638
                                                                                                                                                 |                                                                                                                                              
tCGCCACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTg                                                                                                                                             <  5:135249/149‑1 (MQ=255)
tCGCCACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTg                                                                                                                                             <  5:205292/149‑1 (MQ=255)
      ccTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCg                                                                                                                                       <  2:44579/149‑1 (MQ=255)
      ccTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCg                                                                                                                                       <  1:72393/149‑1 (MQ=255)
          gctgcGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAc                                                                                                                                    <  7:65181/148‑1 (MQ=255)
               gcCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATg                                                                                                                              <  4:222130/149‑1 (MQ=255)
               gcCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATg                                                                                                                              <  3:67120/149‑1 (MQ=255)
                    ggggAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGaa                                                                                                                         <  5:1825/149‑1 (MQ=255)
                     ggggAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAg                                                                                                                        <  1:6691/149‑1 (MQ=255)
                     gggAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAg                                                                                                                        >  3:12228/1‑149 (MQ=255)
                      ggAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAgg                                                                                                                       <  6:84958/149‑1 (MQ=255)
                       gAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAgg                                                                                                                       >  7:21334/1‑148 (MQ=255)
                        aaGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGaa                                                                                                                     <  5:252844/149‑1 (MQ=255)
                                  ccTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCg                                                                                                              >  3:351490/1‑146 (MQ=255)
                                  ccTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCg                                                                                                              <  4:351490/146‑1 (MQ=255)
                                   cTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGcga                                                                                                          >  3:77579/1‑149 (MQ=255)
                                             gtACGTTCCACAACCTGGGCCTGAACATAATCCGCAAGGAACACGAGCGCCTGGGCAATAAGACGGGTTTCTCGATCGTCGACGAGTCCGACATCAAGGCGATGCTGTCGGACATCATGCAGAAGGAATACTCAGGCGACGACGGCATc                                                                                                >  1:294171/1‑149 (MQ=255)
                                                 gTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgac                                                                                             >  5:167844/1‑148 (MQ=255)
                                                  ttCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                           <  8:21334/149‑1 (MQ=255)
                                                  ttCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                           <  7:311519/149‑1 (MQ=255)
                                                  ttCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                           <  3:285071/149‑1 (MQ=255)
                                                   tCCACCGCCGGGGCCTGACCGTCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTGTCTCGATCTTCGACGAGTCCGGCATCAAGGCGGGGCTGTCGGGCATCATGCAGAAGGAGTGCTCCGGCGCGGCCGGCGTCGACgag                                                                                          >  2:350345/1‑149 (MQ=255)
                                                    ccACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgac                                                                                             >  3:42539/1‑145 (MQ=255)
                                                     cacaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAcatcat                                                                                                                               <  2:241940/110‑1 (MQ=255)
                                                     cacaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAcatcat                                                                                                                               >  1:241940/1‑110 (MQ=255)
                                                      acaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGAt                                                                                        >  3:105785/1‑148 (MQ=255)
                                                                    aaCATCATCCGCAAGGAACACGATCGCCTCGGCTATAAGCCGTGTTTCTCGATCTTCGACGAGTCCGAGATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGTCTACGACGGCATCGAAGAGATCAAGAACATTATc                                                                           <  6:89593/147‑1 (MQ=255)
                                                                    aaCATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCgg                                                                         >  1:326276/1‑149 (MQ=255)
                                                                    aaCATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCg                                                                          <  4:12228/148‑1 (MQ=255)
                                                                         catcCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTAGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACCGCATCGACGAGATCAAGAAAATGATCGGTGCCt                                                                    >  5:324724/1‑149 (MQ=255)
                                                                             cGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCACGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTgg                                                                  >  8:249668/1‑147 (MQ=255)
                                                                             cGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGgaa                                                                <  4:42539/149‑1 (MQ=255)
                                                                              gCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTg                                                                   >  2:342966/1‑145 (MQ=255)
                                                                              gCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTg                                                                   <  1:342966/145‑1 (MQ=255)
                                                                                  ggAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAAc                                                            >  6:400/1‑148 (MQ=255)
                                                                                       acGAGCGCCTGGGCTAAAAGCCGGGATTCTCGATCTTCGACTAGTCCGACATCAAGGCGTTGCTGTCGGACATCATCCAGAAGGAATACTCCGGCGACGACGGCAACGACGAGATCAAGAACAGGATAGGTGCCTCGAAGAACGCCCTg                                                      >  4:54314/1‑149 (MQ=255)
                                                                                         gAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGt                                                    >  7:294391/1‑149 (MQ=255)
                                                                                          aGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTc                                                   <  2:326276/149‑1 (MQ=255)
                                                                                                          gCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCAACGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGccc                                   >  3:156495/1‑149 (MQ=255)
                                                                                                           ccGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCt                                  >  3:254399/1‑149 (MQ=255)
                                                                                                              ggTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGa                               <  4:77579/149‑1 (MQ=255)
                                                                                                               gtctctcGATCATCGACCAGTCCGCCCTCAACCCCCTGCAGGCGGACATCATGCAGAAGGAATACTCCCCCGACCACGGCACAGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGaa                              <  3:282222/146‑1 (MQ=255)
                                                                                                               ggttCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTg                                                                   <  3:87350/110‑1 (MQ=255)
                                                                                                               gTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGa                               >  2:96525/1‑148 (MQ=255)
                                                                                                                  tctcTATCTTCGACGAGTCAGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACTGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGaaaa                            <  6:324724/148‑1 (MQ=255)
                                                                                                                  tctcGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCa                                          >  8:283099/1‑134 (MQ=255)
                                                                                                                  tctcGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCa                                          <  7:283099/134‑1 (MQ=255)
                                                                                                                   ctcGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATTCTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAg                           <  7:249668/148‑1 (MQ=255)
                                                                                                                   ctcGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAgg                          >  2:285294/1‑149 (MQ=255)
                                                                                                                      gATCTTCGACGAGTCCGACATCACGGCGTTGCTGTCGGACATCATGCAGAATGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGgc                         <  1:350345/147‑1 (MQ=255)
                                                                                                                      gATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGGGCCTGGAAGAACGACCTGGTca                                                  >  2:358686/1‑121 (MQ=255)
                                                                                                                             gacgaGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAAcc                 <  6:167844/148‑1 (MQ=255)
                                                                                                                             atgatgTCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTTCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAAccc                >  2:99142/6‑149 (MQ=255)
                                                                                                                                gaGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCcgcg             <  4:254399/149‑1 (MQ=255)
                                                                                                                                 aggccGACCTCAAGGCGGTGCCGGCGGACATCATGCAGAAGGAATACTCCGCCGACGACCGCACCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCcgcgc            <  3:316812/146‑1 (MQ=255)
                                                                                                                                  gTCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGgcgc                       >  1:360072/1‑137 (MQ=255)
                                                                                                                                  gTCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAAAACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGgcgc                       <  2:360072/137‑1 (MQ=255)
                                                                                                                                   tCCGCCAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGAACGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCcgcgc            >  4:77462/1‑147 (MQ=255)
                                                                                                                                   tCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCcgcgc            >  6:211623/1‑147 (MQ=255)
                                                                                                                                      gACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAg         >  2:357013/1‑147 (MQ=255)
                                                                                                                                      gACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAg         >  6:315672/1‑147 (MQ=255)
                                                                                                                                           cAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGAccg  >  5:302250/1‑149 (MQ=255)
                                                                                                                                            aaGGCGTTGCTGTCGGACATCATGCAGAAGGATTACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGAc                                                         <  2:139758/93‑1 (MQ=255)
                                                                                                                                            aaGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGAc                                                         <  2:145149/93‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCGCCACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCG  >  NC_002947/6011351‑6011638

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: ‑