Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)7→8	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	91.9%	244.7 / 16.4	74	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base C (32/36);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597208‑1597400                                                                                                     |                                                                                                 gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑a‑c                                                                                                 <  1:2031301/100‑1 (MQ=255)  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑a‑cc                                                                                                <  2:1587053/100‑1 (MQ=255)    gcTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑ccc                                                                                               <  2:1680590/98‑1 (MQ=255)      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccc‑cc                                                                                           <  2:1786890/100‑1 (MQ=255)       gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccc‑cc                                                                                           <  1:709190/99‑1 (MQ=255)       gCGTACCGTGATGACAGAGTCTTGATACTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccagat                                                                                           >  1:1680603/1‑96 (MQ=255)         gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccc‑cc‑cct                                                                                       <  2:2696363/100‑3 (MQ=255)          tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑cccc‑cc‑cct                                                                                       >  2:2571675/1‑97 (MQ=255)               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGaa                                                                                  <  2:2629596/100‑1 (MQ=255)               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGaa                                                                                  <  2:2354811/100‑1 (MQ=255)               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGaa                                                                                  >  1:1729239/1‑100 (MQ=255)                gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAAt                                                                                 <  2:1983046/100‑1 (MQ=255)                gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAAt                                                                                 >  2:1471864/1‑100 (MQ=255)                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAg                                                                               <  1:1085903/100‑1 (MQ=255)                   gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGt                                                                              >  2:2328576/1‑100 (MQ=255)                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCg                                                                            >  1:1855043/1‑100 (MQ=255)                       gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGcc                                                                          <  2:1283678/100‑1 (MQ=255)                         gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGc                                                                           >  2:915538/1‑97 (MQ=255)                         gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTc                                                                        <  2:528594/100‑1 (MQ=255)                         gTCTTGATCCTGCACATCCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATCTTCAGCCGTAGCCA‑ACCCCC‑CC‑CTCGCGAATACTc                                                                             >  2:749992/1‑95 (MQ=255)                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTcccc                                                                     <  2:2044619/100‑1 (MQ=255)                               aTCCTCCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCAGTAACCT‑A‑CCCC‑CCACTCGCGAATAGTCGGCTccaccat                                                                  >  2:1360153/1‑95 (MQ=255)                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTc                                                                <  1:1415997/100‑1 (MQ=255)                                 catGCACATTCAGAGTTCCGATACCGGTGGCATTGCCACCGAGAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTc                                                                >  1:1862097/3‑100 (MQ=255)                                  cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAGTCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTcc                                                               <  2:1766607/100‑1 (MQ=255)                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTcc                                                               >  2:145152/1‑99 (MQ=255)                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              >  2:383041/1‑100 (MQ=255)                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                             <  2:385286/100‑1 (MQ=255)                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              <  1:2519931/99‑1 (MQ=255)                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTAGCCTCCCCAGTTCCtt                                                             >  2:1617325/1‑100 (MQ=255)                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                             >  1:24873/1‑99 (MQ=255)                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                             <  1:2729613/99‑1 (MQ=255)                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  2:1308961/100‑1 (MQ=255)                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            >  1:458046/1‑100 (MQ=255)                                       caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATTGTCGCCTCCCCAGTTCCTTGAt                                                          <  2:1601981/100‑1 (MQ=255)                                       caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                          >  1:1486864/1‑100 (MQ=255)                                       caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                          <  1:2085836/100‑1 (MQ=255)                                       caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                          <  2:1264284/100‑1 (MQ=255)                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  2:2508534/1‑100 (MQ=255)                                          tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                       >  2:708184/1‑100 (MQ=255)                                            aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                                                     >  2:955390/1‑100 (MQ=255)                                              ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                   >  2:644671/1‑100 (MQ=255)                                               gTTCCGATACCGGTGGCCTTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                  <  2:666714/100‑1 (MQ=255)                                                 tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:2048710/100‑1 (MQ=255)                                                  ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:1365734/99‑1 (MQ=255)                                                       aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          <  1:874290/100‑1 (MQ=255)                                                       aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          <  2:1161368/100‑1 (MQ=255)                                                        ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  1:623625/1‑100 (MQ=255)                                                         cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATATGCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                        >  2:1508486/1‑100 (MQ=255)                                                                  tGCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAt                               <  2:1729225/100‑1 (MQ=255)                                                                   gCCACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                              >  1:1538471/1‑100 (MQ=255)                                                                    ccACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                            >  2:2716331/1‑100 (MQ=255)                                                                     cACCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                            <  2:2566495/100‑1 (MQ=255)                                                                      tccGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                           <  2:102405/99‑1 (MQ=255)                                                                      aCCGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                           <  1:603751/100‑1 (MQ=255)                                                                       ccGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATGGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          >  1:1592917/1‑100 (MQ=255)                                                                       ccGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCTCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          <  1:1592310/100‑1 (MQ=255)                                                                       ccGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          <  2:2105561/100‑1 (MQ=255)                                                                       ccGATAATCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          >  1:102407/1‑100 (MQ=255)                                                                            aaTCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     <  1:1272200/100‑1 (MQ=255)                                                                            aaTCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     <  2:789825/100‑1 (MQ=255)                                                                             aTCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  2:100887/100‑1 (MQ=255)                                                                              tCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   >  1:821971/1‑100 (MQ=255)                                                                              tCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:145154/100‑1 (MQ=255)                                                                              tCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   >  2:811909/1‑100 (MQ=255)                                                                              tCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:2392806/100‑1 (MQ=255)                                                                              tCGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  2:33018/100‑1 (MQ=255)                                                                               cGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  >  2:499570/1‑100 (MQ=255)                                                                               cGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  >  2:1167233/1‑100 (MQ=255)                                                                               cGTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  <  2:1666702/100‑1 (MQ=255)                                                                                gcaTTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  1:1978801/3‑100 (MQ=255)                                                                                gTATTTTCAGCCGTAACCACC‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 <  1:1123754/100‑1 (MQ=255)                                                                                gTATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  2:583539/1‑100 (MQ=255)                                                                                gTATTTTCAGCCGTAACCA‑A‑CCCCTCC‑CTCGCGAATAGTCTCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  <  2:2756294/99‑1 (MQ=255)                                                                                 tATTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatta                <  2:2347744/100‑1 (MQ=255)                                                                                  aTTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               <  2:2045188/100‑1 (MQ=255)                                                                                  aTTTTCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               >  1:1297848/1‑100 (MQ=255)                                                                                    tttCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGCGGAATCATTATTTTTTATTAGcc             >  1:550206/1‑100 (MQ=255)                                                                                    tttCAGCCGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              >  1:2585583/1‑99 (MQ=255)                                                                                        aGCCGTAACCA‑A‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAATTTAAATTCAATTGAAGAATCATTATTTTTTATTAGCCACTCg        >  1:1508916/1‑100 (MQ=255)                                                                                          ccGTAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc       >  2:1908082/1‑100 (MQ=255)                                                                                             tAACCA‑ACCCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    >  2:1257001/1‑100 (MQ=255)                                                                                               acca‑accCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTCGCCACTCGCCAccc  >  2:1878782/1‑100 (MQ=255)                                                                                               acca‑accCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  2:113967/100‑1 (MQ=255)                                                                                               acca‑accCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  1:955396/100‑1 (MQ=255)                                                                                               acca‑accCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCCACAccc  >  1:552311/1‑100 (MQ=255)                                                                                                     |                                                                                                 GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑CCCC‑CC‑CTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597208‑1597400

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