Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,140,401	T→C	8.3%	P404P (CCA→CCG)	asmA ←	putative assembly protein AsmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,140,401	0	T	C	8.3%	166.6 / 6.0	60	P404P (CCA→CCG)	asmA	putative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base C (3/2);  total (31/29)
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.52e-01

GTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTCAT  >  NC_000913/2140256‑2140550                                                                                                                                                  |                                                                                                                                                      gtgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCt                                                                                                                                                   <  2:490923/150‑1 (MQ=255) gtgCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGACt                                                                                                                                                   <  2:1297909/150‑1 (MQ=255)    ccTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGa                                                                                                                                                <  2:754908/150‑1 (MQ=255)       gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGt                                                                                                                                             >  1:672060/1‑150 (MQ=255)         ccTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacg                                                                                                                                           <  2:549338/150‑5 (MQ=255)          cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                          >  2:1255059/1‑150 (MQ=255)           ttGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacgtt                                                                                                                                         <  1:163768/150‑7 (MQ=255)            tGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATc                                                                                                                                        >  2:1209956/1‑150 (MQ=255)              ccAGTTGTGGCGGAATGCGTCGGCATCTATGTCAACACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcg                                                                                                                                      >  2:186540/1‑150 (MQ=255)                aGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcgcg                                                                                                                                    <  1:723084/150‑1 (MQ=255)                      ggCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTa                                                                                                                              >  1:1219046/1‑150 (MQ=255)                        cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGTCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                            >  1:1385900/1‑150 (MQ=255)                        cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAtt                                                                                                                            >  2:1075064/1‑150 (MQ=255)                           aaTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGAt                                                                                                                         <  2:1124675/150‑1 (MQ=255)                              gCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTc                                                                                                                      >  2:249912/1‑150 (MQ=255)                                gTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcgc                                                                                                                    >  1:872059/1‑150 (MQ=255)                                            gTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGcc                                                                                                        <  2:1217331/150‑1 (MQ=255)                                              cAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCg                                                                                                      <  2:872059/150‑1 (MQ=255)                                                    cGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTg                                                                                                <  2:1219046/150‑1 (MQ=255)                                                         aaaTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAcc                                                                                           <  2:311616/150‑1 (MQ=255)                                                           aTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTg                                                                                         <  1:1209956/150‑1 (MQ=255)                                                              aCCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTcc                                                                                      <  1:295718/150‑1 (MQ=255)                                                                 agccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAcc                                                                                   <  1:1255059/150‑1 (MQ=255)                                                                  gccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCg                                                                                  >  2:468056/1‑150 (MQ=255)                                                                      gTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                              >  2:1097342/1‑150 (MQ=255)                                                                        gACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGTATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGt                                                                            >  1:1418048/1‑150 (MQ=255)                                                                          cATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCa                                                                                     >  2:1182275/1‑139 (MQ=255)                                                                          cATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCa                                                                                     <  1:1182275/139‑1 (MQ=255)                                                                             ttttCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGcc                                                                       <  1:418221/150‑1 (MQ=255)                                                                              tttCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCATCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCTTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGcgc                                                                      <  2:720686/150‑3 (MQ=255)                                                                                  cGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCa                                                                  >  2:1065701/1‑150 (MQ=255)                                                                                   ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                                 >  2:1009288/1‑150 (MQ=255)                                                                                   ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATAGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                                 >  1:314444/1‑150 (MQ=255)                                                                                       aaCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCTGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                             >  2:53814/1‑150 (MQ=255)                                                                                        aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAg                                                            >  2:1950/1‑150 (MQ=255)                                                                                           aaaTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTaa                                                         <  1:1075064/150‑1 (MQ=255)                                                                                             aTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAtt                                                       >  1:516959/1‑150 (MQ=255)                                                                                                ggATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTcc                                                    >  1:44534/1‑150 (MQ=255)                                                                                                  aTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAg                                                  <  1:1950/150‑1 (MQ=255)                                                                                                   tAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGc                                                 >  1:119342/1‑150 (MQ=255)                                                                                                       tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTa                                                                               <  1:698140/116‑1 (MQ=255)                                                                                                       tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTa                                                                               >  2:698140/1‑116 (MQ=255)                                                                                                             ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACt                                       <  1:1097342/150‑1 (MQ=255)                                                                                                              cTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACtt                                      >  2:838404/1‑150 (MQ=255)                                                                                                                 cAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                   <  1:1065701/150‑1 (MQ=255)                                                                                                                    aaTGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGt                                <  1:468056/150‑1 (MQ=255)                                                                                                                        gTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                            >  2:1035309/1‑150 (MQ=255)                                                                                                                          aCCAATCTCAACGTTTTCCAGCCGTGGATGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTg                          <  2:44534/150‑1 (MQ=255)                                                                                                                               tctcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTgg                     <  2:672060/150‑1 (MQ=255)                                                                                                                                 tcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                   >  1:371242/1‑150 (MQ=255)                                                                                                                                  cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCaa                  >  1:182473/1‑150 (MQ=255)                                                                                                                                  cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                   >  2:647370/1‑149 (MQ=255)                                                                                                                                  cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                   <  1:647370/149‑1 (MQ=255)                                                                                                                                      gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc              <  2:516959/150‑1 (MQ=255)                                                                                                                                      gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc              <  1:1109778/150‑1 (MQ=255)                                                                                                                                      aacttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc              >  1:405326/4‑150 (MQ=255)                                                                                                                                       acttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCt             >  1:1238959/3‑150 (MQ=255)                                                                                                                                        cttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTg            >  1:818527/2‑150 (MQ=255)                                                                                                                                             aGCCGTGGCTGGAAGTTTATCCCCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaaaa       >  2:725867/1‑150 (MQ=255)                                                                                                                                                cGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAattc    <  2:119342/150‑1 (MQ=255)                                                                                                                                                  tGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACGTCTGTAAAattcat  >  2:663831/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTCAT  >  NC_000913/2140256‑2140550

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