Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,521,267	C→A	5.1%	L61L (CTG→CTT)	hofO ←	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,521,267	0	C	A	5.1%	281.9 / 4.9	80	L61L (CTG→CTT)	hofO	DNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base A (2/2);  total (41/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 = 3.40e-01

CGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGTG  >  NC_000913/3521120‑3521415                                                                                                                                                    |                                                                                                                                                     cGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCa                                                                                                                                                     >  3:622511/1‑149 (MQ=255)             ttCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTg                                                                                                                                         >  1:627505/1‑149 (MQ=255)                 cacaGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGc                                                                                                                                     >  3:102519/1‑149 (MQ=255)                  acaGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGTGAGCGGTGTCTACCAGGCGATACAGTGCGGcc                                                                                                                                    >  4:311514/1‑149 (MQ=255)                     gCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCAc                                                                                                                                 >  4:107038/1‑149 (MQ=255)                       gTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCActg                                                                                                                               <  3:124973/149‑1 (MQ=255)                        ttttCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCActgc                                                                                                                              >  3:99726/1‑149 (MQ=255)                            cAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGt                                                                                                                          <  1:592320/149‑1 (MQ=255)                            cAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGt                                                                                                                          <  4:75673/149‑1 (MQ=255)                             aaCGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTggaa                                                                                                                         >  2:48753/1‑147 (MQ=255)                              aCGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGtga                                                                                                                        <  4:173083/149‑1 (MQ=255)                                gCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGtgatg                                                                                                                      >  1:433961/1‑149 (MQ=255)                                gCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGtgatg                                                                                                                      >  3:12350/1‑149 (MQ=255)                                 ccAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGtgatga                                                                                                                     >  1:286832/1‑149 (MQ=255)                                   aaCTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACt                                                                                                                   >  3:408474/1‑149 (MQ=255)                                    aCTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACtt                                                                                                                  >  4:647220/1‑149 (MQ=255)                                       cGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGcc                                                                                                               >  3:166601/1‑149 (MQ=255)                                          cTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGt                                                                                                            >  3:348552/1‑149 (MQ=255)                                          cTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGt                                                                                                            >  3:676626/1‑149 (MQ=255)                                           tCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTa                                                                                                           <  3:647220/149‑1 (MQ=255)                                           tCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTa                                                                                                           >  2:59689/1‑149 (MQ=255)                                             ccTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGt                                                                                                         <  1:68238/149‑1 (MQ=255)                                             ccTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGt                                                                                                         <  4:102519/149‑1 (MQ=255)                                                  gcgGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAAt                                                                                                    >  1:655625/1‑149 (MQ=255)                                                        ggatgCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGc                                                                                              <  4:555046/149‑1 (MQ=255)                                                              cAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTcgcg                                                                                        >  3:564382/1‑149 (MQ=255)                                                               aGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCgcgc                                                                                       >  2:229323/1‑149 (MQ=255)                                                                ggAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGct                                                                                      >  2:331013/1‑149 (MQ=255)                                                                 gAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGctc                                                                                     <  3:15003/149‑1 (MQ=255)                                                                  aaaCCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGctct                                                                                    <  4:12350/149‑1 (MQ=255)                                                                  aaaCCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGctct                                                                                    >  3:616760/1‑149 (MQ=255)                                                                       aGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATg                                                                               <  1:573892/149‑1 (MQ=255)                                                                        gTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTAcc                                                                                                                                                      <  2:269777/77‑1 (MQ=255)                                                                        gTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTAcc                                                                                                                                                      >  1:269777/1‑77 (MQ=255)                                                                        gTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGccc                                                                                                              <  1:272316/117‑1 (MQ=255)                                                                        gTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGccc                                                                                                              >  2:272316/1‑117 (MQ=255)                                                                           gcgcgcCGGATAACTGAAAATCCAGTGGCGAAGAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGgtgt                                                                           <  4:564382/149‑1 (MQ=255)                                                                           gcgcgcCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCTGTAGTCGAATTAATGCGTCGCGCTCTTCATGgtgt                                                                           <  4:662010/149‑1 (MQ=255)                                                                             gcgcCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGgtgtgt                                                                         <  4:440954/149‑1 (MQ=255)                                                                              cgcCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTc                                                                        <  4:540188/149‑1 (MQ=255)                                                                              cgcCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTc                                                                        <  4:605979/149‑1 (MQ=255)                                                                                 cGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGAt                                                                     <  1:85574/149‑1 (MQ=255)                                                                                   gATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCTCGCTCTTCATGGTGTGTCGATGa                                                                   <  4:43388/149‑1 (MQ=255)                                                                                   gATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGa                                                                   <  2:624396/149‑1 (MQ=255)                                                                                    aTAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGAc                                                                  <  4:676626/149‑1 (MQ=255)                                                                                    aTAACTGAAAATCCAGTGGCAAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGAc                                                                  <  2:286832/149‑1 (MQ=255)                                                                                         tGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGaaa                                                             >  2:683197/1‑149 (MQ=255)                                                                                         tGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGaaa                                                             <  4:99726/149‑1 (MQ=255)                                                                                          gAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGaaaa                                                            >  2:385506/1‑149 (MQ=255)                                                                                              aTCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGa                                                        <  3:107038/149‑1 (MQ=255)                                                                                                      gCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaact                                                <  1:651392/149‑1 (MQ=255)                                                                                                        gAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGa                                                                   >  1:432934/1‑128 (MQ=255)                                                                                                        gAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGa                                                                   <  2:432934/128‑1 (MQ=255)                                                                                                        gAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaa                                                  <  1:561418/145‑1 (MQ=255)                                                                                                        gAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaa                                                  >  2:561418/1‑145 (MQ=255)                                                                                                         aaaaGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaactaac                                             <  1:683197/149‑1 (MQ=255)                                                                                                               gCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGc                                       >  2:421859/1‑149 (MQ=255)                                                                                                                 aGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCaa                                     >  1:307554/1‑149 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGa                                                        <  3:141066/129‑1 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGa                                                        >  4:141066/1‑129 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAAc                                    <  2:620253/149‑1 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAAc                                    >  2:6219/1‑149 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAAc                                    <  4:166601/149‑1 (MQ=255)                                                                                                                  gCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAAc                                    >  2:657503/1‑149 (MQ=255)                                                                                                                    gTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCa                                  >  3:194141/1‑149 (MQ=255)                                                                                                                            cTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGccc                          >  4:500907/1‑149 (MQ=255)                                                                                                                             tCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAAt                                                          <  1:32591/116‑1 (MQ=255)                                                                                                                             tCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAAt                                                          >  2:32591/1‑116 (MQ=255)                                                                                                                                 tGAAGGGAGCGGTGTCTACGAGGCGATACAGTGCGGCCCACTGCTGGGGATGGCTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAACAACCAGAATGCCCAGCaa                     >  2:490169/1‑149 (MQ=255)                                                                                                                                         gCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCg             >  1:82889/1‑149 (MQ=255)                                                                                                                                          cGGTGTCTACAAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCa                                      <  2:514223/123‑1 (MQ=255)                                                                                                                                          cGGTGTCTACAAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCa                                      >  1:1190/1‑123 (MQ=255)                                                                                                                                          cGGTGTCTACAAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCa                                      >  1:514223/1‑123 (MQ=255)                                                                                                                                          cGGTGTCTACAAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCa                                      <  2:1190/123‑1 (MQ=255)                                                                                                                                             tgtCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAgg         <  4:408474/149‑1 (MQ=255)                                                                                                                                              gtCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGc        <  4:622511/149‑1 (MQ=255)                                                                                                                                                cTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGc              >  2:118183/1‑141 (MQ=255)                                                                                                                                                cTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGc              <  1:118183/141‑1 (MQ=255)                                                                                                                                                  aCCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGTAAAGCTGGCGGAGGCgggg    >  4:471449/1‑149 (MQ=255)                                                                                                                                                  aCCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCgggg    <  2:627505/149‑1 (MQ=255)                                                                                                                                                  aCCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCgggg    <  1:421859/149‑1 (MQ=255)                                                                                                                                                  aCCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCgggg    >  3:130090/1‑149 (MQ=255)                                                                                                                                                   ccAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGt   <  3:471449/149‑1 (MQ=255)                                                                                                                                                   ccAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGt   <  3:500907/149‑1 (MQ=255)                                                                                                                                                    cAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGTg  >  3:428106/1‑149 (MQ=255)                                                                                                                                                    cAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGTg  >  1:406911/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGCAACCAGAATGCCCAGCAAAGCTGGCGGAGGCGGGGTG  >  NC_000913/3521120‑3521415

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