Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,401,459	0	T	C	26.3%	114.5 / 15.2	57	D625G (GAT→GGT)	csrD	regulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base T (37/5);  new base C (0/15);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 7.03e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCG  >  NC_000913/3401322‑3401598                                                                                                                                          |                                                                                                                                             cGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcc                                                                                                                                            <  2:510450/140‑1 (MQ=255)        aCAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGcc                                                                                                                                    <  2:28859/140‑1 (MQ=255)                  tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAac                                                                                                                          >  2:597869/1‑140 (MQ=255)                       tGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACcgc                                                                                                                     >  1:772477/1‑140 (MQ=255)                          gaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCt                                                                                                                  >  1:971098/1‑140 (MQ=255)                          gaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCt                                                                                                                  >  2:762946/1‑140 (MQ=255)                             gCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGaa                                                                                                               >  1:322577/1‑140 (MQ=255)                              caccccgcATTTTTCTAGGTGAACACAGTTAAAGCGCGTGGTTAAGCGGCGTACCTTTAAAAATATTTTTGCGCGTTAGTATCAAGGGGCGGTGAGGGCGCAAAAAACCCCCC‑TTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                              <  2:131907/135‑1 (MQ=255)                              cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                              >  1:921261/1‑140 (MQ=255)                              cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                              >  1:535525/1‑140 (MQ=255)                                  gggttttattCTAGGAAAAAACGGATTAATCGGCAGTTTAACCGGATTATCTTTGTAATTTTTTTTTCTGCTTAGTATCACTTGGTTGTAGGGCCGAAAAAAACCCCCC‑TTCCCCGCCTGTAACACCGCGCTGAATCAgg                                                                                                          <  2:971098/135‑1 (MQ=255)                                   cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                                         >  1:217007/1‑140 (MQ=255)                                   cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                                         >  1:871984/1‑140 (MQ=255)                                      attattCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTTCCCGCCTGTAACACCGCGCTGAATCAGGGCCt                                                                                                      >  1:487731/1‑140 (MQ=255)                                                  gAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                                          >  1:806533/1‑140 (MQ=255)                                                    aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcgc                                                                                        >  2:45537/1‑140 (MQ=255)                                                                gacagggTAAACCCGAGTTTCTTTGTAAATATTTTTCCACGTTAGTATAAGGGGCCTGTGAGGACGCAAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCg                                                                            <  2:761213/135‑1 (MQ=255)                                                                 ggCAGGTAAAACCGATTATCTTTGTGAATATTTTTTCACGTTATTATCAAGTGGCTGTGAGGCCGCGAAAAAACCCCC‑TGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCgg                                                                           <  1:597869/140‑1 (MQ=255)                                                                      gtaaaaCCGATTATCTTTGTAAATTTTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCg                                                                      <  2:384862/137‑1 (MQ=255)                                                                          acccggatATTTTTGAGAATTTTTTTTCATTTTAGTGTAAAGTGGGTGTG‑GGCGCCGAAAAAAACCCCGTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTaaa                                                                  <  2:806533/133‑1 (MQ=255)                                                                             cGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCt                                                               <  1:762946/140‑1 (MQ=255)                                                                              gAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTg                                                              >  2:764215/1‑140 (MQ=255)                                                                               aGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTgg                                                             >  2:324922/1‑140 (MQ=255)                                                                                gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTgtc                                                            >  1:761213/1‑138 (MQ=255)                                                                                gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                            >  2:678142/1‑140 (MQ=255)                                                                                gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                            >  2:430018/1‑140 (MQ=255)                                                                                gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                            >  2:324387/1‑140 (MQ=255)                                                                                   tttttGTGAATTTTTTTTCCCGTTAGTATCAAGTGGTTGTGAGGACGCAAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                                         <  2:239258/138‑1 (MQ=255)                                                                                    ttttGTAAATATTTTTTCACGTTATTATCAAGGGGCTGTGAGGACGCAAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCt                                                        <  2:884084/139‑1 (MQ=255)                                                                                       tgtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGt                                                     >  2:172811/1‑140 (MQ=255)                                                                                         tgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                                   >  1:966561/1‑140 (MQ=255)                                                                                         taaaTTTTTTTCCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                                   <  2:43675/138‑1 (MQ=255)                                                                                           aaTATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCg                                                 >  2:721752/1‑140 (MQ=255)                                                                                           aaTATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCg                                                 >  1:397932/1‑140 (MQ=255)                                                                                               ttttttCACGTTAGTATAAAGGGGCTGGGGGGACGCAAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGc                                             <  2:139522/140‑1 (MQ=255)                                                                                                   ttCCCGTTAGTATAAGGTGGCTGGGAGGACGCAAAAAACCCCCT‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                                         <  2:217007/140‑1 (MQ=255)                                                                                                    tCACTTTATTATCAGGGGCTTGTGCGGACGCAAAAAAACCCCC‑TTCCCCGCCTGTAACCCCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCt                                        <  2:772477/140‑1 (MQ=255)                                                                                                       cGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTcc                                     >  1:16222/1‑140 (MQ=255)                                                                                                       cGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTcc                                     >  1:772891/1‑140 (MQ=255)                                                                                                          tAGTATCAGGTGGCTGTGAGGACGCAAAAAACCCCCC‑TTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacc                                  <  2:772891/140‑1 (MQ=255)                                                                                                              gacaggTGGCTGTGAGGACGCAAAAAACCCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGc                              <  2:748776/135‑1 (MQ=255)                                                                                                                  aGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                          >  1:765212/1‑140 (MQ=255)                                                                                                                  aGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                          >  2:527973/1‑140 (MQ=255)                                                                                                                  aGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                          >  1:128890/1‑140 (MQ=255)                                                                                                                  aGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                          >  1:760311/1‑140 (MQ=255)                                                                                                                    tGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa                        >  2:685756/1‑140 (MQ=255)                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                      >  2:636926/1‑140 (MQ=255)                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                      >  2:652586/1‑140 (MQ=255)                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                      >  2:903561/1‑140 (MQ=255)                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                      >  2:520183/1‑140 (MQ=255)                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                      >  1:912880/1‑140 (MQ=255)                                                                                                                       cTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                     >  2:660460/1‑140 (MQ=255)                                                                                                                       cTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                     <  1:527973/140‑1 (MQ=255)                                                                                                                       cTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                     <  1:324387/140‑1 (MQ=255)                                                                                                                       cTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                     >  1:814695/1‑140 (MQ=255)                                                                                                                        tgtgAGGACGCGAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcag                    <  2:814695/140‑1 (MQ=255)                                                                                                                         gtgAGGACGCGAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                   <  1:430018/140‑1 (MQ=255)                                                                                                                            gggACGCAAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCagcagc                <  2:397932/139‑1 (MQ=255)                                                                                                                            aGGACGCGAAAAAACCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCagcagc                <  1:324922/140‑1 (MQ=255)                                                                                                                              gACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagcctgtc              >  1:748776/1‑135 (MQ=255)                                                                                                                                 gcaaaaaaaCCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGtt           <  1:394391/137‑1 (MQ=255)                                                                                                                                         aTCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTc   <  1:764215/140‑1 (MQ=255)                                                                                                                                          tCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCg  <  2:921261/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             CGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCC‑TTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCG  >  NC_000913/3401322‑3401598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑