Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,957,110	A→T	7.7%	I592K (ATA→AAA)	ptrA ←	protease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,957,110	0	A	T	7.7%	198.8 / 5.9	65	I592K (ATA→AAA)	ptrA	protease 3
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (3/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGC  >  NC_000913/2956962‑2957258                                                                                                                                                     |                                                                                                                                                     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                                                                                                                                            <  2:58867/97‑1 (MQ=38)                                                               gCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGCATAAGTTTTTCCACCAAc                                                                                                                                            >  1:58867/1‑97 (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                                        >  2:295605/1‑131 (MQ=255)                                                                                                                                 tAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTa                                        <  1:295605/131‑1 (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)                                                                                                                                                    tATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGc  <  1:1153363/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGC  >  NC_000913/2956962‑2957258

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