Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,064,160	T→C	14.5%	L42P (CTG→CCG) ‡	yccE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,160	0	T	C	14.5%	170.7 / 21.3	69	L42P (CTG→CCG) ‡	yccE	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (35/24);  new base C (5/5);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-02

TTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATATTGCGTA  >  NC_000913/1064021‑1064306                                                                                                                                            |                                                                                                                                                   ttATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttt                                                                                                                                           <  3:23221/149‑1 (MQ=255) ttATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttt                                                                                                                                           <  3:462873/149‑1 (MQ=255)  tataGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttt                                                                                                                                           <  1:410690/148‑1 (MQ=255)         gCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                                                                   >  3:188186/1‑149 (MQ=255)             tatatGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTtata                                                                                                                               <  3:728123/149‑1 (MQ=255)              atatGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATAc                                                                                                                              <  3:220517/149‑1 (MQ=255)                 tGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAAt                                                                                                                           <  3:840437/149‑1 (MQ=255)                     tagtaATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                                                                        >  1:622734/1‑148 (MQ=255)                                aTGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTAAGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaa                                                                                                            >  3:92524/1‑149 (MQ=255)                                 tGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                                                                           >  4:380264/1‑149 (MQ=255)                                    aaTTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagta                                                                                                         >  1:577697/1‑148 (MQ=255)                                      ttAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAAt                                                                                                      >  4:721779/1‑149 (MQ=255)                                       tAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACaa                                                                                                                            >  1:305569/1‑126 (MQ=255)                                       tAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACaa                                                                                                                            <  2:305569/126‑1 (MQ=255)                                                tGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGAt                                                                                             >  2:197829/1‑148 (MQ=255)                                                            tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTgg                                                                                >  4:870761/1‑149 (MQ=255)                                                             ttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTgg                                                                                >  1:772733/1‑148 (MQ=255)                                                             ttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGtttt                                                                                                                                          <  2:178411/90‑1 (MQ=255)                                                             ttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGtttt                                                                                                                                          >  1:178411/1‑90 (MQ=255)                                                                 aaCAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                                                                           >  3:82622/1‑149 (MQ=255)                                                                 aaCAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtg                                                                            <  2:622734/148‑1 (MQ=255)                                                                    aGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGa                                                                         >  2:780227/1‑148 (MQ=255)                                                                       ccTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATc                                                                      >  1:322489/1‑148 (MQ=255)                                                                           gAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                                                                 <  3:870761/149‑1 (MQ=255)                                                                            aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                                                                 >  2:865768/1‑148 (MQ=255)                                                                             aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                                                                        >  4:138345/1‑92 (MQ=255)                                                                             aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                                                                        <  3:138345/92‑1 (MQ=255)                                                                              cGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa                                                               >  2:148506/1‑148 (MQ=255)                                                                                atatAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa                                                             <  1:389793/148‑1 (MQ=255)                                                                                atatAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa                                                             <  2:322489/148‑1 (MQ=255)                                                                                   taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGaa                                                                        <  3:94137/134‑1 (MQ=255)                                                                                   taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGaa                                                                        >  4:94137/1‑134 (MQ=255)                                                                                    aaaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaa                                                                                                            >  4:237161/1‑97 (MQ=255)                                                                                    aaaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaa                                                                                                            <  3:237161/97‑1 (MQ=255)                                                                                    aaaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGaa                                                         <  1:795524/148‑1 (MQ=255)                                                                                        aaTGAGTAAGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGa                                                    <  4:92524/149‑1 (MQ=255)                                                                                         aTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg                                                   >  3:149682/1‑149 (MQ=255)                                                                                         aTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg                                                   >  3:340337/1‑149 (MQ=255)                                                                                             gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGaaa                                               >  4:437131/1‑149 (MQ=255)                                                                                              tACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAAt                                              >  4:294517/1‑149 (MQ=255)                                                                                                 gAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCa                                           >  4:258942/1‑149 (MQ=255)                                                                                                    aaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaa                                        <  4:188186/149‑1 (MQ=255)                                                                                                    aaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGa                                         >  2:314296/1‑148 (MQ=255)                                                                                                       aaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatg                                      <  1:125409/148‑1 (MQ=255)                                                                                                        aTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgat                                    >  4:128049/1‑149 (MQ=255)                                                                                                            aaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatg                                >  3:638649/1‑149 (MQ=255)                                                                                                            aaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataat                                 >  2:584543/1‑148 (MQ=255)                                                                                                             aCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTAAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatg                                >  1:159138/1‑148 (MQ=255)                                                                                                              cATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatgat                              <  3:721779/149‑1 (MQ=255)                                                                                                                 atTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTa                           >  3:678528/1‑149 (MQ=255)                                                                                                                  ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatgat                              <  1:597899/145‑1 (MQ=255)                                                                                                                  ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatgat                              >  2:597899/1‑145 (MQ=255)                                                                                                                  ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTa                           <  1:627072/148‑1 (MQ=255)                                                                                                                      aTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCatt                      >  3:175402/1‑149 (MQ=255)                                                                                                                       tCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCatt                      >  2:628317/1‑148 (MQ=255)                                                                                                                         aCGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCattat                    <  1:197829/148‑1 (MQ=255)                                                                                                                            cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCAttatttat                <  3:294517/149‑1 (MQ=255)                                                                                                                            cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCAttatttat                <  4:340337/149‑1 (MQ=255)                                                                                                                            cTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCAttatttat                <  4:149682/149‑1 (MQ=255)                                                                                                                              ttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGc                                            >  1:837190/1‑119 (MQ=255)                                                                                                                              ttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGc                                            <  2:837190/119‑1 (MQ=255)                                                                                                                                  taaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa                                                             <  1:424507/97‑1 (MQ=255)                                                                                                                                  taaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa                                                             >  2:424507/2‑98 (MQ=255)                                                                                                                                     aaCACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAatat        <  1:628317/148‑1 (MQ=255)                                                                                                                                     aaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAatat        >  1:513651/1‑148 (MQ=255)                                                                                                                                     aaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAatat        >  1:901934/1‑148 (MQ=255)                                                                                                                                     aaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAatat        >  2:482096/1‑148 (MQ=255)                                                                                                                                         cacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATATTGCg    >  2:356998/1‑148 (MQ=255)                                                                                                                                          acTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATATTGCGTa  >  3:387118/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                   TTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATATTGCGTA  >  NC_000913/1064021‑1064306

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