Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,027,735	A→G	C13C (TGT→TGC)	yccT ←	UPF0319 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,027,735	0	A	G	95.2%	219.7 / ‑1.3	62	C13C (TGT→TGC)	yccT	UPF0319 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (25/34);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

GCGAAACACTAACTGATGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAACATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGCTCTCATTCATTCGT  >  NC_000913/1027588‑1027882                                                                                                                                                    |                                                                                                                                                    gcgAAACACTAACTGATGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCAt                                                                                                                                                  >  2:300822/1‑151 (MQ=255)      acacTAACTGATGTGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGc                                                                                                                                             <  1:406372/151‑1 (MQ=255)       cacTAACTGATGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGc                                                                                                                                             >  2:322416/1‑150 (MQ=255)              tGATGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAgg                                                                                                                                     >  1:403787/1‑151 (MQ=255)                aTGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAggtg                                                                                                                                   <  2:1109873/151‑1 (MQ=255)                aTGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAggtg                                                                                                                                   <  1:198056/151‑1 (MQ=255)                aTGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAggtg                                                                                                                                   >  2:774881/1‑151 (MQ=255)                 tGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAggtgg                                                                                                                                  >  2:69581/1‑151 (MQ=255)                       ccATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACg                                                                                                                             <  2:1175958/150‑1 (MQ=255)                        cATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGAt                                                                                                                           <  2:186452/151‑1 (MQ=255)                            aTCCGGTTCAATGCGGTCCACACCGCGCAACAGAGGGCTGGAAGCCGGCTTGCCGTCGGGTACGAGAAGATCGCGATCGGTTGAAAGCCGCAGCGTGGTTGCAAAAACGGCTACCGTCAAGCGAAATGCAATCAAGGTCGTCGCGATgg                                                                                                                         <  1:1132427/149‑2 (MQ=255)                              ccAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCgg                                                                                                                      <  2:264774/150‑1 (MQ=255)                                 gTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGtttt                                                                                                                  >  1:899764/1‑151 (MQ=255)                                    cAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCa                                                                                                                >  1:963114/1‑150 (MQ=255)                                    cAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAACATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCAt                                                                                                               >  1:985664/1‑151 (MQ=255)                                     aaTGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATa                                                                                                              <  2:1070860/151‑1 (MQ=255)                                        gCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGt                                                                                                           <  1:594181/151‑1 (MQ=255)                                               gCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCt                                                                                                    <  1:1178807/151‑1 (MQ=255)                                                cACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTg                                                                                                   >  2:1125940/1‑151 (MQ=255)                                                      gcagcaGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCaaaa                                                                                             <  2:820731/151‑1 (MQ=255)                                                      gcagcaGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCaaa                                                                                              >  2:1140346/1‑150 (MQ=255)                                                       cagcagAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCaaa                                                                                              >  2:682274/1‑149 (MQ=255)                                                        agcagAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCaaa                                                                                              >  1:364243/1‑148 (MQ=255)                                                          cagAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGTcc                                                                                         <  1:355404/151‑1 (MQ=255)                                                          cagAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGAcc                                                                                         >  1:983538/1‑151 (MQ=255)                                                             agagCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGcc                                                                                      >  2:761931/1‑151 (MQ=255)                                                                      agccTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAACATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGaataa                                                                               <  1:121815/147‑1 (MQ=255)                                                                        ccTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGaataataat                                                                           <  2:3260/151‑1 (MQ=255)                                                                         cttcttGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAAtg                                                                          >  2:1085988/1‑151 (MQ=255)                                                                           tcttGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAAtgt                                                                         <  2:774940/150‑1 (MQ=255)                                                                           tcttGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAAtg                                                                          <  2:709592/149‑1 (MQ=255)                                                                            cttGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATgtgt                                                                       >  2:794122/1‑151 (MQ=255)                                                                               gCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCa                                                                     <  2:316869/150‑1 (MQ=255)                                                                                 cGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCaa                                                                    <  2:612817/149‑1 (MQ=255)                                                                                    cGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATTCCGGTTTTCATAGGTAGCTCCTGTAAAAAGACCGCCCGAATAATAATGTGTCAAAAttt                                                               >  1:183444/1‑151 (MQ=255)                                                                                     gAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAAtttt                                                              >  1:411996/1‑151 (MQ=255)                                                                                      aGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAAttttt                                                             <  2:136941/151‑1 (MQ=255)                                                                                         aCCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCtt                                                          <  1:1215519/151‑1 (MQ=255)                                                                                          ccAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCttt                                                         <  1:247805/151‑1 (MQ=255)                                                                                            agaagaTCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTa                                                        >  2:868157/1‑150 (MQ=255)                                                                                             gaagaTCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAAt                                                      <  2:156317/151‑1 (MQ=255)                                                                                               agaTCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCa                                                    <  2:433835/151‑1 (MQ=255)                                                                                                gaTCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCAt                                                   >  1:78926/1‑151 (MQ=255)                                                                                                   cGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCAtact                                                <  1:1212964/151‑1 (MQ=255)                                                                                                      cATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCAtacttac                                             <  2:316896/151‑1 (MQ=255)                                                                                                           gTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACa                                        <  1:218336/151‑1 (MQ=255)                                                                                                             tGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACat                                       <  1:649270/150‑1 (MQ=255)                                                                                                              gAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACatat                                     <  1:270334/151‑1 (MQ=255)                                                                                                                aaGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACatat                                     >  1:344215/1‑149 (MQ=255)                                                                                                                    cGCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGAc                                 <  1:435420/149‑1 (MQ=255)                                                                                                                     gCAGCGTGGTGGCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACat                               >  1:1251767/1‑150 (MQ=255)                                                                                                                               ggCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCt                     <  1:636779/150‑1 (MQ=255)                                                                                                                                gCAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCttt                   <  1:539327/151‑1 (MQ=255)                                                                                                                                gCAAAAACGGATACCGGCAAACATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAAT‑ATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCtttt                  <  2:42385/151‑1 (MQ=255)                                                                                                                                 cAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCtttt                  <  2:618063/151‑1 (MQ=255)                                                                                                                                 cAAAAACGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCtttt                  <  2:388664/151‑1 (MQ=255)                                                                                                                                  aaaaaCGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTg                 <  1:314975/151‑1 (MQ=255)                                                                                                                                   aaaaCGGATACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGc                >  2:1037892/1‑151 (MQ=255)                                                                                                                                           tACCGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATCTAGACATATTTCCCGCTTTTGCTCtcattc        >  2:687411/1‑151 (MQ=255)                                                                                                                                             ccGGCAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGCTCtcattcat      <  1:1140282/151‑1 (MQ=255)                                                                                                                                                 cAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGCTCTCATTCATTCGt  >  2:255791/1‑151 (MQ=255)                                                                                                                                                 cAAGCATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGCTCTCATTCATTCGt  >  2:463355/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGAAACACTAACTGATGCGGTCCATTATCCAGTTCAATGCTGTCCGCACCGCGCAGCAGAGAGCTGGAAACCTTCTTGCCGTCGAGTACCAGAAGATCGACATCGGTTGAAAGCCGCAGCGTGGTGGCAAAAACGGATACCGGCAAACATAATGCAATCAAGGTGGTCACGATGCCGGTTTTCATAGGTAGCTCCTGTCAAAAGACCGCCCGAATAATAATGTGTCAAAATTTTTCTTTAATCATACTTACTAACATATAGACATATTTCCCGCTTTTGCTCTCATTCATTCGT  >  NC_000913/1027588‑1027882

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