Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,981,237	A→C	6.0%	F316V (TTT→GTT)	araE ←	arabinose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,981,237	0	A	C	6.0%	205.2 / 4.9	67	F316V (TTT→GTT)	araE	arabinose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base C (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGTAAACTGCTGCATCGCC  >  NC_000913/2981089‑2981385                                                                                                                                                     |                                                                                                                                                     gTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAac                                                                                                                                                     >  1:979221/1‑150 (MQ=255)   cATGCCAACAGAGAGCCAGGACAAGAAACTGGAAGCCGTACCGTAATCAAACTGCATCAGCCAAAAGACAAGCACCAGAGTGCCTAACGCCATCAAGCAGAAACCAATATTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGAAAacac                                                                                                                                                   >  2:568458/1‑150 (MQ=255)      gCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGc                                                                                                                                                <  2:979221/150‑1 (MQ=255)      gCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGc                                                                                                                                                <  2:317898/150‑1 (MQ=255)      gCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGc                                                                                                                                                <  2:568685/150‑1 (MQ=255)         aaCAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAAt                                                                                                                                             <  1:750877/150‑1 (MQ=255)               gagCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGt                                                                                                                                       >  2:1944480/1‑150 (MQ=255)                 gCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTgg                                                                                                                                     >  1:1944828/1‑150 (MQ=255)                  ccAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGc                                                                                                                                    <  2:716724/150‑1 (MQ=255)                     ggACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGaa                                                                                                                                 >  2:1951087/1‑150 (MQ=255)                           gCCACTTGAAGACGTACCGTTATCAAAATGCAACAGGAAATATCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTAAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATaaa                                                                                                                           >  2:847626/1‑150 (MQ=255)                           gCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATaaa                                                                                                                           >  1:134509/1‑150 (MQ=255)                           gCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATaaa                                                                                                                           >  1:61469/1‑150 (MQ=255)                            ccACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAg                                                                                                                          >  2:965416/1‑150 (MQ=255)                            ccACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAg                                                                                                                          <  1:1944480/150‑1 (MQ=255)                             cACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAgg                                                                                                                         <  1:1133059/150‑1 (MQ=255)                               cTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTc                                                                                                                       >  2:164143/1‑150 (MQ=255)                                 ggAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAg                                                                                                                     >  2:387308/1‑150 (MQ=255)                                          aCCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGAc                                                                                                            <  2:1787544/150‑1 (MQ=255)                                           ccGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGAcc                                                                                                           >  2:62468/1‑150 (MQ=255)                                                 tCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTc                                                                                                     <  1:1993707/150‑1 (MQ=255)                                                   aaaCTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGc                                                                                                   >  2:747589/1‑150 (MQ=255)                                                   aaaCTGCATCAGGCAATAGACCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGc                                                                                                   >  2:236362/1‑150 (MQ=255)                                                     aCTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCaa                                                                                                 >  1:1954580/1‑150 (MQ=255)                                                        gCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGGCAGCCCTACGACCAGAGTCGCAatca                                                                                              <  2:699578/150‑1 (MQ=255)                                                          aTCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAatcatc                                                                                            >  1:434822/1‑150 (MQ=255)                                                              ggCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCtgtt                                                                                        >  1:613083/1‑150 (MQ=255)                                                               gCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCtgttg                                                                                       <  2:61469/150‑1 (MQ=255)                                                                    aGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAAc                                                                                                                                >  2:113947/1‑104 (MQ=255)                                                                    aGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAAc                                                                                                                                <  1:113947/104‑1 (MQ=255)                                                                       ccAGAACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCTGTTTACGCCCTGCTTTATCTACAGTAAACACCGAAATAAAGGTGGCGAACATAAAGGTCAGCCATACGACCAGAGTCGCAATCATGTGTTGTTATGTgg                                                                               >  2:553056/1‑150 (MQ=255)                                                                             ccAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTaa                                                                         >  1:602830/1‑150 (MQ=255)                                                                               agagTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAg                                                                       <  2:602830/150‑1 (MQ=255)                                                                                gagTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGc                                                                      >  2:713056/1‑150 (MQ=255)                                                                                    gCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAACCACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGc                                                                  >  1:457414/1‑150 (MQ=255)                                                                                    gCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAACCACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGc                                                                  >  1:464504/1‑150 (MQ=255)                                                                                          cGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCAtttt                                                            <  1:847626/150‑1 (MQ=255)                                                                                          cGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGGCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCAtttt                                                            <  2:1519265/150‑1 (MQ=255)                                                                                          cGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTAATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCAtttt                                                            >  2:1286075/1‑150 (MQ=255)                                                                                           gCCATCCCGCGGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTg                                                           <  2:1254505/150‑1 (MQ=255)                                                                                            ccATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGa                                                          <  1:1670474/150‑1 (MQ=255)                                                                                             cATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGaa                                                         <  1:711886/150‑1 (MQ=255)                                                                                               tCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGa                                                       <  1:62468/150‑1 (MQ=255)                                                                                                   gCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATAcg                                                   <  2:434822/150‑1 (MQ=255)                                                                                                         aCCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGgcgc                                             >  1:330228/1‑150 (MQ=255)                                                                                                            aaTTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTc                                                                                                                       >  2:945060/1‑73 (MQ=255)                                                                                                            aaTTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTc                                                                                                                       <  1:945060/73‑1 (MQ=255)                                                                                                            aaTTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCgta                                          <  1:470595/150‑1 (MQ=255)                                                                                                                ttCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTAc                                      >  1:1623095/1‑150 (MQ=255)                                                                                                                  cAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACat                                    <  1:164143/150‑1 (MQ=255)                                                                                                                  cAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACat                                    >  2:435207/1‑150 (MQ=255)                                                                                                                   agagCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGCCGCgtatcttcaa                                   <  2:892093/150‑8 (MQ=255)                                                                                                                      gCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACatgatg                                >  1:703805/1‑150 (MQ=255)                                                                                                                       ccGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGCCCAGAGTCGCAATCAGCTGTTGTTCTGTGGTCGTAAAGCCCGCCAGTTTGAAGATACGCGGCGCGTAGTACatgagga                               >  2:773397/1‑150 (MQ=255)                                                                                                                           tttACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGtt                           <  2:613083/150‑1 (MQ=255)                                                                                                                            ttACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTc                          <  2:1623095/150‑1 (MQ=255)                                                                                                                             tACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCa                         >  1:1681917/1‑150 (MQ=255)                                                                                                                              aCGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCAt                        <  1:568458/150‑1 (MQ=255)                                                                                                                               cGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATa                       >  2:892073/1‑150 (MQ=255)                                                                                                                               cGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATa                       >  2:1201382/1‑150 (MQ=255)                                                                                                                                 cccTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCgt                                           <  2:575764/128‑1 (MQ=255)                                                                                                                                 cccTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCgt                                           >  1:575764/1‑128 (MQ=255)                                                                                                                                   cTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCgg                   <  1:1471844/150‑1 (MQ=255)                                                                                                                                    tGCTTTATCTACCGTAACCACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGt                  <  2:464504/150‑1 (MQ=255)                                                                                                                                    tGCTTTATCTACCGTAACCACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGt                  <  2:457414/150‑1 (MQ=255)                                                                                                                                            cTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGTAAActgct          >  1:521620/1‑150 (MQ=255)                                                                                                                                              aCCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATAcgc                                                  <  1:448581/108‑1 (MQ=255)                                                                                                                                              aCCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATAcgc                                                  >  2:448581/1‑108 (MQ=255)                                                                                                                                               ccGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGTAAACTGCTGCa       >  1:1428145/1‑150 (MQ=255)                                                                                                                                                    aaCACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGAAACGCGGCGCGTAGTACATGATGATGTTCATACCGGTAAACTGCTGCATCGcc  <  1:387308/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTACCGTTATCAAACTGCATCAGGCAATAGCCCAGCACCAGAGTGCCTAACGCCATCACGCTGAAACCAATTTTCAGAGCCGGTTTACGCCCTGCTTTATCTACCGTAAACACCGCAATAAAGGTGGCGAACATAAAGGTCAGCCCTACGACCAGAGTCGCAATCATCTGTTGTTCTGTGGTCGTAAAGCCCGCCATTTTGAAGATACGCGGCGCGTAGTACATGATGATGTTCATACCGGTAAACTGCTGCATCGCC  >  NC_000913/2981089‑2981385

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