Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,686	A→C	N259H (AAC→CAC)	radA →	DNA repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,686	0	A	C	97.6%	303.8 / ‑1.3	83	N259H (AAC→CAC)	radA	DNA repair protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (46/35);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACG  >  NC_000913/4626540‑4626832                                                                                                                                                   |                                                                                                                                                   cgcgcgGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAGAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCg                                                                                                                                                >  1:547335/1‑151 (MQ=255)   cgcgGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGc                                                                                                                                               <  2:497761/150‑1 (MQ=255)   cgcgGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCt                                                                                                                                              >  2:1279008/1‑151 (MQ=255)    gcgGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCtt                                                                                                                                             <  2:978917/151‑1 (MQ=255)           ggCGATTGTCGTGGTGGGGCTCGTAACCAAAGACGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACGGTATCGACTGTGCGGTGCTTTTGAATGGCGATGCGGACTCCCGTTTCCGCACCTTGAGCAGCCATAAAAACCGCTTCGgcgc                                                                                                                                       <  2:461039/150‑1 (MQ=255)           ggCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGgcgc                                                                                                                                       <  2:174741/150‑1 (MQ=255)           ggCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGgcgc                                                                                                                                       <  1:1218200/150‑1 (MQ=255)               aTTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGa                                                                                                                                  <  1:706044/151‑1 (MQ=255)                 tGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGaa                                                                                                                                 <  2:680032/150‑1 (MQ=255)                  gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATg                                                                                                                               >  2:545812/1‑151 (MQ=255)                       ggtggGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTg                                                                                                                          >  1:870462/1‑151 (MQ=255)                        gtggGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTgg                                                                                                                         <  2:1106524/151‑1 (MQ=255)                            ggCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGt                                                                                                                     <  2:324081/151‑1 (MQ=255)                               aCGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCt                                                                                                                   <  1:1135358/150‑1 (MQ=255)                                 gTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                                                                                                 <  1:1030765/150‑1 (MQ=255)                                   aaCCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcgcg                                                                                                              <  2:1232650/151‑1 (MQ=255)                                       aaaGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATg                                                                                                           <  1:267792/150‑1 (MQ=255)                                        aaGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAc                                                                                                         <  2:1163459/151‑1 (MQ=255)                                         aGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                                                                                                        >  1:190676/1‑151 (MQ=255)                                              gTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAAc                                                                                                    >  2:20370/1‑150 (MQ=255)                                                 cGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAggg                                                                                                <  1:846375/151‑1 (MQ=255)                                                 cGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAgg                                                                                                 <  1:975218/150‑1 (MQ=255)                                                  gctggctgGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAACCACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAgggg                                                                                               >  1:684429/1‑151 (MQ=255)                                                    tggctggcCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCt                                                                                             >  2:655680/1‑151 (MQ=255)                                                    tggctggcCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCt                                                                                             >  2:1292193/1‑151 (MQ=255)                                                     ggctggcCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCCCCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTg                                                                                            <  2:321652/151‑1 (MQ=255)                                                        tggcCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGt                                                                                         >  2:115044/1‑151 (MQ=255)                                                            ccGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGa                                                                                       >  2:1296086/1‑149 (MQ=255)                                                            ccGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAg                                                                                     <  2:983421/151‑1 (MQ=255)                                                             cGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGt                                                                                    <  2:842374/151‑1 (MQ=255)                                                              gAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTc                                                                                   >  1:1234373/1‑151 (MQ=255)                                                              gAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTAGCGATGACCGAACAGGGGCTGCGTGAAGTc                                                                                   >  2:710339/1‑151 (MQ=255)                                                                  gTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                >  1:71863/1‑150 (MQ=255)                                                                  gTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                >  2:818990/1‑150 (MQ=255)                                                                    gCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAc                                                                             >  2:1148567/1‑151 (MQ=255)                                                                       ggAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAccc                                                                           <  2:113942/150‑1 (MQ=255)                                                                       ggAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAcc                                                                            <  1:1289044/149‑1 (MQ=255)                                                                       ggAACACGGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAc                                                                             <  2:706335/148‑1 (MQ=255)                                                                          acacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCg                                                                       >  2:330749/1‑151 (MQ=255)                                                                           cacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                        <  1:74279/149‑1 (MQ=255)                                                                            acTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGc                                                                     >  1:991390/1‑151 (MQ=255)                                                                            acTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGc                                                                     >  2:638839/1‑151 (MQ=255)                                                                              tGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCa                                                                    >  1:1290079/1‑150 (MQ=255)                                                                                  tCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtttt                                                               >  1:1140089/1‑151 (MQ=255)                                                                                       tGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTa                                                           >  2:1101582/1‑150 (MQ=255)                                                                                        gTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTaa                                                          <  1:294481/150‑1 (MQ=255)                                                                                        gTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTa                                                           >  1:47089/1‑149 (MQ=255)                                                                                          tCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGt                                                        <  2:500701/150‑1 (MQ=255)                                                                                          tCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                                                       >  1:803591/1‑151 (MQ=255)                                                                                           cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcg                                                      >  2:535834/1‑151 (MQ=255)                                                                                           cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcg                                                      <  2:1004574/151‑1 (MQ=255)                                                                                           cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                                                       <  2:174835/150‑1 (MQ=255)                                                                                            ggTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcgc                                                     <  2:1016246/151‑1 (MQ=255)                                                                                            ggTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcg                                                      <  1:446201/150‑1 (MQ=255)                                                                                              tGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgg                                                   >  2:100406/1‑151 (MQ=255)                                                                                                cTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcggcg                                                 >  1:28337/1‑151 (MQ=255)                                                                                                     ggATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGaa                                            <  1:242568/151‑1 (MQ=255)                                                                                                     ggATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGa                                             >  1:325697/1‑150 (MQ=255)                                                                                                     ggATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGa                                             >  1:843988/1‑150 (MQ=255)                                                                                                      gATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGa                                             >  1:57587/1‑149 (MQ=255)                                                                                                         ggCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTg                                         <  1:249710/150‑1 (MQ=255)                                                                                                            gATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCt                                     >  1:968104/1‑151 (MQ=255)                                                                                                             aTGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTc                                    >  1:426202/1‑151 (MQ=255)                                                                                                              tGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTc                                    >  1:139570/1‑150 (MQ=255)                                                                                                              tGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTcc                                   >  1:881381/1‑151 (MQ=255)                                                                                                              tGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTcc                                   >  1:280452/1‑151 (MQ=255)                                                                                                                ccGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCgg                                 <  2:392866/151‑1 (MQ=255)                                                                                                                 cGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGCGCTAGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCg                                  >  1:771650/1‑149 (MQ=255)                                                                                                                 cGACTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGAGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGt                                >  2:780151/1‑151 (MQ=255)                                                                                                                   aCTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGt                                >  2:364257/1‑149 (MQ=255)                                                                                                                    cTCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGc                             >  2:426399/1‑151 (MQ=255)                                                                                                                     tCCCGTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCt                            >  2:143889/1‑151 (MQ=255)                                                                                                                         gTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAg                         >  1:1205879/1‑150 (MQ=255)                                                                                                                         gTTTTCGCACCTTGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGt                        <  1:1025852/151‑1 (MQ=255)                                                                                                                                     tGCGCAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATg            >  1:34446/1‑151 (MQ=255)                                                                                                                                       cgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATggg          >  1:24831/1‑151 (MQ=255)                                                                                                                                        gcAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATggg          >  1:1027760/1‑150 (MQ=255)                                                                                                                                        gcAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGgga         <  2:991450/151‑1 (MQ=255)                                                                                                                                         cAGCCATAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGgga         >  1:824951/1‑150 (MQ=255)                                                                                                                                              aTAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAAc   <  2:1191627/151‑1 (MQ=255)                                                                                                                                              aTAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAAc   >  2:1118059/1‑151 (MQ=255)                                                                                                                                              aTAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAAc   <  1:1278930/151‑1 (MQ=255)                                                                                                                                               tAAACACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTTGTATGGGAAGGAAcg  >  2:1163108/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACG  >  NC_000913/4626540‑4626832

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