Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	42.4%	116.4 / 18.6	67	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (35/3);  new base C (0/28);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.32e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAA  >  NC_000913/3443130‑3443416                                                                                                                                                        |                                                                                                                                            gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTAccc                                                                                                                                          >  1:55356/1‑151 (MQ=255) gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTAccc                                                                                                                                          <  1:396534/151‑1 (MQ=255) gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTAccc                                                                                                                                          >  1:1016886/1‑151 (MQ=255)   cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTAccccc                                                                                                                                        <  1:736294/151‑1 (MQ=255)       tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGcc                                                                                                                                     >  2:1066003/1‑150 (MQ=255)              acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCAc                                                                                                                             >  2:82269/1‑151 (MQ=255)                  aTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTc                                                                                                                         >  2:1073554/1‑151 (MQ=255)                   tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTc                                                                                                                         >  2:880326/1‑150 (MQ=255)                    aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTc                                                                                                                         >  2:302276/1‑149 (MQ=255)                     aaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAg                                                                                                                       >  1:63797/1‑150 (MQ=255)                       gCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCgg                                                                                                                    >  1:564073/1‑151 (MQ=255)                       gCGGTTGCCAAGGCTGAAAATACAGCGAAATTGTTGTCAGCCAGTCCCACCCAGTCCCGCCCCCG‑ACCCATGACGCGA‑TGGTCGCCGGAAAAAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTCCCCCCGCCATATTCACTTTCAGCgg                                                                                                                    <  1:955040/151‑1 (MQ=255)                         ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGGACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTTACTTTCAGCGGt                                                                                                                   >  1:585679/1‑150 (MQ=255)                            tGCCCAGGCTGAAAAAAAAGAAAATTTGTTGCCAGCAAGTCCCACCAGGCCCCGCCCCCA‑AACAAGGCCGCGA‑TGGCCCCGGAAAACAAAAAAATCCGGAAG‑GCGAAGATGGG‑CGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                  <  1:217956/148‑1 (MQ=255)                                  gggtgaaaaaACAGCGAAATTGTTCGGAACCATTTACAACC‑GTGCCCCCCCAGAAACCAGGGC‑CGATTGGTCCCCGGGAACAAAATAGTTC‑GGAA‑AGGAAGATTCCGGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTaa                                                                                                                 <  1:575631/137‑1 (MQ=255)                                     taaaaaTACAGCAAAATTTTTGGCAGCCAGTCCCAACCATCCCCCCCCCCG‑ACCCAGGCCGCGA‑TGGCCCCCGGAAACAAAAAAATGCGGAAA‑GGAAAGATGGC‑CGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                        <  1:32236/146‑1 (MQ=255)                                     tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                      >  2:522623/1‑151 (MQ=255)                                        aaaTACAGCAAATTTTTTTGCGCCCAGTCCCACCCATTCCCCCCCCCA‑ACCCAGGGCGCGG‑TGGCCGCGGGGAAAAAAATATTCCTGGAA‑GCAAAGATGGC‑CGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                   <  1:242380/151‑1 (MQ=255)                                        aaaTACAGCAAATTTGTTGTCAGCCAGTCCAAACCAGTCCCCCCCCCG‑ACCCATGACGCGA‑GGGTCGCCGGAAAAAAAATAATACTGGAA‑GCGAAGATGGC‑CGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgt                                                                                                    <  2:867460/150‑1 (MQ=255)                                             cAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                               >  2:655037/1‑150 (MQ=255)                                                 gAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                             >  2:416220/1‑148 (MQ=255)                                                  aagttgttgTCAGCCATTCCCAACCAGTCCCGCCCCCG‑AACAATGACGGGT‑TGGTCGCCGGGAACAGAATAATCCGGGAA‑GGGAAGTTTGC‑GGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                         <  2:681045/148‑1 (MQ=255)                                                     ttgttgTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg                                                                                      >  1:435790/1‑151 (MQ=255)                                                     ttgttgTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg                                                                                      >  1:249534/1‑151 (MQ=255)                                                         tgtCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                   >  2:412925/1‑150 (MQ=255)                                                         tgtCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                   >  2:218701/1‑150 (MQ=255)                                                               ccAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                            >  2:831054/1‑151 (MQ=255)                                                               ccAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                            >  1:755243/1‑151 (MQ=255)                                                                cAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                           >  1:837351/1‑151 (MQ=255)                                                                 aGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGAGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                           >  2:317895/1‑150 (MQ=255)                                                                   tcccAACCAGCACCGCCCCCG‑AACCATGACGGGA‑TGGTCGCGGGAAAAAAAATAATACTGGAA‑GCAAAGATTGC‑GGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                          <  2:910912/147‑1 (MQ=255)                                                                     ccAACCAGCCCCCCCCCCG‑AACCAGGCCGCGT‑TGGCCGCCGGAAACAAAAAAATACTGGAA‑GCGAAGATTGC‑GGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                           <  2:730319/146‑1 (MQ=255)                                                                     caacccggccccccccccA‑AACAATGACGCGA‑TGTTCGCCGGGAAAAAAATAATCTGGGAA‑GAGAAGATTGC‑CGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                        <  2:972992/136‑1 (MQ=255)                                                                         ccAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                  >  2:559130/1‑151 (MQ=255)                                                                               ccccccccA‑AACCAGGACGCGA‑TGGTCGCGGGAAAAAAAATAATGCGGGAA‑GAGAAGATTGC‑CGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAg                                                             <  2:197246/147‑1 (MQ=255)                                                                                  cccccc‑AACCATGACGCGA‑TGGGCGCCGGGGACAGAATAATACTTGAA‑GCGAAAATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGGAAAAGTGTGCTCTGTGCAGCATAGACACGAAGACCCTTCTGACGTTTCGCGGAGttt                                                          >  1:663240/1‑150 (MQ=255)                                                                                  cccccG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                          >  2:825637/1‑150 (MQ=255)                                                                                           ccATGACGCGA‑TGGTCGCCGGAAAAAAAATAATACTGGAA‑GCGAAGTTTGC‑CGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACaa                                                   <  2:662856/149‑1 (MQ=255)                                                                                            gacgACGGGT‑TGGTCGCCGGAAAAAAAATAATACGCGAA‑GCGAAGTTTGC‑CGTTTTCACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                 <  1:171017/147‑1 (MQ=255)                                                                                              tGACGCGA‑TGGTCGCGGGAAAAAAAATAATACTGAAG‑GCAAATATTGC‑GGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                  <  2:423527/147‑1 (MQ=255)                                                                                               gACGGGA‑TGGTCGCGGGGAACAGAATAATACTGGAA‑GCGAAGATGGC‑GGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACaa                                                   <  2:564173/145‑1 (MQ=255)                                                                                                aCGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                            >  1:516364/1‑151 (MQ=255)                                                                                                 cactgttGGTCGTCGGTAACACAAT‑ATGCTTGAACGTGAAGTCTGC‑CCGCATTTCCCCCGCCATAGTCACTTTCAGCGGTAAATGTGTGCTCTGTGCCGCATAGACACGCCGACCTTGCTGTCGTTTCGCGTAGTTTACCACAAAgc                                                <  1:1102421/141‑1 (MQ=255)                                                                                                    gA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                         >  1:660370/1‑150 (MQ=255)                                                                                                        ggccGCCGGAAAAAAATTAAAACTGGAA‑GCAAAGATTGC‑GGTTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                     <  1:32350/148‑1 (MQ=255)                                                                                                        ggTCGCGGGAAACAAAAAAATACTGGAG‑GCGAAGTTGGC‑GGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                       <  2:121117/149‑1 (MQ=255)                                                                                                          tCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                   >  2:193565/1‑151 (MQ=255)                                                                                                             cggggAAAAGAATAATACTGGAG‑GCAAAGATTGC‑GGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                 <  1:1062157/148‑1 (MQ=255)                                                                                                              cGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                >  2:380709/1‑150 (MQ=255)                                                                                                              cGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                >  2:1050444/1‑150 (MQ=255)                                                                                                               gggAACAAAATAATACAGGAG‑GGGAAGTTTGC‑GGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                               <  2:132650/150‑1 (MQ=255)                                                                                                               ggaaaaaaaaaaataCGGGAA‑GCAAAGATGGC‑GGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                 <  2:67530/137‑1 (MQ=255)                                                                                                                   aCAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                           >  2:425444/1‑150 (MQ=255)                                                                                                                     aaaataatACTGGAA‑GCGAATAGTGC‑CGGTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                        <  1:492614/149‑1 (MQ=255)                                                                                                                     aGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                         >  1:722090/1‑150 (MQ=255)                                                                                                                       aataatACTGGAA‑GCGAAGATTGC‑GGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                        <  1:218672/149‑1 (MQ=255)                                                                                                                       aataatACGGAAA‑GCGAAGATGGC‑GGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                      <  2:508492/151‑1 (MQ=255)                                                                                                                          aataCTGGAA‑GCGAAGATTGC‑GGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                   <  1:287286/151‑1 (MQ=255)                                                                                                                          aataCTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                    >  2:116773/1‑150 (MQ=255)                                                                                                                           ataCTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                   >  2:194788/1‑150 (MQ=255)                                                                                                                            taCTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                  >  2:484717/1‑150 (MQ=255)                                                                                                                             aCTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                  >  2:434187/1‑149 (MQ=255)                                                                                                                                ggAA‑GCGAAGATTCC‑CGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc             <  2:1041327/151‑1 (MQ=255)                                                                                                                                gaaa‑GCGAAGTGGGC‑GGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc             <  2:315602/149‑1 (MQ=255)                                                                                                                                gaaa‑GCGAAGATTGC‑CGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa               <  2:127071/147‑1 (MQ=255)                                                                                                                                   agtgagtgATTGG‑CGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg            <  1:892537/143‑1 (MQ=255)                                                                                                                                               c‑cGTGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaaa  <  2:657943/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                            GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCG‑AACCATGACGCGA‑TGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGC‑CGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAA  >  NC_000913/3443130‑3443416

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