Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,272	T→A	16.1%	pseudogene (716/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,272	0	T	A	16.1%	208.7 / 30.6	87	pseudogene (716/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (8/6);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGG  >  NC_000913/1597174‑1597368                                                                                                   |                                                                                                 tatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTACGATACCGGTGGCAtt                                                                                                 >  2:2390956/1‑100 (MQ=255)     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                             >  2:530363/1‑100 (MQ=255)     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                             >  1:340581/1‑100 (MQ=255)     gACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                             <  2:2159674/100‑1 (MQ=255)        tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                          <  2:1108249/100‑1 (MQ=255)        tattgactgtgccattaccgaaatatccattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                          >  1:608705‑M2/63‑100 (MQ=255)        tattgactgtgccattaccgaaatatccattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                          >  1:924092‑M2/63‑100 (MQ=255)          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGTTCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                        <  1:925688/100‑3 (MQ=255)          tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGAt                                                                                        >  2:455816/1‑100 (MQ=255)             cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatn                                                                                     <  1:577121/100‑3 (MQ=255)             cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                     <  2:1437424/100‑3 (MQ=255)             actgtgccattaccgaaatatccattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGgtatc                                                                                     <  2:1345344‑M2/43‑3 (MQ=255)               ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                     <  1:1172623/98‑3 (MQ=255)               tcaCCGAGATATCCATTATAGAGTCTGCTTACCGTGACGAGAGGCGCTTGCTCCTGCACATTCAGCGTTCCGATACCGGTGTCAATGCCACCGGTatcg                                                                                    <  1:2418330/97‑4 (MQ=255)                    gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                              <  1:115995/100‑1 (MQ=255)                              ttATAGGGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCGGCCGTa                                                                     <  1:1988995/99‑1 (MQ=255)                               tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                   <  2:2202334/100‑1 (MQ=255)                                 tccattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca                                                                 <  1:1787952‑M2/63‑1 (MQ=255)                                  ccattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                >  1:2257228‑M2/37‑100 (MQ=255)                                  agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                >  1:2118084/1‑100 (MQ=255)                                   gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaac                                                               <  1:1939686/100‑1 (MQ=255)                                   cattatagagtctgcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaac                                                               >  1:1142594‑M2/36‑100 (MQ=255)                                     nccTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccc                                                             >  2:2398172/3‑100 (MQ=255)                                      tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccc                                                            <  2:590308/100‑1 (MQ=255)                                        tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccccc                                                          <  2:1167642/100‑1 (MQ=255)                                         gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccccc                                                         >  2:10777/1‑100 (MQ=255)                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCt                                                        <  2:1395763/100‑1 (MQ=255)                                            tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcg                                                      >  2:460635/1‑100 (MQ=255)                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                    >  2:1493390/1‑100 (MQ=255)                                               cGTGATGCCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                   >  2:1362939/1‑100 (MQ=255)                                                ntgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                  >  1:403715/2‑100 (MQ=255)                                                gcgtaccgtgatgacagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                  <  2:753451‑M2/78‑1 (MQ=255)                                                gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                  >  1:720506/1‑100 (MQ=255)                                                    tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTGTCAGCCGTAACCAACCCCCCCTCGCGACTAGt                                              >  2:41340/1‑100 (MQ=255)                                                    accgtggtgacagggtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCGCCGATAATGGTATTTTCAGCCGTAACCAGCCCCCCCTCGCGAATAg                                               <  1:474085‑M2/81‑1 (MQ=255)                                                     gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTc                                             <  1:13816/100‑1 (MQ=255)                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                           >  1:633936/1‑100 (MQ=255)                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                           >  2:223953/1‑100 (MQ=255)                                                          agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTc                                        <  2:821229/100‑1 (MQ=255)                                                           gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                       >  2:64780/1‑100 (MQ=255)                                                              cagagtctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCa                                    >  2:1343622‑M2/9‑100 (MQ=255)                                                                   tctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGCAACCAA‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                              <  1:1178414‑M2/97‑1 (MQ=255)                                                                   ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                               >  2:755037/1‑100 (MQ=255)                                                                    ctGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACTCACCCTCGCGAATAGTCGCCTCCCCAGTTCCt                              >  1:1308874‑M2/3‑100 (MQ=255)                                                                     tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                             >  2:1135862‑M2/2‑100 (MQ=255)                                                                     tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                              <  1:739261‑M2/98‑1 (MQ=255)                                                                      gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                            <  2:586733/100‑1 (MQ=255)                                                                      gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                            >  2:1459904/1‑100 (MQ=255)                                                                      gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAACAGTCGCCTCCCCAGTTcc                               >  1:2247086/1‑97 (MQ=255)                                                                       cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                           <  1:410578/100‑1 (MQ=255)                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                         >  1:2222194/1‑100 (MQ=255)                                                                         caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                         <  1:1099688/100‑1 (MQ=255)                                                                          aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                        >  1:866889/1‑100 (MQ=255)                                                                           ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                       >  2:1187078/1‑100 (MQ=255)                                                                           ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                       <  2:403528/100‑1 (MQ=255)                                                                           ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATGTTCAGCCGTGAGCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTGGATTAc                       >  1:894198/1‑100 (MQ=255)                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                      <  2:1466457/100‑1 (MQ=255)                                                                              aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                    <  2:309734/100‑1 (MQ=255)                                                                              aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCCTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                    <  1:282003/100‑1 (MQ=255)                                                                                ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                  >  1:2150478/1‑100 (MQ=255)                                                                                ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                  <  2:689567/100‑1 (MQ=255)                                                                                   tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa               <  1:425322/100‑1 (MQ=255)                                                                                    cgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa              >  2:2327076/8‑100 (MQ=255)                                                                                    cgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa              >  2:325906/8‑100 (MQ=255)                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa               >  2:812603/1‑99 (MQ=255)                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCGCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa              >  1:1812929/1‑100 (MQ=255)                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt            <  2:258023/100‑1 (MQ=255)                                                                                       ttatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCATCCCCAGTTCCTTGATTACCAATTTGAAAtt           >  1:451710/5‑100 (MQ=255)                                                                                       ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt           >  2:1985589/5‑100 (MQ=255)                                                                                        tACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc          <  2:1898006/100‑1 (MQ=255)                                                                                         atcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa         >  1:1421387/3‑100 (MQ=255)                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa         >  1:149611/3‑100 (MQ=255)                                                                                         atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa         >  2:2090969/3‑100 (MQ=255)                                                                                         aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa         <  1:2171058/100‑1 (MQ=255)                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCTAAGAGTCGCCTCCCCAGTTCCTTTATTACCAATTTGAAATTCaa        <  1:1251454/99‑1 (MQ=255)                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt       <  1:460633/100‑1 (MQ=255)                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt       >  2:2005079/1‑100 (MQ=255)                                                                                            ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt      <  2:2175084/100‑1 (MQ=255)                                                                                            ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt      >  1:508055/1‑100 (MQ=255)                                                                                            ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt      >  1:1621204/1‑100 (MQ=255)                                                                                             gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt       >  1:419312/1‑98 (MQ=255)                                                                                             gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg     <  2:924095/100‑1 (MQ=255)                                                                                             gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg     <  2:340583/100‑1 (MQ=255)                                                                                             gTGGCATTGCCACCGATAATCGTATTTTCAGCCATAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg     <  1:1184525/100‑1 (MQ=255)                                                                                               ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg   >  2:1934465/1‑100 (MQ=255)                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg  >  1:179870/1‑100 (MQ=255)                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg  <  1:1311308/100‑1 (MQ=255)                                                                                                 cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg  >  1:785537/1‑99 (MQ=255)                                                                                                   |                                                                                                 TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGG  >  NC_000913/1597174‑1597368

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