Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,313:1	(C)7→8	100%	pseudogene (675/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,306	1	.	C	81.6%	210.5 / 46.7	76	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (5/9);  new base C (25/36);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400                                                                                                      |                                                                                                 agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTATCAGCCGTaacc‑aa                                                                                                   >  2:968404/1‑100 (MQ=255)  gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacc‑aa‑‑c                                                                                                <  2:1280841/100‑1 (MQ=255)  gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacc‑aa‑‑c                                                                                                >  2:816401/1‑100 (MQ=255)  gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacc‑aa‑‑c                                                                                                >  2:2282072/1‑100 (MQ=255)  gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacc‑aa‑‑c                                                                                                <  2:1950162/100‑1 (MQ=255)    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑cc                                                                                             >  2:930905/1‑100 (MQ=255)       tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccc                                                                                          >  1:1398969/1‑100 (MQ=255)       tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccc                                                                                          <  1:1228082/100‑1 (MQ=255)       tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccc                                                                                          >  2:2599306/1‑100 (MQ=255)        gCGTACCGGGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑cccccc                                                                                         <  2:2452239/100‑1 (MQ=255)         cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccccc                                                                                        <  1:1952615/100‑2 (MQ=255)          gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccccct                                                                                       >  1:1653178/1‑98 (MQ=255)           tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccccctc                                                                                      >  1:2414037/1‑97 (MQ=255)            aCCGTGATGCGAGAATCTTGATCCGGCGCCTTCAGGGGTCCGGTATGGGTGGCATGGCCACCGATAAGCGTATTTTCACCCGTAACC‑AA‑‑c‑ccccccctc                                                                                      >  1:1850027/1‑96 (MQ=255)            aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccccctcg                                                                                     <  2:1456554/100‑5 (MQ=255)            aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑c‑ccccccctc                                                                                      >  1:2566040/1‑96 (MQ=255)              cGTGATGCCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGTCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTcgcg                                                                                    <  1:966136/100‑1 (MQ=255)                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATa                                                                                >  2:1609034/1‑99 (MQ=255)                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAg                                                                               >  1:2372689/1‑100 (MQ=255)                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAg                                                                               >  1:1585255/1‑100 (MQ=255)                    gaCAGAGACTTGACCCTGCACATTCAGGGTTCCGATACCGCAGGCACCTACACCCATAATCGCATTTCCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAg                                                                               <  1:280475/98‑1 (MQ=255)                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            <  2:638893/100‑1 (MQ=255)                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            <  1:2885441/100‑1 (MQ=255)                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            >  2:972566/1‑100 (MQ=255)                       agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑ATCCC‑CCCCCCTCGCGAATAGTCg                                                                            <  1:259015/100‑1 (MQ=255)                                 tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGtt                                                                 <  2:2911233/100‑1 (MQ=255)                                   cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                               <  2:2403447/100‑1 (MQ=255)                                    tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              <  2:12877/100‑1 (MQ=255)                                     gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                             >  2:2901397/1‑100 (MQ=255)                                      cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCTAC‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:1793041/100‑1 (MQ=255)                                      cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:263612/100‑1 (MQ=255)                                        caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                          >  2:1772409/1‑100 (MQ=255)                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  2:2346882/100‑1 (MQ=255)                                         aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  2:1561910/100‑1 (MQ=255)                                          ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AACCC‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  1:2805652/100‑1 (MQ=255)                                            cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                      >  2:703275/1‑100 (MQ=255)                                              gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACG‑AAGCC‑CCCCCCTCGCGAGTAGGCGCCTCCCCAGATCCTTGAtttc                                                       >  2:217480/1‑96 (MQ=255)                                              gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                   >  1:2513730/1‑100 (MQ=255)                                               ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCAAATAATCGCCTCCCCAGTTCCTTGATTACCa                                                     >  1:1436702/1‑98 (MQ=255)                                                gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                  <  2:716405/100‑1 (MQ=255)                                                 ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 <  1:2282057/100‑1 (MQ=255)                                                 ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 <  1:816398/100‑1 (MQ=255)                                                   ccGCTACCGGATGCTTTGCCACCGACAAGCGCATTTTCAGCCGGAACC‑AA‑‑C‑CCGCAATCGCAAAGAGTCTCCTCCCCAGTTCCTTGATCACCAAGATGaa                                              <  1:2626177/100‑1 (MQ=255)                                                    cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              >  2:1430057/1‑100 (MQ=255)                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  1:1198758/100‑1 (MQ=255)                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  1:1772401/100‑1 (MQ=255)                                                     attACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  2:2736920/3‑100 (MQ=255)                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          >  1:626257/4‑100 (MQ=255)                                                       tACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                           <  2:460693/100‑1 (MQ=255)                                                        tccgatctcATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          <  1:2651478/92‑1 (MQ=255)                                                        atcGGTGGCAATGCCACCGATAATCGAATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTGGCGAATAGTCGCCTCCCCAGATCCTTGATTACCAATTGGCAATTc                                          >  1:354859/3‑99 (MQ=255)                                                         tcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  1:847554/2‑100 (MQ=255)                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAAATCa                                         >  2:161364/1‑100 (MQ=255)                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                      <  2:2414052/100‑1 (MQ=255)                                                              ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                   >  2:1120910/1‑100 (MQ=255)                                                               gCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                   >  2:2721022/1‑100 (MQ=255)                                                                caatGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:734215/97‑1 (MQ=255)                                                                cATTGCCACCGATTATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTTCCAATTTGAAATTCAATTGAg                                   <  2:217044/99‑1 (MQ=255)                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCagatcgga                                  >  2:2651494/1‑93 (MQ=255)                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:920408/100‑1 (MQ=255)                                                                   tGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                              <  1:402535/100‑1 (MQ=255)                                                                     ccACCGATAATCGTATTTTCAGCCGTAACC‑AACCC‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                              <  1:703466/100‑1 (MQ=255)                                                                      cACCGATAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                            <  1:407968/100‑1 (MQ=255)                                                                      cACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                           >  1:2441468/1‑100 (MQ=255)                                                                         cGATAATCGTATTTTCAGCCGTAACC‑AACCC‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          >  2:2138547/1‑100 (MQ=255)                                                                            tAATCGTATTTTCAGCCGTAACC‑AACCC‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttt                       <  2:2817037/100‑1 (MQ=255)                                                                            tAATCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                      <  1:317506/100‑1 (MQ=255)                                                                             aaTCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  2:599073/100‑1 (MQ=255)                                                                             aaTCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  2:91359/100‑1 (MQ=255)                                                                              aTCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  2:832186/1‑100 (MQ=255)                                                                              aTCGTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  2:1116490/1‑100 (MQ=255)                                                                               tCGTATTTTCAGTCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:1931999/100‑1 (MQ=255)                                                                                 gTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAAATCAATTGAGGAATCATTATTTTttatt                 >  2:732459/1‑100 (MQ=255)                                                                                 gTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 <  2:579704/100‑1 (MQ=255)                                                                                 gTATTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 <  2:2720166/100‑1 (MQ=255)                                                                                   aTTTTCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               >  1:2103376/1‑100 (MQ=255)                                                                                   aTTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:2329438/100‑1 (MQ=255)                                                                                    ttttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:1398974/100‑1 (MQ=255)                                                                                    ttttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:1793075/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             <  2:1143247/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             <  1:201093/100‑1 (MQ=255)                                                                                     tttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             <  1:860494/100‑1 (MQ=255)                                                                                      ttCAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa            >  2:2196313/1‑100 (MQ=255)                                                                                       tCAGCCGTAACC‑AA‑‑C‑CTCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt          <  1:2783186/100‑1 (MQ=255)                                                                                       ncAGCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCCTTATTTTTTATTAGCCAc           >  2:2298834/2‑100 (MQ=255)                                                                                          gCCGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGACTTCCAATTTGAACTTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  1:2328873/1‑99 (MQ=255)                                                                                          gCCGTAACA‑AC‑‑CGCCCCCGTAGCGAATAGTGGCCTCCCCAGTTCCTTGTTTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        <  1:1995723/100‑1 (MQ=255)                                                                                           ccGTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTAGACTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  1:1956235/1‑98 (MQ=255)                                                                                             gTAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCa     >  2:1348657/1‑100 (MQ=255)                                                                                              tAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    <  2:749152/100‑1 (MQ=255)                                                                                              tAACC‑AAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    <  2:382427/100‑1 (MQ=255)                                                                                              tAACC‑AAC‑C‑CCCCCCTCGCGAAGAGTCGCCTCCCCAGTTCCTTGATTACCAAGTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    >  1:1301176/1‑100 (MQ=255)                                                                                               aacc‑aac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc   >  2:13810/1‑100 (MQ=255)                                                                                               aacc‑aac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTCCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtttttatcttct          >  1:85677/1‑87 (MQ=255)                                                                                                acc‑aac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  >  1:1951977/1‑100 (MQ=255)                                                                                                acc‑aac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  2:1397475/100‑1 (MQ=255)                                                                                                acc‑aac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  1:1116485/100‑1 (MQ=255)                                                                                                      |                                                                                                 AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACC‑AA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400

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