Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,930,409	T→A	9.0%	T188T (ACA→ACT) ‡	cheZ ←	protein phosphatase CheZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,930,409	0	T	A	9.0%	168.3 / 4.2	61	T188T (ACA→ACT) ‡	cheZ	protein phosphatase CheZ
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (4/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCG  >  NC_002947/4930265‑4930557                                                                                                                                                 |                                                                                                                                                     cGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                                      <  4:161698/145‑1 (MQ=255) cGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTgtgg                                                                                                                                                  <  5:194538/149‑1 (MQ=255)     ccTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGAc                                                                                                                                              <  2:41199/149‑1 (MQ=255)       atccccccGAGTCGTCTGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGAcgc                                                                                                                                            <  8:139322/145‑1 (MQ=255)          aCCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCACACTGGTGAccgagg                                                                                                                                         <  3:771/149‑6 (MQ=255)             ccGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGa                                                                                                                                      >  8:28232/1‑149 (MQ=255)              cGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGa                                                                                                                                      >  5:52181/1‑148 (MQ=255)               gAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGAt                                                                                                                                     >  1:114790/1‑148 (MQ=255)                    ggATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTg                                                                                                                               >  6:247949/1‑149 (MQ=255)                    ggATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTg                                                                                                                               >  3:184367/1‑149 (MQ=255)                      aTGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGc                                                                                                                             <  5:148356/149‑1 (MQ=255)                         tttttCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCg                                                                                                                           <  7:28232/148‑1 (MQ=255)                         tttttCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCg                                                                                                                           <  4:263716/148‑1 (MQ=255)                           tttCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGt                                                                                                                         <  4:301763/148‑1 (MQ=255)                                cGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAgg                                                                                                                     <  7:58040/147‑1 (MQ=255)                                cGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAgg                                                                                                                     >  8:58040/1‑147 (MQ=255)                                  aTCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCt                                                                                                                 >  7:253000/1‑149 (MQ=255)                                     attttCAGCACGCAAATGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGa                                                                                                               <  8:71557/147‑1 (MQ=255)                                      ttttCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTCGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGa                                                                                                               >  7:224351/1‑147 (MQ=255)                                       tttCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTcc                                                                                                                  <  4:184367/143‑1 (MQ=255)                                       tttCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGa                                                                                                               <  7:218941/146‑1 (MQ=255)                                       tttCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATa                                                                                                             <  6:52181/148‑1 (MQ=255)                                             cACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCt                                                                                                        <  1:73378/147‑1 (MQ=255)                                             cACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTgg                                                                                                      <  1:301104/149‑1 (MQ=255)                                                gCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAg                                                                                                                      >  7:263930/1‑130 (MQ=255)                                                gCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAg                                                                                                                      <  8:263930/130‑1 (MQ=255)                                                    ttGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGa                                                                                                               <  1:288294/133‑1 (MQ=255)                                                    ttGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGa                                                                                                               >  2:288294/1‑133 (MQ=255)                                                           tGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTCAGCAGGTTGCTTTACACCTCGGTCACCAGAGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCATGGGCCAGCAAAATGTc                                                                                        >  8:175449/1‑149 (MQ=255)                                                             gTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                                      >  5:241395/1‑85 (MQ=255)                                                             gTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAgt                                                                                                                                                      <  6:241395/85‑1 (MQ=255)                                                                aTGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTCGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTg                                                                                    <  8:224351/148‑1 (MQ=255)                                                                 tGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAg                                                                                  >  3:226794/1‑149 (MQ=255)                                                                      tataCCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGc                                                                              <  5:247949/148‑1 (MQ=255)                                                                          ccGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTcc                                                                                                                  >  5:240611/1‑108 (MQ=255)                                                                          ccGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTcc                                                                                                                  <  6:240611/108‑1 (MQ=255)                                                                          ccGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCgg                                                                          >  7:234106/1‑148 (MQ=255)                                                                           cGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGa                                                                                   >  6:285848/1‑138 (MQ=255)                                                                           cGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGa                                                                                   <  5:285848/138‑1 (MQ=255)                                                                           cGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGa                                                                         >  5:26515/1‑148 (MQ=255)                                                                            ggCAAAGCGGTCGACCTGGCTGGCGATCAGCACCAGCTTGAGCAGGTTGCTTTCCTCCTCGGTCACCAGTGTGGTGACGCGGTTGATCACGTGGCGGGTCAGGTGCTGATAGTCCTGGGGCATCAAATTGTCGTTGAGGTGCCCGGaaa                                                                       <  7:175449/149‑1 (MQ=255)                                                                                aaGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCtt                                                                   <  3:215752/149‑1 (MQ=255)                                                                                         aCCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTAcc                                                          >  8:36601/1‑149 (MQ=255)                                                                                                tggcCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGgcg                                                   >  3:149536/1‑149 (MQ=255)                                                                                                   ccATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCtg                                                <  8:210143/149‑1 (MQ=255)                                                                                                    cATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCtgt                                               >  6:36912/1‑149 (MQ=255)                                                                                                      tCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGc                                             <  5:36912/149‑1 (MQ=255)                                                                                                      tCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGc                                             <  8:253000/149‑1 (MQ=255)                                                                                                      tCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGc                                             <  3:234482/149‑1 (MQ=255)                                                                                                         gCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGt                                           >  1:136096/1‑148 (MQ=255)                                                                                                         gCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGt                                           >  1:140766/1‑148 (MQ=255)                                                                                                          cACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCa                                         <  6:26515/149‑1 (MQ=255)                                                                                                           aCCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAg                                        >  2:284988/1‑149 (MQ=255)                                                                                                                     gCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGc                                                                 >  5:25153/1‑114 (MQ=255)                                                                                                                     gCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGc                                                                 <  6:25153/114‑1 (MQ=255)                                                                                                                        ggTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACt                           <  2:114790/149‑1 (MQ=255)                                                                                                                            gCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCt                       >  1:78102/1‑149 (MQ=255)                                                                                                                                  caagcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGAc                   >  3:3231/8‑147 (MQ=255)                                                                                                                                   aCCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGAc                   >  6:5983/1‑146 (MQ=255)                                                                                                                                   aCCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAg                <  2:78102/149‑1 (MQ=255)                                                                                                                                    ccTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAg                >  4:263955/1‑148 (MQ=255)                                                                                                                                     gcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCa                                         <  7:25391/118‑1 (MQ=255)                                                                                                                                     gcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCa                                         >  8:25391/5‑122 (MQ=255)                                                                                                                                     gcgcgTCACCACACTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGAc                   >  6:173044/5‑144 (MQ=255)                                                                                                                                                 tgtgGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCg  <  4:226794/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGGACCTTCACCCCGAGTCGGATGTTTTTCTCGATCTTTTTCTGCACGCAATTGATCGTGGTCATGTTTTATACCGGCAAAGCGGTCGACCTGGCTGGCCATCAGCACCAGCTTGAGCAGGTTGCTTTCCACCTCGGTCACCAGTGTGGTGACGCGCTTGATCACCTGGCCGGTCAGGTCCTGATAGTCCTGGGCCAGCAAAATGTCGTTGAGGTGCCCGGAAACCTTGCGGTTACCCTCGGCGCTGTGCGTCAGGAAACTGTCGACTCGCTTGACCAGATCACGAAATTCCG  >  NC_002947/4930265‑4930557

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