Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,781,767	T→C	F125L (TTT→CTT)	ydiS →	putative oxidoreductase with FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,781,767	0	T	C	97.0%	217.0 / ‑3.9	66	F125L (TTT→CTT)	ydiS	putative oxidoreductase with FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (37/27);  minor base A (1/1);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGTTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCG  >  NC_000913/1781641‑1781901                                                                                                                               |                                                                                                                                       gaaAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGt                                                                                                                           >  2:51791/1‑140 (MQ=255)    aGCGCCGTTACCCTCGATTTTCACCGCGAGCAACCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTcg                                                                                                                        >  1:264103/1‑140 (MQ=255)     gcgcCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTcgc                                                                                                                       >  2:768608/1‑140 (MQ=255)      cgcCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTcgcg                                                                                                                      >  1:794059/1‑140 (MQ=255)              cccTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGAcgcg                                                                                                              <  2:858407/140‑1 (MQ=255)                   gATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGt                                                                                                         >  2:503790/1‑140 (MQ=255)                        tCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCAGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTg                                                                                                    >  2:452333/1‑140 (MQ=255)                          aCCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGaa                                                                                                  <  2:424336/140‑1 (MQ=255)                             gcgAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGa                                                                                               >  2:72070/1‑140 (MQ=255)                             gcgAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGa                                                                                               >  2:1033940/1‑140 (MQ=255)                                agcagcCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAAc                                                                                            <  1:268550/140‑1 (MQ=255)                                agcagcCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCATGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAAc                                                                                            <  2:1006020/140‑1 (MQ=255)                                  cagcCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACaa                                                                                          >  2:531689/1‑140 (MQ=255)                                         aTGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACt                                                                                   >  2:840512/1‑140 (MQ=255)                                          tGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTg                                                                                  >  1:734507/1‑140 (MQ=255)                                               cGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTg                                                                             <  1:531689/140‑1 (MQ=255)                                               cGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTg                                                                             <  1:840512/140‑1 (MQ=255)                                                   acacGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAgg                                                                         <  1:363032/140‑1 (MQ=255)                                                     acGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGAAAACAAGGTCACTGGCGTGCAGGCt                                                                       >  1:591022/1‑140 (MQ=255)                                                      cGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCa                                                                                     <  1:363333/125‑1 (MQ=255)                                                      cGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCa                                                                                     >  2:363333/1‑125 (MQ=255)                                                         aTCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTgggg                                                                   >  2:435229/1‑140 (MQ=255)                                                          tCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGga                                                                  <  2:80612/140‑1 (MQ=255)                                                             tataCCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCAGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGgatga                                                               <  1:452333/140‑1 (MQ=255)                                                              ataCCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTggg                                                                    <  2:833603/134‑1 (MQ=255)                                                              ataCCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTggg                                                                    >  1:833603/1‑134 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTgg                                                                     >  2:247592/1‑126 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTgg                                                                     <  1:247592/126‑1 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGgatga                                                               <  1:662775/132‑1 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGgatga                                                               >  2:662775/1‑132 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGga                                                                  >  1:425209/1‑129 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGga                                                                  <  2:252433/129‑1 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGga                                                                  >  1:252433/1‑129 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGga                                                                  <  2:425209/129‑1 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCg                                                       >  2:808354/1‑140 (MQ=255)                                                                     aTTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCg                                                       >  2:1102367/1‑140 (MQ=255)                                                                      ttGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGa                                                      <  1:1102367/140‑1 (MQ=255)                                                                      ttGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGATTATCCCAGGCGTTCGCGTCGACGCGCTGGGTCGTGAAGGAAACAAGGTCACTGGCGTGGAGGGTGGGGGTGATAATCTCGa                                                      >  2:1143426/1‑140 (MQ=255)                                                                       tGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGaa                                                     >  2:803978/1‑140 (MQ=255)                                                                       tGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGAAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGaa                                                     <  2:591022/140‑1 (MQ=255)                                                                        gCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAg                                                    <  1:808354/140‑1 (MQ=255)                                                                        gCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAg                                                    <  1:435229/140‑1 (MQ=255)                                                                         cGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGc                                                   >  2:910988/1‑140 (MQ=255)                                                                           tAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGa                                                 >  2:670133/1‑140 (MQ=255)                                                                           tAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGa                                                 >  1:1118174/1‑140 (MQ=255)                                                                           tAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGa                                                 <  1:1033940/140‑1 (MQ=255)                                                                            aaTCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGaa                                                <  2:264103/140‑1 (MQ=255)                                                                              tCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATtct                                                         <  1:658117/129‑1 (MQ=255)                                                                              tCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATtct                                                         >  2:658117/1‑129 (MQ=255)                                                                                gTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGTGTGCAGGCTGGGGATGATATTCTCGAAGCGAAtgtg                                            <  1:503790/140‑1 (MQ=255)                                                                                    ggACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGaa                                                <  2:391957/132‑1 (MQ=255)                                                                                    ggACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGaa                                                >  1:391957/1‑132 (MQ=255)                                                                                          gTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGGTGATATTCTCGAAGCGAATGTGGTGATTCTgg                                  >  2:615170/1‑140 (MQ=255)                                                                                           tGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGc                                 >  2:986695/1‑140 (MQ=255)                                                                                           tGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGc                                 >  1:1061608/1‑140 (MQ=255)                                                                                           tGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGc                                 >  1:213899/1‑140 (MQ=255)                                                                                                gATGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATg                            >  1:882255/1‑140 (MQ=255)                                                                                                  tGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGc                          >  2:465288/1‑140 (MQ=255)                                                                                                  tGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGc                          >  2:459182/1‑140 (MQ=255)                                                                                                  tGGAACAAGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGc                          >  1:587137/1‑140 (MQ=255)                                                                                                        aaGCCGAGCAGGCGGGCGCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAAc                    >  1:583251/1‑140 (MQ=255)                                                                                                                    cgggcgCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGa                <  2:77382/132‑1 (MQ=255)                                                                                                                    cgggcgCGCAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGa                >  1:77382/1‑132 (MQ=255)                                                                                                                       gcgcgcAGATTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCcgc     <  1:1143426/140‑1 (MQ=255)                                                                                                                         gcgcAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCcgctc   <  1:72070/140‑1 (MQ=255)                                                                                                                          cgcAGCTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCcgctcg  <  2:1061608/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                       GAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCTGGACCCGTGGTTGATGGAACAAGCCGAGCAGGCGGGCGCGCAGTTTATCCCAGGCGTTCGCGTCGACGCGCTGGTTCGTGAAGGAAACAAGGTCACTGGCGTGCAGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCG  >  NC_000913/1781641‑1781901

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