Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,842,718	T→C	T60A (ACA→GCA)	b1762 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,718	0	T	C	99.0%	293.5 / ‑4.4	95	T60A (ACA→GCA)	b1762	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (43/51);  total (43/52)
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.88e-01

CGGCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGTTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTAATC  >  NC_000913/1842628‑1842806                                                                                           |                                                                                         cGGCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTa                                                                                        <  2:4297309/93‑1 (MQ=255)  ggCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGTTAt                                                                                       <  1:4322618/93‑1 (MQ=255)  ggCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTAt                                                                                       <  2:4480197/93‑1 (MQ=255)   gcacGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATc                                                                                      <  1:2425752/90‑1 (MQ=255)   gCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATc                                                                                      <  2:2741559/93‑1 (MQ=255)   gCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATc                                                                                      <  1:2326931/93‑1 (MQ=255)   gCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATc                                                                                      <  2:4275397/93‑1 (MQ=255)    cccGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCg                                                                                     >  1:4014481/1‑93 (MQ=255)     ccGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGc                                                                                    <  2:2576511/93‑1 (MQ=255)     ccGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGc                                                                                    <  1:3411433/93‑1 (MQ=255)        aaCAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGc                                                                                 <  1:785832/93‑1 (MQ=255)         aCAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCg                                                                                >  2:2003693/1‑93 (MQ=255)            aTAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAac                                                                             >  1:4172795/1‑93 (MQ=255)             tAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCCTCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAaca                                                                            <  1:3141984/93‑1 (MQ=255)             tAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAaca                                                                            <  1:930944/93‑1 (MQ=255)             tAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAaca                                                                            <  1:1073496/93‑1 (MQ=255)               aGCCCCATGATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAAcaca                                                                          <  1:1303271/93‑1 (MQ=255)                 ccccATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACaaa                                                                        <  2:514884/93‑1 (MQ=255)                  cccATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAAc                                                                       <  2:3218227/93‑1 (MQ=255)                  cccATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAAc                                                                       >  2:3478113/1‑93 (MQ=255)                   ccATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACa                                                                      <  1:727135/93‑1 (MQ=255)                     atcatAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACAtt                                                                    <  1:3389865/93‑1 (MQ=255)                      tcatAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTa                                                                   >  1:546960/1‑93 (MQ=255)                      tcatAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTa                                                                   >  2:3811375/1‑93 (MQ=255)                      tcatAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTa                                                                   >  2:2702004/1‑93 (MQ=255)                        atAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACa                                                                 >  1:2671675/1‑93 (MQ=255)                         taattaatAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAg                                                                >  1:1446020/1‑93 (MQ=255)                          aattaatAGAAAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAACCCGTTGGATGCTATCGCTGCGAACACAAACATTACAGc                                                               <  1:1300146/93‑1 (MQ=255)                             taatAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGac                                                            >  1:1795426/1‑93 (MQ=255)                             taatAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGac                                                            >  2:536498/1‑93 (MQ=255)                               atAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGacaa                                                          <  2:486043/93‑1 (MQ=255)                                tAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGacaac                                                         >  1:1660631/1‑93 (MQ=255)                                 aGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGacaaca                                                        <  2:1804601/93‑1 (MQ=255)                                 aGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGacaaca                                                        >  1:1435629/1‑93 (MQ=255)                                   aTAAAGAAAAGAGCATCCCGCGGTACTCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAAcaca                                                      <  1:4528470/93‑1 (MQ=255)                                   aTAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAAcaca                                                      >  1:236813/1‑93 (MQ=255)                                     aaagaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATg                                                    >  2:1404795/1‑93 (MQ=255)                                     aaagaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATg                                                    >  2:614764/1‑93 (MQ=255)                                     aaagaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATg                                                    <  1:789789/93‑1 (MQ=255)                                      aagaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGt                                                   >  1:3825925/1‑93 (MQ=255)                                      aagaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGt                                                   >  1:2167321/1‑93 (MQ=255)                                       agaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGtt                                                  >  2:4262531/1‑93 (MQ=255)                                       agaaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGtt                                                  <  2:1124505/93‑1 (MQ=255)                                         aaaagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTtaa                                                >  2:2954102/1‑93 (MQ=255)                                           aagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTtaaat                                              >  2:850965/1‑93 (MQ=255)                                           aagaGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTtaaat                                              >  2:4278875/1‑93 (MQ=255)                                            agagCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTtaaata                                             >  1:2867394/1‑93 (MQ=255)                                            agagCACCCCGCGGTAATCAGGCCTGTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTtaaata                                             <  1:854472/93‑1 (MQ=255)                                               gCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATaaaa                                          >  1:2457402/1‑93 (MQ=255)                                               gCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATaaaa                                          <  2:2522869/93‑1 (MQ=255)                                                 aTCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAAcc                                        <  2:1605331/93‑1 (MQ=255)                                                  tCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAAgcc                                       >  1:2846148/1‑93 (MQ=255)                                                     cgcgGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCa                                    <  2:3222710/93‑1 (MQ=255)                                                       cgGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACt                                  <  1:1014130/93‑1 (MQ=255)                                                        ggTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACtt                                 <  2:1606388/93‑1 (MQ=255)                                                         gTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTa                                >  2:1194970/1‑93 (MQ=255)                                                         gTAATCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTa                                <  1:4621306/93‑1 (MQ=255)                                                           aaTCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACa                              >  2:3026341/1‑93 (MQ=255)                                                            aTCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACaa                             >  2:2773141/1‑93 (MQ=255)                                                             tCAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACaaa                            <  1:908128/93‑1 (MQ=255)                                                              cAGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAag                           <  1:1006141/93‑1 (MQ=255)                                                               aGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagc                          >  1:2228838/1‑93 (MQ=255)                                                               aGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagc                          >  2:902420/1‑93 (MQ=255)                                                               aGGCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAaga                          >  1:3477644/1‑92 (MQ=255)                                                                ggCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagca                         <  2:2420692/93‑1 (MQ=255)                                                                ggCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagca                         <  1:2938687/93‑1 (MQ=255)                                                                ggCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagca                         <  1:2163141/93‑1 (MQ=255)                                                                ggCCTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagca                         >  2:4237198/1‑93 (MQ=255)                                                                  ccTCTTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAagcagc                       <  2:4305481/93‑1 (MQ=255)                                                                   ctctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACCACACATGTTAAATAAAACCCCCACTTACAAAgcagca                      >  2:2019703/1‑93 (MQ=255)                                                                   ctctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAgcagca                      >  2:4094676/1‑93 (MQ=255)                                                                   ctctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAgcagca                      >  2:3973654/1‑93 (MQ=255)                                                                   ctctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAAACCCCACTTACAAAgcagca                      >  2:503159/1‑93 (MQ=255)                                                                    tctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAATCCCCCACTTACAAAGCAGCAt                     <  2:364286/93‑1 (MQ=255)                                                                     ctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATg                    <  2:4002509/93‑1 (MQ=255)                                                                     ctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATg                    <  2:927006/93‑1 (MQ=255)                                                                     ctTTCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATg                    <  2:2421567/93‑1 (MQ=255)                                                                      tttCATTCAAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATgg                   >  1:1605886/1‑93 (MQ=255)                                                                         cattcaAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCt                <  2:3531333/93‑1 (MQ=255)                                                                         cattcaAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCt                <  2:4411379/93‑1 (MQ=255)                                                                         cattcaAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCt                <  2:2271430/93‑1 (MQ=255)                                                                           ttcaAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTa              >  2:3861830/1‑93 (MQ=255)                                                                           ttcaAACCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTa              <  2:155405/93‑1 (MQ=255)                                                                              aaaCCGTTGGATGCTATCGCTGCGATCTCAATCATTACAGCGACAACACTTGTTAAATAAAACCCCCTCTTACAAAGCAGCATGGCCTTAttt           <  2:1026531/93‑1 (MQ=255)                                                                              aaaCCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTAttt           <  1:2228866/93‑1 (MQ=255)                                                                               aaCCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTAtttt          <  1:2063908/93‑1 (MQ=255)                                                                               aaCCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTAtttt          >  2:753140/1‑93 (MQ=255)                                                                               aaCCGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTAtttt          >  1:1764036/1‑93 (MQ=255)                                                                                  cGTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCtt       <  1:967843/93‑1 (MQ=255)                                                                                   gTTGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCttt      >  1:1085769/1‑93 (MQ=255)                                                                                    ttGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTa     <  2:2299277/93‑1 (MQ=255)                                                                                    ttGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTa     <  1:4520887/93‑1 (MQ=255)                                                                                     tGGATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTaa    >  2:3306926/1‑93 (MQ=255)                                                                                       gATGCTATCGCTGCGAACACAACCATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTAatc  <  2:3398299/93‑1 (MQ=255)                                                                                       gATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTAatc  >  2:1366691/1‑93 (MQ=255)                                                                                       gATGCTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTAatc  <  1:2037086/93‑1 (MQ=255)                                                                                           |                                                                                         CGGCCCGAACAATAAGCCCCATCATAATTAATAGATAAAGAAAAGAGCATCCCGCGGTAATCAGGCCTCTTTCATTCAAACCGTTGGATGTTATCGCTGCGAACACAAACATTACAGCGACAACACATGTTAAATAAAACCCCCACTTACAAAGCAGCATGGCCTTATTTTCTTTAATC  >  NC_000913/1842628‑1842806

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