Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,952,569	C→T	100%	A29V (GCG→GTG)	ilvE →	branched‑chain amino acid aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,952,569	0	C	T	95.9%	291.8 / 2.2	97	A29V (GCG→GTG)	ilvE	branched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (41/52);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

GGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCG  >  NC_000913/3952473‑3952662                                                                                                 |                                                                                              ggAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCt                                                                                            >  1:1068063/1‑100 (MQ=255)   aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGtgctgc                                                                                          <  2:2101922/100‑1 (MQ=255)   aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGtgctgc                                                                                          <  1:1983828/100‑1 (MQ=255)    atatAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCa                                                                                         >  2:888293/1‑100 (MQ=255)      ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGGACGTGCTGCAc                                                                                        >  1:2378833/1‑99 (MQ=255)      ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                       >  2:226692/1‑100 (MQ=255)      ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                       >  2:2186911/1‑100 (MQ=255)      ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACt                                                                                       <  1:959807/100‑1 (MQ=255)       taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGCAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCAc                                                                                        >  1:926332/1‑98 (MQ=255)       taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTa                                                                                      >  2:538878/1‑100 (MQ=255)          aaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATgg                                                                                   <  2:711563/100‑1 (MQ=255)           gatctccACGAAGAAAGCTGATTACATTTGGTTCTATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCTCTATGGc                                                                                  <  1:1358139/95‑1 (MQ=255)              gACCACGAAGAAAGCTGCGTACATTGGGTTCAATGGGGAGATGGCTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTCTGGc                                                                                  <  1:1454338/97‑1 (MQ=255)                ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTc                                                                             >  1:824459/1‑100 (MQ=255)                ccACGAACAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCagatc                                                                             >  2:1358147/1‑96 (MQ=255)                 cACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCg                                                                            <  2:599954/100‑1 (MQ=255)                   cGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGt                                                                          >  2:590075/1‑100 (MQ=255)                   cGAAGAAAGCTGATTACATTAGGTTCAATGGGGAGATGGTTCGCTAGGAAGACGCGATGGTGCATGTGATGTCGCATGTGCTGCACTATGGCTCTTCGGt                                                                          >  1:2441681/1‑100 (MQ=255)                       gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttttt                                                                      <  1:1316181/100‑1 (MQ=255)                       gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttttt                                                                      >  1:2226053/1‑100 (MQ=255)                        aaaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTg                                                                     <  1:29249/100‑1 (MQ=255)                        aaaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTg                                                                     <  2:2019742/100‑1 (MQ=255)                         aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGa                                                                    >  2:118975/1‑100 (MQ=255)                           ncTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                  >  2:1411534/2‑100 (MQ=255)                           gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                  <  2:538978/100‑1 (MQ=255)                           gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                  <  2:690486/100‑1 (MQ=255)                            cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                 <  2:2179720/100‑1 (MQ=255)                            cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                 <  1:609223/100‑1 (MQ=255)                            cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                 <  1:1078086/100‑1 (MQ=255)                            cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                 <  1:2212347/100‑1 (MQ=255)                             tGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGc                                                                <  2:2138474/100‑1 (MQ=255)                               aTTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCa                                                               <  2:2277208/99‑1 (MQ=255)                                ttACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATc                                                             >  2:2279423/1‑100 (MQ=255)                                tgaCATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATc                                                             >  2:1757070/3‑100 (MQ=255)                                     tttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTg                                                        <  1:1494572/100‑1 (MQ=255)                                        ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTa                                                     >  1:706621/1‑100 (MQ=255)                                         gTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTAc                                                    >  2:2151809/1‑100 (MQ=255)                                          ttCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACg                                                   <  2:1032088/100‑1 (MQ=255)                                           tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGa                                                  >  2:1898781/1‑100 (MQ=255)                                           tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGa                                                  >  1:2062301/1‑100 (MQ=255)                                             aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACt                                                >  1:1691405/1‑100 (MQ=255)                                              aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTc                                               >  1:2033294/1‑100 (MQ=255)                                                 gggAGATGGTTCGCTGGGAAGACGCGAAGGTGTATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGc                                             <  1:856316/99‑1 (MQ=255)                                                 gggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGca                                            <  1:2533595/100‑1 (MQ=255)                                                  ggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcac                                           <  2:644303/100‑1 (MQ=255)                                                  ggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcac                                           >  2:211948/1‑100 (MQ=255)                                                  ggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcac                                           >  1:778772/1‑100 (MQ=255)                                                   gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcaca                                          <  1:2515511/100‑1 (MQ=255)                                                    agaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACaa                                         <  1:2268295/100‑1 (MQ=255)                                                      aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAg                                       <  2:1717902/100‑1 (MQ=255)                                                      aTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAg                                       <  1:1325962/100‑1 (MQ=255)                                                       tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                      <  2:2338008/100‑1 (MQ=255)                                                         gTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAc                                    >  1:760831/1‑100 (MQ=255)                                                         gTTCGCTGGGAAGAAGCGAAGGTGAATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAc                                    <  1:492258/100‑1 (MQ=255)                                                          ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                   <  2:2033306/100‑1 (MQ=255)                                                          ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                   <  2:446168/100‑1 (MQ=255)                                                          ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                   >  2:2010520/1‑100 (MQ=255)                                                           tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCg                                  >  1:1447296/1‑100 (MQ=255)                                                            cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                 >  2:2024968/1‑100 (MQ=255)                                                             gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                >  2:772566/1‑100 (MQ=255)                                                             gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                <  2:784719/100‑1 (MQ=255)                                                             gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                >  2:1132242/1‑100 (MQ=255)                                                               tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                              <  1:690502/100‑1 (MQ=255)                                                               tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                              <  2:1332552/100‑1 (MQ=255)                                                                gggAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                             <  1:643208/100‑1 (MQ=255)                                                                gggAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                             <  1:1948049/100‑1 (MQ=255)                                                                 ggAAGACGCGAAGGTGCATGGGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                            <  1:594976/100‑1 (MQ=255)                                                                  gAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                           >  2:2392143/1‑100 (MQ=255)                                                                   aaGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAtt                          >  1:40728/1‑100 (MQ=255)                                                                      aCGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                       <  1:1898771/100‑1 (MQ=255)                                                                       ctcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                      >  2:1759570/3‑100 (MQ=255)                                                                       cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                      <  1:984242/100‑1 (MQ=255)                                                                       cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                      <  2:2263960/100‑1 (MQ=255)                                                                       cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                      <  2:104018/100‑1 (MQ=255)                                                                       cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGATTGTATTCCGc                      >  1:1665619/1‑100 (MQ=255)                                                                         cgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                    >  1:818492/1‑100 (MQ=255)                                                                         cgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATCCCGCCa                    >  1:2059532/1‑100 (MQ=255)                                                                         cgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGGATTCCGCCa                    >  1:1572346/1‑100 (MQ=255)                                                                           aaGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATc                  <  2:1691230/100‑1 (MQ=255)                                                                            aGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAATGGACCGGTTGTATTCCGCCATCg                 <  1:984636/100‑1 (MQ=255)                                                                             ggTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGt                <  2:2101043/100‑1 (MQ=255)                                                                             ggTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGt                <  2:503934/100‑1 (MQ=255)                                                                                  aTGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCa           <  1:863605/100‑1 (MQ=255)                                                                                  aTGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCa           <  2:511711/100‑1 (MQ=255)                                                                                  aTGTGATGGCGCATGTGCTGCACTATGGAACTTCGGTTTGTAAAGTTATCCGTTGCTAACGCTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc            <  1:805125/99‑1 (MQ=255)                                                                                   tgtgATGTCGCACGTGCTGCACTATGGCACTTCGGTCTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCat          <  2:472316/100‑1 (MQ=255)                                                                                    gtgATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata         <  2:2393874/100‑1 (MQ=255)                                                                                    gtgATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata         <  1:717948/100‑1 (MQ=255)                                                                                     ttaTGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCatat        >  2:1742859/3‑100 (MQ=255)                                                                                       aTGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc      <  2:2082146/100‑1 (MQ=255)                                                                                       aTGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCatattc      >  1:1569560/1‑98 (MQ=255)                                                                                         gTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAg    <  1:728456/100‑1 (MQ=255)                                                                                         gTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAg    >  1:2529219/1‑100 (MQ=255)                                                                                          tCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGc      >  1:994369/1‑97 (MQ=255)                                                                                           cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg  >  2:2360319/1‑100 (MQ=255)                                                                                           cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg  <  1:2147016/100‑1 (MQ=255)                                                                                           cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg  >  2:792197/1‑100 (MQ=255)                                                                                           cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg  <  1:2406027/100‑1 (MQ=255)                                                                                             cACGTGCTGCACTATGGCACATCGGTTTTTGAAGGTATCCGTTGCTACGAGTCGCCCAAAGGACCGGTTGTATTCCGACATCGTGAGCATATGCAGCg  <  1:1267033/98‑1 (MQ=255)                                                                                                 |                                                                                              GGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCG  >  NC_000913/3952473‑3952662

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: ‑