Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,794	0	A	C	41.0%	104.0 / 17.0	61	G397G (GGT→GGG)	tap	methyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (0/25);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 3.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATC  >  NC_000913/1969644‑1969938                                                                                                                                                       |                                                                                                                                                 gcgGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCa                                                                                                                                                  <  2:898332/151‑1 (MQ=255)                gACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGc                                                                                                                                    <  1:801900/150‑1 (MQ=255)                gACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGc                                                                                                                                    >  2:10384/1‑150 (MQ=255)                   aaTATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaaaa                                                                                                                                <  1:131240/151‑1 (MQ=255)                    atatCGATCATGGTCGCGGGGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGCCTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCCCGCTACGGCTGGCAAGATTGCCCACTTCCCCTGCCACTACCGCaaaa                                                                                                                                <  2:859484/150‑1 (MQ=255)                       tCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCAc                                                                                                                            >  1:875128/1‑151 (MQ=255)                        cGATCATGGTCGCGGCGGCTTTATTCCCCAGTTCCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAAcc                                                                                                                              <  2:972379/148‑1 (MQ=255)                        cGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACg                                                                                                                           >  1:858668/1‑151 (MQ=255)                         gATCATGGTCGCGGCGGCGTTATTCCCCAGTTTCAACCCCTGCTGGACCCGATTGCCTGCCTCTTCGATCACCCTTTTGTTCTCTTTTGCTCCCTGCGCGCTACGGCTGGCAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGc                                                                                                                          <  2:103408/151‑1 (MQ=255)                            cagggTCGCGGCGGCGTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCACCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGccc                                                                                                                        <  2:663882/147‑1 (MQ=255)                               ggTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACGGACTCTTCGATCACCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGc                                                                                                                    <  2:332307/151‑1 (MQ=255)                               ggTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                     >  2:39324/1‑150 (MQ=255)                                      gcggctTTATTCCCCTGTTCCGACCCCTGCTGAACACGTTTGCCTGACTCTTCGATCAGCCCTTTGATCTCTTTTCCTGCCTGCGCGCTACGCCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                               <  2:460564/149‑1 (MQ=255)                                      gcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCgg                                                                                                             >  1:654494/1‑151 (MQ=255)                                        ggctTTATTCCCCATTTCCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCTTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                            <  2:844450/150‑1 (MQ=255)                                        ggcgTTATTCCCCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                           <  1:813531/151‑1 (MQ=255)                                        ggcgTTATTCACCAGTTCCGAACCCTGCGGGACACGTTTGACTGCCTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGGGCGCTAGGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCg                                                                                                              <  2:858061/148‑1 (MQ=255)                                           gTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                         >  1:67801/1‑150 (MQ=255)                                            ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                         >  1:451314/1‑149 (MQ=255)                                               ttCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc                                                                                                      <  1:628018/149‑1 (MQ=255)                                               ttCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCg                                                                                                    >  2:656259/1‑151 (MQ=255)                                                tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                                   >  2:421131/1‑151 (MQ=255)                                                  aCCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCtt                                                                                                 >  2:68654/1‑151 (MQ=255)                                                   ccAGTTCCGACCCCTGCTGGACACGATTGACTGCCTCTTCGACCAGCCCTTTGATCTCTTTTGCGCCCTGCGCGCTACGGCTGGAAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCg                                                                                                    <  2:603833/147‑1 (MQ=255)                                                           gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                        >  1:54903/1‑151 (MQ=255)                                                           gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                        >  1:354085/1‑151 (MQ=255)                                                            accccTGCTGAACACGTTTGACTGACTCTTCGATCAGCCCTTTGTTCTTTTTTGCTCCCGGCGCGCTACGGCTGGCAAGATTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTg                                                                                           <  1:666651/145‑1 (MQ=255)                                                             aCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                       >  1:102633/1‑150 (MQ=255)                                                              cccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                     >  2:492707/1‑151 (MQ=255)                                                                     ggacccatttGACTAACTTTTCGACCACCCTTTTGCTTTTTTTTTCTGCTGGCGCCCTAGGGCGGGCAATTTTGCGCCTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                   <  2:396788/138‑1 (MQ=255)                                                                     ggACCCGATGGACTGCCTCTTCGATCAGCCTTTTGATCTTTTTTGCGGCCTGCGCCCTACGGCGGGAAAGATTGCGCCCTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGc                                                                              <  1:957740/151‑1 (MQ=255)                                                                        cgcgtttgccggactCTCCGATCACCCCTTTGTTCTTTTTTGCTGCCTGCGCCCTACGGCGGGCAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg                                                                           <  2:903160/140‑1 (MQ=255)                                                                           gATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                        >  1:374213/1‑151 (MQ=255)                                                                               gactgactCTCCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                    <  1:733597/151‑1 (MQ=255)                                                                                        ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                            >  1:515515/1‑150 (MQ=255)                                                                                        ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                            >  1:335655/1‑150 (MQ=255)                                                                                         tCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                            >  1:307069/1‑149 (MQ=255)                                                                                          cgaccacccttttactcttTTTTGCGGCCCCCCCGCTAGGGCTGGAAAGTTTCCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                          <  2:43344/135‑1 (MQ=21)                                                                                              cAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                                     >  1:164560/1‑151 (MQ=255)                                                                                                gCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                                   >  2:993152/1‑151 (MQ=255)                                                                                                  ccTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGt                                                 >  1:569923/1‑151 (MQ=255)                                                                                                      tGATCCCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTACGCACTTCACCTTCCACTACCGCAAAACCACACCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGTCGCATACAGGGCCAGAATATTGGACTGGAAAGCAATTCCGt                                                 >  2:592985/1‑147 (MQ=25)                                                                                                         tctcTTTTGCTGCCGGCGCGCTACGGCTGGAAAGATTGCGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                             <  1:912453/148‑1 (MQ=255)                                                                                                          ctctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                         >  1:135015/1‑151 (MQ=255)                                                                                                             ttttGCTCCCTGCGCGCTCCGGCGGGCAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                         <  2:202963/148‑1 (MQ=37)                                                                                                               ttGCTCCCTGCGCGCACCGGCGGGCAAGTTGCCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa                                     <  2:248606/150‑1 (MQ=25)                                                                                                                  ctgccggcGCGCTGGGGCGGGCAAGTTTGCTCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                                   <  1:304822/149‑1 (MQ=25)                                                                                                                   tgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                                >  2:919875/1‑151 (MQ=255)                                                                                                                   tgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                                >  1:271963/1‑151 (MQ=255)                                                                                                                       tgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                            >  1:657401/1‑151 (MQ=255)                                                                                                                           cgcTCCGGCTGCAAAGATTCCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAt                        <  1:422498/151‑1 (MQ=25)                                                                                                                           cgcTCCGGCGGGAAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCaa                         <  1:10384/150‑1 (MQ=255)                                                                                                                             cTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttt                      >  2:579021/1‑151 (MQ=255)                                                                                                                              tACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                     >  1:578265/1‑151 (MQ=255)                                                                                                                                 ggCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                   <  1:453853/150‑1 (MQ=255)                                                                                                                                 ggCTGGCAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttt                      <  1:737143/147‑1 (MQ=255)                                                                                                                                 gccgggCAAGTTTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTg                  <  2:395053/147‑1 (MQ=25)                                                                                                                                   cTGGCAAGATTCCGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtt                       <  1:278018/144‑1 (MQ=255)                                                                                                                                       cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt            >  1:917869/1‑151 (MQ=255)                                                                                                                                            ttccgcACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                <  1:1804/139‑1 (MQ=255)                                                                                                                                                 cACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc  <  1:988218/151‑1 (MQ=255)                                                                                                                                                 cACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc  <  1:307380/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GCGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATC  >  NC_000913/1969644‑1969938

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