Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,065 A→G 9.5% intergenic (+26/‑472) 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,0650AG 9.5% 199.6 / 8.9 73intergenic (+26/‑472)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (3/4);  total (35/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 = 3.67e-01

CTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCA  >  NC_002947/4899917‑4900206
                                                                                                                                                      |                                                                                                                                              
cTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCa                                                                                                                                                <  6:358668/149‑1 (MQ=255)
   cATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGGGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCa                                                                                                                                                <  4:176522/146‑1 (MQ=255)
   cATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGc                                                                                                                                             <  5:117952/149‑1 (MQ=255)
    aTGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGc                                                                                                                                             >  8:286460/1‑148 (MQ=255)
    aTGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGc                                                                                                                                             <  7:286460/148‑1 (MQ=255)
    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                                                                                                                             <  3:360826/136‑1 (MQ=255)
                                ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                             >  4:360826/1‑136 (MQ=255)
                                 cTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAg                                                                                                                                               >  3:317271/1‑117 (MQ=255)
                                 cTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAg                                                                                                                                               <  4:317271/117‑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:262680/149‑1 (MQ=255)
                                              tGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                 <  3:216826/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)
                                                                 cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCGGGcac                                                                                                                                           <  8:136698/89‑3 (MQ=255)
                                                                 cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCGGGcac                                                                                                                                           >  7:136698/1‑87 (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                                        <  8:67539/126‑1 (MQ=255)
                                                                                                                              gCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcc                                        >  7:67539/1‑126 (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)
                                                                                                                                                      |                                                                                                                                              
CTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGT‑‑GCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCA  >  NC_002947/4899917‑4900206

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