Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,069 G→A 7.9% intergenic (+30/‑468) PP_4308 → / → PP_4309 AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,900,0690GA 7.9% 182.0 / 4.8 64intergenic (+30/‑468)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/3);  total (31/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.76e-01

ATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAC  >  NC_002947/4899921‑4900216
                                                                                                                                                      |                                                                                                                                                    
aTGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCg                                                                                                                                                      <  5:112084/149‑1 (MQ=255)
     cgcTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGccc                                                                                                                                                 >  4:282731/1‑149 (MQ=255)
       cTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCgc                                                                                                                                                     >  8:36234/1‑142 (MQ=255)
       cTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCgc                                                                                                                                                     <  7:36234/143‑2 (MQ=255)
       cTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                               >  8:247744/1‑149 (MQ=255)
       cTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCCGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                               <  4:144244/149‑1 (MQ=255)
           aTGGAGCAGCTTGACACCCTTCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGTTTGATTGGCTGCTGT‑‑GCTGGCCTCTTCGCGGGCACGCGTGCGAag                                                                                                                                           <  4:161347/149‑1 (MQ=255)
            tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                          <  8:38500/149‑1 (MQ=255)
              gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑ACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                         <  3:352828/148‑1 (MQ=39)
                  agcTGGACACCCTGCTCCACCTCATTAATGTGCTGGCCGAGACAAAAACCCCTATCAACAACAAGACCCCGGTAAAGCGGCGGTTGCCGCCAATGGATTTTTTGGCTCCTGT‑‑GCTGGCCTCTTCGCGGGCACGCCTGCTAAGAGGCCAg                                                                                                                                    <  2:363514/149‑1 (MQ=25)
                        acacCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAgg                                                                                                                              <  3:247732/149‑1 (MQ=255)
                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                       >  4:360826/1‑136 (MQ=255)
                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                       <  3:360826/136‑1 (MQ=255)
                             cTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                                         <  4:48628/149‑1 (MQ=255)
                              tGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAg                                                                                                                                    <  1:146245/137‑1 (MQ=255)
                                    aCCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTa                                                                                                                 >  4:134539/1‑149 (MQ=255)
                                      cTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtt                                                                                                               >  4:188866/1‑149 (MQ=255)
                                       tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                    >  7:73077/1‑143 (MQ=255)
                                       tCATGAATGGGCTGGCCGAGACCAATACCGCGATACTCAATAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                    <  8:73077/143‑1 (MQ=255)
                                         aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAg                                                                                                                                    <  1:163594/126‑1 (MQ=255)
                                         aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAg                                                                                                                                    >  2:163594/1‑126 (MQ=255)
                                          tGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                           <  3:216826/149‑1 (MQ=255)
                                          tGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                           <  3:262680/149‑1 (MQ=255)
                                             aTGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                                        <  3:188866/149‑1 (MQ=255)
                                                tgCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTaca                                                                                                      >  7:217513/1‑148 (MQ=255)
                                                  cTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTAcaca                                                                                                    <  1:68459/148‑1 (MQ=255)
                                                   tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaa                                                                                                   >  3:240357/1‑148 (MQ=255)
                                                   tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCAAATACCAGGCCTGCCCCCGGAGAGGCCCGTACAAGAAAAGCATAGTATTTTTTTTAAACaa                                                                                                   >  2:184694/1‑148 (MQ=255)
                                                     gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaaa                                                                                                <  3:282731/149‑1 (MQ=255)
                                                      ccGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCTAAGAGGCCTGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAAc                                                                                               <  1:66142/149‑1 (MQ=255)
                                                       cGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAACGCTGAGTATTTTTTGTACACAAAACCg                                                                                              >  1:78175/1‑149 (MQ=255)
                                                          gaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAg                                                                                                                               >  8:11720/1‑114 (MQ=255)
                                                              aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCAt                                                                                       <  7:247744/149‑1 (MQ=255)
                                                               aTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATa                                                                                      <  1:184694/149‑1 (MQ=255)
                                                                 aCCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAgg                                                                                    >  8:125139/1‑149 (MQ=255)
                                                                  ccGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                                   <  6:300191/149‑1 (MQ=255)
                                                                   cgcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTg                                                                                  >  7:167659/1‑149 (MQ=255)
                                                                      gATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGAGATGCCATCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACTTAACCATAGGTg                                                                                  >  1:291769/1‑148 (MQ=255)
                                                                          atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAATGTAACCATAGGTGCCAtttt                                                                           >  6:132332/1‑149 (MQ=255)
                                                                             aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtg                                                                         <  4:240357/148‑1 (MQ=255)
                                                                               gaagaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                      <  3:134539/149‑1 (MQ=255)
                                                                               gaagaCCCCGGTCAAGCGGCGGGTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTTTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                      <  6:288470/149‑1 (MQ=255)
                                                                                         gTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttctt                                                            <  6:105221/149‑1 (MQ=255)
                                                                                                    gTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTACTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGccc                                                 <  7:241819/149‑1 (MQ=255)
                                                                                                          gccAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTc                                           >  5:51192/1‑149 (MQ=255)
                                                                                                             aaTGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                                         >  1:294409/1‑148 (MQ=255)
                                                                                                                 gattgattGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCGGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTacac                                    >  3:3259/1‑149 (MQ=255)
                                                                                                                       ttgGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTc                              <  2:291769/149‑1 (MQ=255)
                                                                                                                          gCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACATGAAAAGCTGCAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                          <  1:70646/103‑1 (MQ=255)
                                                                                                                          gCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcc                                                  >  7:67539/1‑126 (MQ=255)
                                                                                                                          gCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcc                                                  <  8:67539/126‑1 (MQ=255)
                                                                                                                          gCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTAAACAAAAACGTAACCATAGGTGCCAttttt                                                                          >  2:70646/1‑102 (MQ=255)
                                                                                                                            tCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACTATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAccc                         >  5:1238/1‑149 (MQ=255)
                                                                                                                              cTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctc                       >  3:234573/1‑149 (MQ=255)
                                                                                                                                      tGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                 >  6:308006/1‑149 (MQ=255)
                                                                                                                                        gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                 <  8:217513/147‑1 (MQ=255)
                                                                                                                                         ccTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAATGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc              <  5:132332/149‑1 (MQ=255)
                                                                                                                                         ccTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc              >  2:175629/1‑149 (MQ=255)
                                                                                                                                           tctTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc             >  5:189186/1‑148 (MQ=255)
                                                                                                                                           tctTCGCGGGCACGCCTGAGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACAAAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCACAAATTCGTTACAAACTTTAGACACTCTAttttcaac             >  4:370067/1‑144 (MQ=255)
                                                                                                                                           tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCCCCCTCTATTTTGAcc             >  6:224915/1‑148 (MQ=255)
                                                                                                                                            ctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAGTTTTGTGTGCAGTTTCTTCAGTGGCGACCCAATTCGTTACACAGTTTGCACCCTCTATTTTGACCa            >  7:169061/1‑148 (MQ=255)
                                                                                                                                             ttCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttctt                                                            >  6:165039/1‑99 (MQ=255)
                                                                                                                                             ttCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttctt                                                            <  5:165039/99‑1 (MQ=255)
                                                                                                                                                gCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTg                                                       >  2:162790/1‑101 (MQ=255)
                                                                                                                                                gCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTg                                                       <  1:162790/101‑1 (MQ=255)
                                                                                                                                                     cGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAc  <  8:167659/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAC  >  NC_002947/4899921‑4900216

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: