Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,957,104	A→T	7.9%	F594Y (TTT→TAT)	ptrA ←	protease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,957,104	0	A	T	7.9%	192.7 / 5.9	63	F594Y (TTT→TAT)	ptrA	protease 3
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCG  >  NC_000913/2956957‑2957253                                                                                                                                                    |                                                                                                                                                      tGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGaaaaa                                                                                                                                                     <  1:887143/150‑1 (MQ=255)  ggTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAAc                                                                                                                                                    >  1:1076924/1‑150 (MQ=255)      ataCCATGAATTCGCCTGGTTAAGCTGCTTCTACGTATGCGTATATCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTAt                                                                                                                                                <  1:1075602/150‑1 (MQ=255)       taCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATg                                                                                                                                               <  2:1393438/150‑1 (MQ=255)       taCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATg                                                                                                                                               <  1:75901/150‑1 (MQ=255)        aCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATg                                                                                                                                               <  2:791664/149‑1 (MQ=255)        aCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATg                                                                                                                                               >  1:791664/1‑149 (MQ=255)        aCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGc                                                                                                                                              >  2:577258/1‑150 (MQ=255)          cAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGcca                                                                                                                                            >  2:566583/1‑150 (MQ=255)               cTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGCATAAGTTTTTCCACCAAc                                                                                                                                       <  2:179403/150‑1 (MQ=255)                ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAAcc                                                                                                                                      <  2:1362211/150‑1 (MQ=255)                 tCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCg                                                                                                                                     <  2:142140/150‑1 (MQ=255)                 tCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGATGGAAAAACTTATGCCACCAACCg                                                                                                                                     <  1:54417/150‑1 (MQ=255)                   gCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGAc                                                                                                                                   <  2:1224786/150‑1 (MQ=255)                   gCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGAc                                                                                                                                   >  1:723259/1‑150 (MQ=255)                       gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCt                                                                                                                               <  1:566583/150‑1 (MQ=255)                            aGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTg                                                                                                                          >  1:550652/1‑150 (MQ=255)                             gCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGc                                                                                                                         <  2:723259/150‑1 (MQ=255)                               tGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCtt                                                                                                                       <  2:550652/150‑1 (MQ=255)                                gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTa                                                                                                                      >  1:765384/1‑150 (MQ=255)                                   cTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACt                                                                                                                   >  2:225777/1‑150 (MQ=255)                                        gTAGCGGTAAAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCa                                                                                                              <  1:1268290/150‑1 (MQ=255)                                           gCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgc                                                                                                           >  2:309617/1‑150 (MQ=255)                                                     tAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTg                                                                                                 >  1:150958/1‑150 (MQ=255)                                                           aCCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                                                           <  2:781666/150‑1 (MQ=255)                                                               cTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATc                                                                                       >  1:814711/1‑150 (MQ=255)                                                                tCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCa                                                                                      <  2:524771/150‑1 (MQ=255)                                                                  gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCAtt                                                                                    <  2:1076924/150‑1 (MQ=255)                                                                   aGCAATGCCTGGAACAGCTGCGGCAGACGCTGTGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTg                                                                                   <  2:293596/150‑1 (MQ=255)                                                                   aGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTg                                                                                   <  2:110960/150‑1 (MQ=255)                                                                    gCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGCATAAGTTTTTCCACCAAc                                                                                                                                       >  1:58867/1‑97 (MQ=38)                                                                    gCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGCATAAGTTTTTCCACCAAc                                                                                                                                       <  2:58867/97‑1 (MQ=38)                                                                      aaTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACt                                                                                                                   >  2:872113/1‑115 (MQ=255)                                                                      aaTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACt                                                                                                                   <  1:872113/115‑1 (MQ=255)                                                                             ggAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                                                         <  1:309617/150‑1 (MQ=255)                                                                                aCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCa                                                                      <  2:618043/150‑1 (MQ=255)                                                                                    cTGCCGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTg                                                                  <  2:656974/150‑1 (MQ=255)                                                                                          cAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCg                                                            >  2:940913/1‑150 (MQ=255)                                                                                          cAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCg                                                            <  1:225777/150‑1 (MQ=255)                                                                                                    gtgtAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTcc                                                  >  1:76875/1‑150 (MQ=255)                                                                                                       tAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCa                                                                             >  1:1047440/1‑120 (MQ=255)                                                                                                       tAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCa                                                                             <  2:1047440/120‑1 (MQ=255)                                                                                                             ttAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTc                                         <  2:765384/150‑1 (MQ=255)                                                                                                               aGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCgg                                       >  1:1100580/1‑150 (MQ=255)                                                                                                                gCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGa                                      >  2:1153363/1‑150 (MQ=255)                                                                                                                     aaCCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTCCGGATTACg                                 >  2:1343713/1‑150 (MQ=255)                                                                                                                       ccATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCa                               >  2:1385458/1‑150 (MQ=255)                                                                                                                        cATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCaa                              >  1:1148850/1‑150 (MQ=255)                                                                                                                         aTAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCaaa                             <  2:589262/150‑1 (MQ=255)                                                                                                                             ggCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGCGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCa                         >  2:1313241/1‑150 (MQ=255)                                                                                                                               ccGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAgg                       >  2:69414/1‑150 (MQ=255)                                                                                                                               cagtttttAGAGTTGGTGGAAAAACTTATGCCACAAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCATCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTACATGGCTTTAGGATTACGCAAAATCAgg                       >  2:643050/3‑150 (MQ=255)                                                                                                                                cGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGc                      <  2:1148850/150‑1 (MQ=255)                                                                                                                                  ttgttAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTg                    <  1:940913/150‑1 (MQ=255)                                                                                                                                     ttAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCATGCTGACa                 <  1:1385458/150‑1 (MQ=255)                                                                                                                                      tAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTa                                   <  1:295605/131‑1 (MQ=255)                                                                                                                                      tAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTa                                   >  2:295605/1‑131 (MQ=255)                                                                                                                                         ggTTGGTGGCATAAGTTTTTCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAg             >  1:470729/2‑150 (MQ=255)                                                                                                                                         ggTTGGTGGCATAAGTTTTTCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAg             >  1:465777/2‑150 (MQ=255)                                                                                                                                           ttGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCt           >  2:95234/1‑150 (MQ=255)                                                                                                                                             ggtggAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCttt         >  2:857706/1‑150 (MQ=255)                                                                                                                                                ggAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTggg      >  2:1259522/1‑150 (MQ=255)                                                                                                                                                 gAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGc     <  2:150958/150‑1 (MQ=255)                                                                                                                                                 gAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGc     <  2:1100580/150‑1 (MQ=255)                                                                                                                                                  aaaaaCTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCt    >  2:197009/1‑150 (MQ=255)                                                                                                                                                    aaaCTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCg  <  2:1147260/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCG  >  NC_000913/2956957‑2957253

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