Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,323,286	C→T	16.5%	intergenic (+35/‑161)	atoC → / → atoD	fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,323,286	0	C	T	16.5%	232.3 / 33.4	97	intergenic (+35/‑161)	atoC/atoD	fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (48/33);  new base T (10/6);  total (58/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.11e-01

CATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAT  >  NC_000913/2323188‑2323378                                                                                                   |                                                                                             cATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcg                                                                                             >  2:17644/1‑100 (MQ=255)    cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCaa                                                                                          >  1:79990/1‑100 (MQ=255)      gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                        >  2:1166490/1‑100 (MQ=255)        cGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                      >  2:1478554/1‑100 (MQ=255)         gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                     <  1:549509/100‑1 (MQ=255)         gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                     >  1:360401/1‑100 (MQ=255)          ccGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCt                                                                                    >  2:467546/1‑100 (MQ=255)           cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                   <  2:2085436/100‑1 (MQ=255)           cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                   <  2:578972/100‑1 (MQ=255)           cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                   <  1:535708/100‑1 (MQ=255)             tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCa                                                                                 <  1:1848396/100‑1 (MQ=25)             tGCATTGATATATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCa                                                                                 <  1:1068818/100‑1 (MQ=25)              gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCAt                                                                                >  1:1799845/1‑100 (MQ=255)               cATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATa                                                                               <  1:2104985/100‑1 (MQ=25)                aTTGATGTATAAACTCCAGGAATACGGTATCTATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                              <  1:451800/100‑1 (MQ=25)                aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                              <  2:366917/100‑1 (MQ=25)                 ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                             >  1:1877163/1‑100 (MQ=255)                 ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                             >  2:2521241/1‑100 (MQ=255)                   gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                           >  2:802084/1‑100 (MQ=255)                    aTGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGcaag                                                                          <  2:401866/100‑4 (MQ=21)                     tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                         <  2:595220/100‑1 (MQ=255)                      gTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTc                                                                        >  2:2341162/1‑100 (MQ=255)                        ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                      <  1:1450447/100‑1 (MQ=255)                        ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                      >  2:1229916/1‑100 (MQ=255)                         taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAtt                                                                     >  2:300239/1‑100 (MQ=255)                         taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAtt                                                                     >  2:2514787/1‑100 (MQ=255)                          aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCATAACTTGTTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCGTtc                                                                    >  1:2422502/1‑100 (MQ=255)                          aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                    <  2:1129607/100‑1 (MQ=255)                           aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                   >  1:1537810/1‑100 (MQ=255)                             cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                 <  1:694061/100‑1 (MQ=255)                             cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                 >  2:828180/1‑100 (MQ=255)                              tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                <  1:1448145/100‑1 (MQ=255)                              tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                <  2:96556/100‑1 (MQ=255)                               ccAGGAATACGGTATCGATCCTGCGGATGTATAACACCACAACCTGCTATGCAGTAATTTGCACAGTGCGCAATTGTCTGCATAGCCGCTc                                                                        >  1:1674854/1‑91 (MQ=255)                               ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                               <  1:300236/100‑1 (MQ=255)                                    aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaaa                                                          >  2:288233/1‑100 (MQ=255)                                     aTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAAt                                                         >  2:1216361/1‑100 (MQ=255)                                       aCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcc                                                       <  2:463563/100‑1 (MQ=255)                                        ccgTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                      >  2:1751027/3‑100 (MQ=255)                                        cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                      <  1:153065/100‑1 (MQ=255)                                         ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcccc                                                     >  1:1775594/1‑100 (MQ=255)                                           tATCGATCCGTCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCAt                                                   <  2:2012075/100‑1 (MQ=255)                                           tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCAt                                                   >  2:1693377/1‑100 (MQ=255)                                                  ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                            >  2:1453211/1‑100 (MQ=255)                                                    ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                          <  2:2542325/100‑1 (MQ=255)                                                     gcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATc                                         <  1:1807798/100‑1 (MQ=255)                                                     gcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATc                                         <  2:1722771/100‑1 (MQ=255)                                                      cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                        >  2:2031727/1‑100 (MQ=255)                                                       ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                       <  2:867865/100‑1 (MQ=255)                                                       ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                       <  2:358602/100‑1 (MQ=255)                                                         aTGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCtt                                     >  2:710636/1‑100 (MQ=255)                                                            tataACACCAAAACTTGCTATTCAGAAATTTGCACAGTACGCAATTTTCTGCATAGACGCTCATTCTCCTTATAAATCCCAATCCAATTTATcctt                                      >  1:2104465/1‑94 (MQ=255)                                                            tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTACCCCTTCat                                  >  2:2398340/1‑100 (MQ=255)                                                             ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCata                                 >  2:2154008/1‑100 (MQ=255)                                                                acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                              <  2:2142496/100‑1 (MQ=255)                                                                acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                              >  2:2480236/1‑100 (MQ=255)                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                             >  2:687505/1‑100 (MQ=255)                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                             >  1:112396/1‑100 (MQ=255)                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATCTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                             >  1:583360/1‑100 (MQ=255)                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                            >  1:1683835/1‑100 (MQ=255)                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                            >  1:2014009/1‑100 (MQ=255)                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                            >  2:623312/1‑100 (MQ=255)                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                            >  1:879707/1‑100 (MQ=255)                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                           >  2:1560804/1‑100 (MQ=255)                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                           <  2:575744/100‑1 (MQ=255)                                                                     aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                           <  1:1714278/98‑1 (MQ=255)                                                                     aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTa                         <  1:1096250/100‑1 (MQ=255)                                                                      aaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                        <  1:1447344/100‑1 (MQ=255)                                                                      aaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                        <  1:2480217/100‑1 (MQ=255)                                                                       agcTTGTTATGCAAAAATTTGCACAGTACGCTATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCATGCATATTCAAt                           <  1:974275/94‑1 (MQ=255)                                                                       aaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCGTAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                        <  1:1817038/99‑1 (MQ=255)                                                                        aCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                      <  2:318204/100‑1 (MQ=255)                                                                         cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTTATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                     <  2:644145/100‑1 (MQ=255)                                                                          ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTCATCCCTTCATATTCAATTAGTTaa                    >  1:1790322/1‑100 (MQ=255)                                                                            gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                  <  1:1166486/100‑1 (MQ=255)                                                                            gCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATGAAACCCCATACAATTTATCCCTTCATATTCAATTAGTTAAAt                  >  1:1814622/1‑100 (MQ=21)                                                                             cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                 >  2:757696/1‑100 (MQ=255)                                                                                tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAACTAGTTAAATAACt              >  2:1503934/1‑100 (MQ=25)                                                                                 gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa             >  1:2422026/1‑100 (MQ=255)                                                                                 gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa             >  2:2043001/1‑100 (MQ=255)                                                                                 gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa             <  1:1649998/100‑1 (MQ=255)                                                                                 gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa             >  2:2120648/1‑100 (MQ=25)                                                                                 gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa             >  1:35934/1‑100 (MQ=25)                                                                                  cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa            >  2:1876336/1‑100 (MQ=255)                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt              >  1:2159364/1‑98 (MQ=37)                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa            >  2:72641/1‑100 (MQ=37)                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa            >  2:2426894/1‑100 (MQ=37)                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa            >  2:1499141/1‑100 (MQ=37)                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa            >  1:574620/1‑100 (MQ=37)                                                                                    gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt          >  2:1074261/1‑100 (MQ=255)                                                                                    gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt          >  2:566612/1‑100 (MQ=255)                                                                                    gAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt          >  1:2372057/1‑100 (MQ=37)                                                                                     aaaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc         <  1:687502/100‑1 (MQ=255)                                                                                      aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc        >  1:2243170/1‑100 (MQ=255)                                                                                       aTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa       >  2:1382150/1‑100 (MQ=255)                                                                                        tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa      >  2:1942099/1‑100 (MQ=255)                                                                                         ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat     <  2:723638/100‑1 (MQ=255)                                                                                             acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat  >  2:1512827/1‑99 (MQ=255)                                                                                                   |                                                                                             CATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAT  >  NC_000913/2323188‑2323378

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