Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,131,310	A→T	9.8%	N98I (AAT→ATT)	flgB →	flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,131,310	0	A	T	9.8%	181.8 / 9.5	61	N98I (AAT→ATT)	flgB	flagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base T (3/3);  total (31/30)
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.86e-01

TTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTT  >  NC_000913/1131163‑1131458                                                                                                                                                    |                                                                                                                                                     ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTaa                                                                                                                                                      >  6:227895/1‑148 (MQ=255) ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGCTaa                                                                                                                                                      >  5:347920/1‑148 (MQ=255)   aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                                                                                                                                   >  8:66498/1‑149 (MQ=255)   aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                                                                                                                                   >  5:70053/1‑149 (MQ=255)     aaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGAATTCCGGACCAGCCTTCGCTTGACGGTAATACCg                                                                                                                                                 >  1:70136/1‑149 (MQ=255)      aaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGt                                                                                                                                                >  7:254105/1‑149 (MQ=255)            aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGatat                                                                                                                                          <  5:227895/149‑1 (MQ=255)            aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGatat                                                                                                                                          <  6:70053/149‑1 (MQ=255)               cAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTCGACGGTATTACCGTcaagcgaa                                                                                                                                        <  2:120305/148‑8 (MQ=255)                   gTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATcgc                                                                                                                                   <  7:23008/149‑1 (MQ=255)                       aCGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                                                                                                                               >  2:55693/1‑149 (MQ=255)                       aCGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                                                                                                                               <  3:335063/149‑1 (MQ=255)                           gATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc                                                                                                                           >  1:222227/1‑149 (MQ=255)                            aTGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc                                                                                                                          <  5:87225/149‑1 (MQ=255)                                aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGt                                                                                                                      >  5:157408/1‑149 (MQ=255)                                          gTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATaa                                                                                                            >  7:355850/1‑149 (MQ=255)                                           ttGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCg                                                                                                                                              <  3:308657/114‑1 (MQ=255)                                           ttGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCg                                                                                                                                              >  4:308657/1‑114 (MQ=255)                                                   aCGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGc                                                                                                    >  6:56470/1‑148 (MQ=255)                                                    cGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                                                                                  <  7:66498/149‑1 (MQ=255)                                                     gATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                                                                                  >  3:203465/1‑148 (MQ=255)                                                     gATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                                                                                  >  4:234993/1‑148 (MQ=255)                                                         aCCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCa                                                                                             <  6:160477/149‑1 (MQ=255)                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATAc                                                                                               >  3:293873/1‑145 (MQ=255)                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           <  4:203465/149‑1 (MQ=255)                                                            tCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGAt                                                                                          >  4:164497/1‑149 (MQ=255)                                                                    acacATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTa                                                                                  <  2:222227/149‑1 (MQ=255)                                                                    acacATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTa                                                                                  <  6:157408/149‑1 (MQ=255)                                                                        aTTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAgcgc                                                                              >  4:183585/1‑149 (MQ=255)                                                                        aTTACGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAgcgc                                                                              >  2:316246/1‑149 (MQ=255)                                                                             ggCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGa                                                                         <  8:269973/149‑1 (MQ=255)                                                                                  aGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCggg                                                                    <  3:183585/149‑1 (MQ=255)                                                                                          aCGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaa                                                                                                                                             >  8:120487/1‑66 (MQ=255)                                                                                          aCGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaa                                                                                                                                             <  7:120487/68‑3 (MQ=255)                                                                                          aCGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCa                                                             >  5:580/1‑148 (MQ=255)                                                                                              ctcctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAgg                                                         >  2:166956/1‑148 (MQ=255)                                                                                                 ctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCa                                                       <  1:55693/147‑1 (MQ=255)                                                                                                 ctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCATCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGCTGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCat                                                      <  5:56470/148‑1 (MQ=255)                                                                                                  tACCGCAGAACTGCAATACCTTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAAAACCAGATGAGCCTTAGCCCGTTGAGCGGGCAAATCAAAGTAatga                                                    >  3:44157/1‑149 (MQ=255)                                                                                                    ccGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgatg                                                  <  1:316246/149‑1 (MQ=255)                                                                                                     cGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACACACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgat                                                   <  8:137278/147‑1 (MQ=255)                                                                                                     cGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACACACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgat                                                   >  7:137278/1‑147 (MQ=255)                                                                                                       cAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAAc                                               <  8:254105/149‑1 (MQ=255)                                                                                                            cTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                          <  3:164497/149‑1 (MQ=255)                                                                                                            cTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                          <  4:293873/149‑1 (MQ=255)                                                                                                             tGCAATACCGTATTCCGGACCAGCCTTCGCTTGTCGGATATACCGTCGATATGGATCGCGAACGCTCCCAGTTTGCCGATAACACCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                         <  6:347920/149‑1 (MQ=255)                                                                                                                   aCCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                          <  8:128531/142‑1 (MQ=255)                                                                                                                   aCCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                          >  7:128531/1‑142 (MQ=255)                                                                                                                   aCCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGc                                   >  4:88693/1‑149 (MQ=255)                                                                                                                             ggACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAAc                         <  8:355850/149‑1 (MQ=255)                                                                                                                                   gCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATgg                   >  2:345046/1‑149 (MQ=255)                                                                                                                                     cTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                          <  7:214847/124‑1 (MQ=255)                                                                                                                                     cTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCa                 <  7:126930/149‑1 (MQ=255)                                                                                                                                     cTTCGCTTGACGCTAATACCGTCGATATGGATCGCGAACGCAACCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGATGAGCGGGCAAATCAAAGGCATCATGAACGattt                                          >  8:214847/1‑124 (MQ=255)                                                                                                                                      ttCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGa                                                                                                               <  5:428/54‑1 (MQ=255)                                                                                                                                      ttCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGa                                                                                                               >  6:428/1‑54 (MQ=255)                                                                                                                                      ttCGCTTGACGGTAATACCGTCGAAATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGCGCAAATCAAAGGCATGATGAACGTTTTACAGAGCTGAAAAAAACGGATGGCa                 >  4:38777/1‑148 (MQ=255)                                                                                                                                         tatcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGtt                                                                           <  5:327256/83‑1 (MQ=255)                                                                                                                                         tatcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGtt                                                                           >  6:327256/5‑87 (MQ=255)                                                                                                                                         tatcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCActgc             >  1:240953/5‑149 (MQ=255)                                                                                                                                         gCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCActgc             >  7:260497/1‑149 (MQ=255)                                                                                                                                              aCGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAat        <  3:234993/149‑1 (MQ=255)                                                                                                                                                ggTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAatat      >  2:200897/1‑149 (MQ=255)                                                                                                                                                   aaTACCGTCGATATGGATCGCGAACGCCCCCAGTTTGCCGATAACACCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCAGTGCTGAATAtttt   >  8:231839/1‑149 (MQ=255)                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  1:166956/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTT  >  NC_000913/1131163‑1131458

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