Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,518,082	C→G	8.1%	noncoding (35/76 nt)	alaX ←	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,518,082	0	C	G	8.1%	191.6 / 7.5	62	noncoding (35/76 nt)	alaX	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (21/36);  new base G (2/3);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTAT  >  NC_000913/2517933‑2518225                                                                                                                                                      |                                                                                                                                                tCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgc                                                                                                                                                 >  1:896534/1‑150 (MQ=255) tCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgc                                                                                                                                                 >  2:216523/1‑150 (MQ=255)      tgatgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatgc                                                                                                                                            <  2:521077/150‑1 (MQ=255)       gatgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCa                                                                                                                                           <  2:359411/150‑1 (MQ=255)       gatgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCCGGATCGAACCGCTGACTTCTTGCATGCCATGCa                                                                                                                                           <  1:889862/150‑1 (MQ=255)           cTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                       >  1:1249801/1‑150 (MQ=255)            tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg                                                                                                                                      <  1:940402/150‑4 (MQ=255)               cACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct                                                                                                                                   >  2:202591/1‑150 (MQ=255)                aCCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCCAACCGCTGACCTCTTGCATGCCATGCAAGCGCtctc                                                                                                                                  <  1:855125/150‑1 (MQ=255)                aCCAGGCGGCACTAACCCGCACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGACCCGCTGACCTCTTGCATGGCATGCAAgaggtcag                                                                                                                                  <  1:1616090/150‑8 (MQ=38)                 ccAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg                                                                                                                                         <  2:645961/142‑1 (MQ=255)                 ccAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg                                                                                                                                         >  1:645961/1‑142 (MQ=255)                 ccAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCGCTTGCATGGCATGCAAgaggtc                                                                                                                                    >  1:1524544/1‑142 (MQ=255)                 ccAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCGCTTGCATGGCATGCAAgaggtc                                                                                                                                    <  2:1524544/147‑6 (MQ=255)                   aGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCa                                                                                                                               <  2:382435/150‑1 (MQ=255)                       ggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTg                                                                                                                           <  1:719777/150‑1 (MQ=255)                        gcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                          >  2:465248/1‑150 (MQ=255)                        gcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                          >  2:1298051/1‑150 (MQ=255)                         cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGCCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                         >  1:722108/1‑150 (MQ=255)                         cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                         <  2:595200/150‑1 (MQ=255)                         cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                         <  1:1298051/150‑1 (MQ=255)                          aCTAACCCGGCCGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAACACCCCCGCGTTGCAGCGAGGGTTGTAAATTTTGTGGAGCTAATCGTGATCGAAACGCTGACCTCTTGCATGCCAAGCATGCGCTCTCCCAGCTGAGc                                                                                                                        <  2:1433267/150‑1 (MQ=25)                           cTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGCCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCt                                                                                                                       >  1:693737/1‑150 (MQ=255)                                 cGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcc                                                                                                                 <  2:896534/150‑1 (MQ=255)                                    aCGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCa                                                                                                              <  1:807578/150‑1 (MQ=255)                                     cGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGccccgt                                                                                                             >  1:1367366/1‑148 (MQ=255)                                     cGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGccccga                                                                                                             >  2:294462/1‑148 (MQ=255)                                        aGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGAt                                                                                                          >  2:1124097/1‑150 (MQ=255)                                         gACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATACAAAGCAAGAGCACACCAATCTGATCAATACACCCAAGaag                                                                                                         >  1:1100991/1‑148 (MQ=25)                                          aCTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGc                                                                                                        >  1:790365/1‑150 (MQ=255)                                               aTTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGCACCGCTGCCCTCTTGCAGGCCAGGCACGCTCTCTCCCCGGTGGGCGATAGCCCCACGATGAGTTTa                                                                                                   >  1:1125397/1‑150 (MQ=25)                                               aTTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTa                                                                                                   <  1:790942/150‑1 (MQ=255)                                               aTTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTa                                                                                                   <  2:1228592/150‑1 (MQ=255)                                                 tACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCCGGAGGGGTTTagt                                                                                                 >  2:231754/1‑148 (MQ=255)                                                 tACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACg                                                                                                 <  1:1145255/150‑1 (MQ=255)                                                   cAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTa                                                                                               <  2:1420145/150‑1 (MQ=255)                                                    aaaaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                               <  1:80278/101‑1 (MQ=255)                                                    aaaaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                               >  2:80278/1‑101 (MQ=255)                                                      aaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCa                                                                                            >  2:955439/1‑150 (MQ=255)                                                            cAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTg                                                                                      >  2:71345/1‑150 (MQ=255)                                                               aaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTg                                                                                   >  1:726210/1‑150 (MQ=255)                                                                 aGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTggg                                                                                 <  1:71345/150‑1 (MQ=255)                                                                           aaaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAttt                                                                       <  1:771582/150‑1 (MQ=255)                                                                            aaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtg                                                                      <  2:722391/150‑1 (MQ=255)                                                                              aaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                    <  2:959377/150‑1 (MQ=255)                                                                              aaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                    <  1:216523/150‑1 (MQ=255)                                                                              aaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                    <  1:955439/150‑1 (MQ=255)                                                                              aaCCCCCGCGTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGACCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                    <  2:726068/150‑1 (MQ=255)                                                                               aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                                   <  1:1297459/150‑1 (MQ=255)                                                                               aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                                   <  1:1567226/150‑1 (MQ=255)                                                                                cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtgg                                                                  <  1:1124097/150‑1 (MQ=255)                                                                                cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtgg                                                                  <  2:1395904/150‑1 (MQ=255)                                                                                cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtgg                                                                  <  2:1240967/150‑1 (MQ=255)                                                                                cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtgg                                                                  <  1:465248/150‑1 (MQ=255)                                                                                  cccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                                <  2:605259/150‑1 (MQ=255)                                                                                   ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGATCTATAGCGCCCCCGTCGGATAACGAGGCAAGGTTGCAGGGGGAAAAATTTGGTGGAgg                                                               >  1:774585/1‑149 (MQ=17)                                                                                   ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGc                                                               <  2:726210/150‑1 (MQ=255)                                                                                    cGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCCTGCAAGCGCTCTCCCAGCGGAGCTATAGCCCCACGATGCGCTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCt                                                              <  2:1531175/150‑1 (MQ=255)                                                                                              ctaGGGTTGGAAATTTTGTGTATCTAATCGTGATCGAACCGCTTAACTCTTTCATTCCATGCAAGCGCTTTCCCATATGAGCTATATCCCCACTATGCGTTTACGTACCAAGTTTGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATc                                                    <  2:790365/148‑1 (MQ=255)                                                                                                       gAAATTTGGTGGAGCTAAGCGGGATCTAACCTCTGACCTCTTTCATGCCATGCAAGCGCTCTCCCATCTGAGCTATATCCCCACTATTAGTTTACGTACCAAGTTTGCTGGGTTGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTg                                           <  2:1249801/150‑1 (MQ=21)                                                                                                                      tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGa                                                      >  2:1271781/1‑124 (MQ=255)                                                                                                                      tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGa                                                      <  1:1271781/124‑1 (MQ=255)                                                                                                                      tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGAGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGAGGGATCGAACCGCTGACCTCTTGcatgcca                            >  2:1291278/1‑150 (MQ=255)                                                                                                                                       gagagcgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAg           >  2:1240714/8‑150 (MQ=255)                                                                                                                                        tGACCTCTACCATGCACGGCAAGCGTTCTCCAAGCTGTGCAATACGCCGACCAAGACTTTACGGACCAAGTATGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGc          <  1:321516/150‑1 (MQ=17)                                                                                                                                             cctTGCATGCCATGCAAGCGCTCTCCCATCTGAGCTATATCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGc     <  1:1240647/149‑1 (MQ=255)                                                                                                                                                tGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCAAGATGCGTTTACGTACCAAGTTTGCTGCGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAACTCTTGCATGCCATGCAAGCGCTCTACCAGCTGAGCtat  >  2:1463966/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                TCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTAT  >  NC_000913/2517933‑2518225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑