Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	922,723	G→A	intergenic (‑133/‑190)	cspD ← / → clpS	inhibitor of DNA replication, cold shock protein homolog/regulatory protein for ClpA substrate specificity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	922,723	0	G	A	98.9%	367.9 / ‑4.2	94	intergenic (‑133/‑190)	cspD/clpS	inhibitor of DNA replication, cold shock protein homolog/regulatory protein for ClpA substrate specificity
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (49/44);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCCTTTTCCATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCA  >  NC_000913/922580‑922870                                                                                                                                                |                                                                                                                                                    cccTTTTCCATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACaa                                                                                                                                              <  2:742117/151‑1 (MQ=255) cccTTTTCCATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACaa                                                                                                                                              <  2:1005116/151‑1 (MQ=255)     tttCCATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGacaagggacaaggg                                                                                                                                           <  1:1388898/150‑1 (MQ=255)        ccATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCa                                                                                                                                        <  2:95871/150‑1 (MQ=255)           tGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGt                                                                                                                                     >  2:277661/1‑150 (MQ=255)           tGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGt                                                                                                                                     >  1:72215/1‑150 (MQ=255)            gCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTAt                                                                                                                                   >  1:679892/1‑151 (MQ=255)             cTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTAt                                                                                                                                   <  1:142600/150‑1 (MQ=255)               tCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGa                                                                                                                                 <  1:1128620/150‑1 (MQ=255)                  aCATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATc                                                                                                                              >  1:943967/1‑150 (MQ=255)                  aCATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATc                                                                                                                              >  1:417891/1‑150 (MQ=255)                    aTCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCaaa                                                                                                                           <  2:20798/151‑1 (MQ=255)                    aTCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCaaa                                                                                                                           <  1:523833/151‑1 (MQ=255)                      ccTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCaaaa                                                                                                                          >  2:1355230/1‑150 (MQ=255)                      ccTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCaaaa                                                                                                                          >  1:1466116/1‑150 (MQ=255)                        ttCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAAtt                                                                                                                       >  2:904378/1‑151 (MQ=255)                          cGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAAttt                                                                                                                      <  2:353855/150‑1 (MQ=255)                          cGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTg                                                                                                                     >  1:1453175/1‑151 (MQ=255)                              aaTCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGc                                                                                                                 <  1:673367/151‑1 (MQ=255)                                   tataCAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAg                                                                                                               <  2:956434/148‑1 (MQ=255)                                    ataCAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTaa                                                                                                           <  2:1283992/151‑1 (MQ=255)                                         aGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCt                                                                                                       >  1:1185024/1‑150 (MQ=255)                                           tAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTAtt                                                                                                    >  1:441927/1‑151 (MQ=255)                                               aTGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGAc                                                                                                 >  1:620846/1‑150 (MQ=255)                                               aTGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACa                                                                                                <  1:799692/151‑1 (MQ=255)                                                  gAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAgg                                                                                              >  2:464811/1‑150 (MQ=255)                                                   aaTAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGaa                                                                                            <  1:534299/151‑1 (MQ=255)                                                     tAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtg                                                                                          >  1:306202/1‑151 (MQ=255)                                                       aaCCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtgtg                                                                                        >  2:617235/1‑151 (MQ=255)                                                       aaCCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtgtg                                                                                        >  2:677360/1‑151 (MQ=255)                                                       aaCCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtgtg                                                                                        >  2:386409/1‑151 (MQ=255)                                                         ccGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCATCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGAc                                                                                      <  2:746187/151‑1 (MQ=255)                                                             ggTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGAt                                                                                  <  2:407898/151‑1 (MQ=255)                                                              gTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATg                                                                                 >  1:1116619/1‑151 (MQ=255)                                                                 agagagGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCg                                                                              <  2:1176606/151‑1 (MQ=255)                                                                  gagagGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGc                                                                             >  2:1337838/1‑151 (MQ=255)                                                                    gagGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTg                                                                           >  1:316774/1‑151 (MQ=255)                                                                          gCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCa                                                                     <  1:388974/151‑1 (MQ=255)                                                                             gTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCaa                                                                    <  2:382445/149‑1 (MQ=255)                                                                               tCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa                                                                 >  2:172923/1‑150 (MQ=255)                                                                               tCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGAt                                                                >  1:513626/1‑151 (MQ=255)                                                                                cAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGAt                                                                >  1:1103202/1‑150 (MQ=255)                                                                                cAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGAt                                                                >  1:667615/1‑150 (MQ=255)                                                                                cAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATa                                                               <  1:390921/151‑1 (MQ=255)                                                                                   aaCCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGAt                                                                >  2:498096/1‑147 (MQ=255)                                                                                    aCCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatg                                                           >  2:1465533/1‑151 (MQ=255)                                                                                    aCCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatg                                                           <  1:1392229/151‑1 (MQ=255)                                                                                    aCCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatg                                                           >  1:1331025/1‑151 (MQ=255)                                                                                     ccTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatg                                                           >  2:1300239/1‑150 (MQ=255)                                                                                      cTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatgat                                                         >  1:1410278/1‑151 (MQ=255)                                                                                        cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATa                                                        <  1:865881/150‑1 (MQ=255)                                                                                        cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAg                                                       <  2:543417/151‑1 (MQ=255)                                                                                         gCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGt                                                      >  1:835302/1‑151 (MQ=255)                                                                                          ccAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGtt                                                     <  2:1041680/151‑1 (MQ=255)                                                                                            aaCTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtat                                                   <  2:1329516/151‑1 (MQ=255)                                                                                                ctAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatct                                                 <  1:911821/149‑1 (MQ=255)                                                                                                ctAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatct                                                 <  2:684752/149‑1 (MQ=255)                                                                                                 tAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatctatc                                              >  1:154523/1‑151 (MQ=255)                                                                                                 tAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatctatc                                              >  2:754347/1‑151 (MQ=255)                                                                                                  aGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatcta                                                <  2:786300/148‑1 (MQ=255)                                                                                                   gAAATACAAGTTAGAGAATTAGGGCGAGTCGGCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGATATTGACAGCAATGTGACAGATGTCGCTGATGCCAACGACAGATGATAGTtattta                                                <  2:821700/147‑1 (MQ=255)                                                                                                     aaTACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCAtg                                           >  2:270542/1‑150 (MQ=255)                                                                                                        aCAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTgg                                        <  2:1451464/150‑1 (MQ=255)                                                                                                        aCAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGa                                       <  2:815541/151‑1 (MQ=255)                                                                                                            tatAGAGAATTAGTGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGAGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  2:1431273/3‑149 (MQ=255)                                                                                                             ttAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  2:419519/1‑148 (MQ=255)                                                                                                             ttAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAg                                   <  2:445750/150‑1 (MQ=255)                                                                                                              tAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  1:524753/1‑147 (MQ=255)                                                                                                              tAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  2:244835/1‑147 (MQ=255)                                                                                                              tAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  2:1361767/1‑147 (MQ=255)                                                                                                              tAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGt                                     >  1:1434313/1‑147 (MQ=255)                                                                                                               agagAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGAtt                                <  1:464865/151‑1 (MQ=255)                                                                                                               agagAATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGAt                                 <  1:1195917/150‑1 (MQ=255)                                                                                                                 agaATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTgg                              >  1:1274979/1‑151 (MQ=255)                                                                                                                 agaATTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTgg                              <  1:457184/151‑1 (MQ=255)                                                                                                                   aaTTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTc                            >  1:810299/1‑151 (MQ=255)                                                                                                                    aTTAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTc                            >  1:1132323/1‑150 (MQ=255)                                                                                                                     ttAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGt                             >  1:643559/1‑148 (MQ=255)                                                                                                                     ttAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTc                            >  2:1354113/1‑149 (MQ=255)                                                                                                                      tAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAg                          <  1:32236/150‑1 (MQ=255)                                                                                                                      tAGGGCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAg                          <  1:811876/150‑1 (MQ=255)                                                                                                                         ggCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCaa                      <  2:310623/151‑1 (MQ=255)                                                                                                                          gCGAGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCaaa                     <  1:617308/151‑1 (MQ=255)                                                                                                                          ccgacgcGTCAAGCATTTGATTAGGGTCAAGGGGCAGGTGTGAATCAAAAATTTGAACCAGTTAACGCTATTGACAGGGATGTGTCAGATGTCGCTGATGCGAACGATAGATGATCGTTATCTATCATGTGGAGTAGATTGGTCAGGCa                       >  2:283395/7‑149 (MQ=255)                                                                                                                             aGCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATa                   >  1:997865/1‑150 (MQ=255)                                                                                                                              gCCGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATa                   >  2:17265/1‑149 (MQ=255)                                                                                                                               ccGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGc                <  1:591979/151‑1 (MQ=255)                                                                                                                                cGTCAAGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGct               >  2:1195507/1‑151 (MQ=255)                                                                                                                                     aGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt          <  1:498154/151‑1 (MQ=255)                                                                                                                                     aGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt          <  2:679818/151‑1 (MQ=255)                                                                                                                                     aGCATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt          <  2:1466008/151‑1 (MQ=255)                                                                                                                                       cATTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGta         <  2:84163/150‑2 (MQ=255)                                                                                                                                        aTTTGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTc         >  2:534184/1‑149 (MQ=255)                                                                                                                                          ttGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGc      <  2:810820/150‑1 (MQ=255)                                                                                                                                           tGACAAGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCg     <  1:904463/150‑1 (MQ=255)                                                                                                                                            gacaagggacaaggGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGc   <  2:1116520/151‑1 (MQ=255)                                                                                                                                             acaagggacaagggGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCa  <  1:465786/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CCCTTTTCCATGCTTCGACATCCTTCGCAAATCTTATACAAGTAAGATGGAATAAACCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCA  >  NC_000913/922580‑922870

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