Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,023,018	G→T	6.1%	S289* (TCA→TAA)	yccS ←	putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,023,018	0	G	T	6.1%	240.9 / 2.5	65	S289* (TCA→TAA)	yccS	putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base T (2/2);  total (29/36)
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.82e-01

AGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGA  >  NC_000913/1022880‑1023156                                                                                                                                           |                                                                                                                                           aGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGa                                                                                                                                           <  1:487693/140‑1 (MQ=255)           ccAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTgg                                                                                                                                 <  1:1610547/140‑1 (MQ=255)            cAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGc                                                                                                                                <  2:1101825/140‑1 (MQ=255)            cAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGc                                                                                                                                <  2:1101654/140‑1 (MQ=255)                    tGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGGGCGCTCAAAATGCGGAACATGTTGTTATGGCTGACGCAACAATATACAGCGTGACAGTTGCTGGCACGCCTgg                                                                                                                        >  2:1210914/1‑140 (MQ=255)                           aTCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTTACAGTTGCTGGCACGCCTGGCCCTGCa                                                                                                                 >  2:599654/1‑140 (MQ=255)                           aTCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTTACAGTTGCTGGCACGCCTGGCCCTGCa                                                                                                                 >  2:1586252/1‑140 (MQ=255)                                             aTCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGaa                                                                                               <  1:785753/140‑1 (MQ=255)                                               cgcgCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACACAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGaaaa                                                                                             <  2:1474469/140‑1 (MQ=255)                                               cgcgCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGaaaa                                                                                             <  2:359596/140‑1 (MQ=255)                                                 cgcATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACg                                                                                           <  1:320925/140‑1 (MQ=255)                                                   cATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGa                                                                                         >  1:620825/1‑140 (MQ=255)                                                   cATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCAGCGACATCAGCCGCTGAAAACGGa                                                                                         >  1:620811/1‑140 (MQ=255)                                                    aTCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGaa                                                                                        <  2:877318/140‑1 (MQ=255)                                                     tCCGCTCCAGCGCATCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGATGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAAc                                                                                       <  2:256163/140‑1 (MQ=255)                                                     tCCGCTCCAGCGCAGCATCATTATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATAGACATCAGCCGCTGAAAACGGAAc                                                                                       <  2:256126/140‑1 (MQ=255)                                                     tCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAAc                                                                                       >  1:884077/1‑140 (MQ=255)                                                      cagCTCCAGCGCAGCATCAATAGGCGTAAAAGCGCGCTCAAAATGCGGATCATGTAGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACa                                                                                      <  2:1069107/138‑1 (MQ=255)                                                                  agcaTCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCtgt                                                                          >  2:1303839/1‑140 (MQ=255)                                                                    caTCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTgg                                                                        >  2:807474/1‑140 (MQ=255)                                                                              cGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGtt                                                              >  1:1666894/1‑140 (MQ=255)                                                                              cGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGtt                                                              >  2:65273/1‑140 (MQ=255)                                                                              cGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGtt                                                              >  1:488917/1‑140 (MQ=255)                                                                                tAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCa                                                            <  2:1468003/140‑1 (MQ=255)                                                                                  aaaGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACg                                                          <  1:90276/140‑1 (MQ=255)                                                                                  aaaGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACg                                                          <  1:1551158/140‑1 (MQ=255)                                                                                  aaaGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACg                                                          <  1:110739/140‑1 (MQ=255)                                                                                   aaGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGc                                                         <  1:1641406/140‑1 (MQ=255)                                                                                    aGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCa                                                        <  1:1475117/140‑1 (MQ=255)                                                                                    aGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCa                                                        <  1:126244/140‑1 (MQ=255)                                                                                        cgcTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGt                                                    >  1:742214/1‑140 (MQ=255)                                                                                        cgcTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGt                                                    >  1:807641/1‑140 (MQ=255)                                                                                        cgcTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGt                                                    >  1:812141/1‑140 (MQ=255)                                                                                        cgcTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGt                                                    >  1:1174647/1‑140 (MQ=255)                                                                                            cAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGa                                                <  2:884077/140‑1 (MQ=255)                                                                                              aaaTGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTTACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATa                                              <  1:599654/140‑1 (MQ=255)                                                                                              aaaTGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTTACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATa                                              <  1:1586252/140‑1 (MQ=255)                                                                                                   cGGATCATGTTGATAAGGCTGACTCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGaa                                         <  2:908288/140‑1 (MQ=255)                                                                                                         aTGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATgag                                   <  2:1556523/140‑1 (MQ=255)                                                                                                              gATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAagag                              <  2:735844/140‑1 (MQ=255)                                                                                                                    gCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCa                        >  1:1370340/1‑140 (MQ=255)                                                                                                                    gCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCa                        <  2:277574/140‑1 (MQ=255)                                                                                                                      tGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACg                      <  1:1282127/140‑1 (MQ=255)                                                                                                                      tGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACg                      <  1:1632823/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  1:1658943/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    >  1:592768/1‑140 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  2:1141454/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  1:1453559/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  1:305810/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  1:191035/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  2:1539588/140‑1 (MQ=255)                                                                                                                        aCGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCt                    <  1:1594102/140‑1 (MQ=255)                                                                                                                         cGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAAAATGAGAAGAGCTGGCACGCTc                   >  2:1458874/1‑140 (MQ=255)                                                                                                                              caaAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCTCGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGaa              >  2:530306/1‑140 (MQ=255)                                                                                                                                aaaTACACCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTTAAAACGGAACAGCACGTAGCTGTGGCGAAAATGTTCACGCAATGTTAGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAata            <  1:190460/140‑1 (MQ=255)                                                                                                                                      aGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCtgt      >  1:1236926/1‑140 (MQ=255)                                                                                                                                        cctGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGACGCTGAGGCGAAAATGTTCACGCAATGTTTGATACTGAAAATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  1:301181/3‑140 (MQ=255)                                                                                                                                        cGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGATAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  1:704041/1‑140 (MQ=255)                                                                                                                                        cGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  2:1582846/1‑140 (MQ=255)                                                                                                                                        cGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  1:360442/1‑140 (MQ=255)                                                                                                                                        cGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  1:390115/1‑140 (MQ=255)                                                                                                                                        cGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGc    >  1:655279/1‑140 (MQ=255)                                                                                                                                          tGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGa  >  1:749796/1‑140 (MQ=255)                                                                                                                                          tGACAGTTGCTGGAACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGa  >  1:1674412/1‑140 (MQ=255)                                                                                                                                          tGACAGTTGCTGGAACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGa  >  1:1674413/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                           AGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGA  >  NC_000913/1022880‑1023156

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