Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,491 C→T 5.9% intergenic (+61/‑46) ycaC‑I → / → ppkB putative hydrolase/polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,4910CT 5.9% 200.0 / 2.0 68intergenic (+61/‑46)ycaC‑I/ppkBputative hydrolase/polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/2);  total (34/34)
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.84e-01

TCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAACCTGCATGTAGAGTTGG  >  NC_002947/827348‑827633
                                                                                                                                                |                                                                                                                                              
tCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTaa                                                                                                                                           <  7:174659/149‑1 (MQ=255)
tCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTaa                                                                                                                                           >  2:241313/1‑149 (MQ=255)
         tGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGTAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCAt                                                                                                                                   >  4:36491/1‑148 (MQ=255)
          gggTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATac                                                                                                                                 >  2:123433/1‑149 (MQ=255)
                ccTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATactcactc                                                                                                                           >  4:188316/1‑149 (MQ=255)
                ccTGTTCTCCAACCATATCCCGGACTAACCCAACCTGATGACCAGCTACAAAAAGCTCGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGCCATAAATCGACCATcctaaatc                                                                                                                           >  8:269442/1‑144 (MQ=38)
                  tGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCActct                                                                                                                          >  6:4851/1‑148 (MQ=255)
                         cAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGAcc                                                                                                                                     >  1:57781/1‑130 (MQ=255)
                         cAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGAcc                                                                                                                                     <  2:57781/130‑1 (MQ=255)
                         cAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGc                                                                                                                  >  6:211225/1‑149 (MQ=255)
                               tatCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAg                                                                                                            <  8:314114/149‑1 (MQ=255)
                               tatCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAg                                                                                                            >  1:316086/1‑149 (MQ=255)
                                 tCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAg                                                                                                            >  4:226672/1‑147 (MQ=255)
                                   ccGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGc                                                                                                        >  3:252551/1‑149 (MQ=255)
                                     ggACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCgt                                                                                                      >  8:75668/1‑149 (MQ=255)
                                       aCTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCgtgt                                                                                                    <  5:211225/149‑1 (MQ=255)
                                         tACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGCGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATact                                                                                                                                <  7:137212/119‑1 (MQ=255)
                                         tACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGACTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTCCAAGCGAAGCGGCAAACTGGCATTAATCGACCATact                                                                                                                                >  8:137212/1‑119 (MQ=255)
                                           ccGGAACCTGATGACCAGCTACAACAAGCTTGCCAAATAGTGCCTTCTGGTGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATg                                                                                                <  7:326130/149‑1 (MQ=255)
                                           ccGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTACTGGGGGCTGCACGACCGCCCCCACATCAACCTGCAAGCGAAGCGGCAAACTGGCATTAATCGCCCATACTCACTCTACATTAGCCTTCAGGAGCGTGTCAt                                                                                                 >  1:250156/1‑149 (MQ=255)
                                           ccGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCtt                                                                                                               >  2:36206/1‑134 (MQ=255)
                                           ccGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCtt                                                                                                               <  1:36206/134‑1 (MQ=255)
                                           ccGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTc                                                                                                              <  6:337222/135‑1 (MQ=255)
                                           ccGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTc                                                                                                              >  5:337222/1‑135 (MQ=255)
                                           acGCAAACTGATGACCAGCTAGAACAAGCTGGCCAAATAGTGCCTTCTGGGGGGTCCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCATCATTAGCCTTCAGGAGCGTGTCa                                                                                                  <  8:127636/146‑1 (MQ=255)
                                             gtaaCCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTc                                                                                              <  3:36491/147‑1 (MQ=255)
                                                aCCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGaa                                                                                           <  7:75668/149‑1 (MQ=255)
                                                aCCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGaa                                                                                           <  4:252551/149‑1 (MQ=255)
                                                   tgatgaCCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTc                                                                                                              >  4:207030/1‑127 (MQ=255)
                                                   tgatgaCCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTc                                                                                                              <  3:207030/127‑1 (MQ=255)
                                                        aCCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCa                                                                                   <  3:188316/149‑1 (MQ=255)
                                                         ccAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAg                                                                                  <  7:169782/149‑1 (MQ=255)
                                                          cAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAg                                                                                  >  2:114769/1‑148 (MQ=255)
                                                            gCTACAACCAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGaa                                                                               <  4:271694/149‑1 (MQ=255)
                                                            gCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGaa                                                                               <  5:4851/149‑1 (MQ=255)
                                                                caacaaGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGccc                                                                           <  7:269442/149‑1 (MQ=255)
                                                                caacaaGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGcc                                                                            <  2:316086/148‑1 (MQ=255)
                                                                           ccAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGTATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCaaaaaa                                                                  >  5:321349/1‑147 (MQ=255)
                                                                           ccAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGTATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACa                                                                >  5:41742/1‑149 (MQ=255)
                                                                           ccAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGGGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACa                                                                >  5:87594/1‑149 (MQ=255)
                                                                           ccAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGCAGCCCCCCAAAAAACa                                                                >  5:53110/1‑149 (MQ=255)
                                                                            cAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAg                                                               >  3:317823/1‑149 (MQ=255)
                                                                            cAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAg                                                               >  1:185148/1‑149 (MQ=255)
                                                                             aaaTAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGc                                                              <  5:171827/149‑1 (MQ=255)
                                                                               aTAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCa                                                            <  1:123433/149‑1 (MQ=255)
                                                                                   tGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGaaa                                                        <  2:185148/149‑1 (MQ=255)
                                                                                           tGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGTATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTggg                                                 <  6:41742/149‑1 (MQ=255)
                                                                                             ggggCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAATTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTgg                                                  <  1:114769/146‑1 (MQ=255)
                                                                                             ggggCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGTATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTggg                                                 <  6:321349/147‑1 (MQ=255)
                                                                                             ggggCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTgg                                                  <  8:254066/146‑1 (MQ=255)
                                                                                             ggggCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCg                                               >  2:89927/1‑149 (MQ=255)
                                                                                             ggggCTGCACGACAGCCCCCACATCAAGCTGCAAGCGAAGCGGAAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTggg                                                 <  2:121743/147‑1 (MQ=255)
                                                                                                  tGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAg                                          >  3:135572/1‑149 (MQ=255)
                                                                                                    cACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAg                                          <  3:226672/147‑1 (MQ=255)
                                                                                                    cACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGc                                        >  4:131469/1‑149 (MQ=255)
                                                                                                                      aaccTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCaa                      <  2:250156/146‑1 (MQ=255)
                                                                                                                      aaTCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGt                                             >  7:140969/1‑126 (MQ=255)
                                                                                                                      aaTCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGt                                             >  1:101886/1‑126 (MQ=255)
                                                                                                                      aaTCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGt                                             <  2:101886/126‑1 (MQ=255)
                                                                                                                      aaTCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGt                                             <  8:140969/126‑1 (MQ=255)
                                                                                                                         cTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGAGTCATGACGACGCCGTCCAGGAAGCCGCCCAAAAACCACCCAGAAAACCTGGGAGGAAAGGCGTACAAAAACCCGATCAAGaa                   >  6:329517/1‑149 (MQ=255)
                                                                                                                          tGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAAc                  <  3:131469/149‑1 (MQ=255)
                                                                                                                          tGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAAc                  >  8:165192/1‑149 (MQ=255)
                                                                                                                           gCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAAcc                 >  1:251993/1‑149 (MQ=255)
                                                                                                                                 gaagcgGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGaa                                                                                           <  3:157234/69‑1 (MQ=255)
                                                                                                                                 gaagcgGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGaa                                                                                           >  4:157234/1‑69 (MQ=255)
                                                                                                                                        cAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAACCTGCATGTAGAGtt    >  4:137487/1‑149 (MQ=255)
                                                                                                                                          aaCTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAACCTGCATGTAGAGTTgg  <  2:338451/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCAATCTGCAAGCGAAGCGGCAAACTGGCATTAATCGACCATACTCACTCTTCATTAGCCTTCAGGAGCGTGTCATGTCGAAGCCGTCCAGGAAGCCCCCCAAAAAACAGCCAGAAAAACTGGGCGGTAAGGCGTACGAAAAAGCGCTCAAGAACCTGCATGTAGAGTTGG  >  NC_002947/827348‑827633

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: