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

CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAA  >  NC_000913/3443128‑3443416                                                                                                                                                        |                                                                                                                                             cAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑c                                                                                                                                            <  2:1056726/151‑1 (MQ=255) cAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑c                                                                                                                                            >  2:213740/1‑151 (MQ=255)   gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑ccc                                                                                                                                          >  1:835062/1‑151 (MQ=255)    aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑cccc                                                                                                                                         >  2:162520/1‑151 (MQ=255)         tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGcc                                                                                                                                     >  1:339755/1‑150 (MQ=255)          agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCat                                                                                                                                   >  1:445764/1‑151 (MQ=255)          agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCat                                                                                                                                   >  2:768929/1‑151 (MQ=255)            agTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCatat                                                                                                                                 >  2:729206/1‑151 (MQ=255)               aactgagaaaGCGGTTGCTCAGGGTGCAAGAACAGCGAACTTGTTGTCACCCTGTTCCAACCAGAACCGCGCCCGAACC‑AGGGCGCGATGGTCGCCGGGAAC‑AGAATAATTCTGGAACCGAAGAAGGCCGCGATTG‑CCCCCGCCa                                                                                                                                    <  2:671769/138‑1 (MQ=255)                       aaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGc                                                                                                                      >  1:743184/1‑151 (MQ=255)                           ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGt                                                                                                                   >  2:758705/1‑150 (MQ=255)                            gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                  >  2:779844/1‑150 (MQ=255)                                     gCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                        >  1:261974/1‑151 (MQ=255)                                      cTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                       >  1:606256/1‑151 (MQ=255)                                          aaatacaaatattttgtttTCAGCCCGTCCCAACTATTCCCGCCCCAAAACA‑ATGCCGTGTTGGTCGCCGCAAAC‑AAAAAATTACAAGAGGGGAAGATGGTGGGGTTCC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                     <  2:836308/136‑1 (MQ=255)                                          aaaTACAGCAAAATTGTTGCCAGCCAGTCCCACCCAGTCCCCCCCCCAAACC‑ATGACGGGATGGTCGCGGGAAAC‑AAAATAATACTGGAAGCGAAGATTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                      <  1:688742/148‑1 (MQ=255)                                            aaaCAGCAAAATTTTTTTCGCCAAGTTCAAACAAG‑ACCCCCCCCCAACCAAGGACGCGATGGCCCCGGGAAAC‑GAAAAAATAGTGAAGGAAAAGAGTGCCGGTTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgca                                                                                                 <  2:682816/149‑1 (MQ=255)                                            aTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                  >  1:874306/1‑150 (MQ=255)                                                   aaatttgtttTCGCCCAGTCCAAACCAGTCCCCCCCCGAAACC‑GTGCGGTGTGGGTCGGCGGAAAA‑AAAATAATACGGGAAGAAAAAATGGCCGGATTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                             <  2:992662/144‑1 (MQ=255)                                                    aatttgttgTCAGCCCTTCCCACCCAGTCCCCCCCCCAACCC‑AGGACGCGATGGTCGCCGGAAAC‑AAAATAATACTGGAAGCGAAGTTTGCCGGGTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAg                                                                                           <  1:561682/146‑1 (MQ=255)                                                    aaaTTTTTGTCAGCCTTTCCCAACCAGTCCCCCCCCCAAACC‑AGGACGCGAGGGGCGCGGGGAAC‑AAAATAATACTGGAAGGGAAGATGGCCGGGTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                          <  1:430778/150‑1 (MQ=255)                                                     aaTTGTTGTCCCCCAGTTCCAACCATTACCGCCCCCGAACC‑ATGGCGCGGTGGTCGCCGAAAAC‑AAAATAATACTGAAAGCGAAAATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGaca                                                                                        <  2:81933/151‑1 (MQ=255)                                                        ttttcccaccctgTTCCAACCAGTCCCCCCCCCAACCA‑AGGACGGGGGGGCCGCCGGGAAA‑AAAAAAATTCTGGAGGAGAAGTTTGCCGGTTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAg                                                                                           <  2:1013935/133‑1 (MQ=255)                                                        tgtggtCAGCCGTTTCCACCCGGCCCCCCCCGCAACCC‑ATGGCGCGAGGCTCGCGGGAAAA‑AGAAAAACAGTGGGAGAGAATATTGGGGGGTTCC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                          <  2:928507/146‑1 (MQ=255)                                                               aGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGACGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                              >  2:287/1‑151 (MQ=255)                                                                 ccAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                            >  2:776046/1‑151 (MQ=255)                                                                   aGTCCCACCCAGCCCCGCCCCCGAACC‑AGGACGGGATGGTCGCGGGAAAA‑AAAAAAATACGGGAAGGGAAGATTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                           <  1:250299/150‑1 (MQ=255)                                                                    gTTCCAACCATTCCCCCCCCCGAACC‑AGGACGGGATGTTCGCGGGGACC‑AAAATAATGCTGGAGGCGAAGATCGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGGCCTTGCt                                                                           <  1:313829/149‑1 (MQ=255)                                                                     tcccACCCATTCCCCCCCCCGACCC‑ATGACGCGATGGTCGCGGGAAAA‑AAAATAATACTGAAGGCAAAGTTGGCCGGGATCC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                         <  1:608321/148‑1 (MQ=255)                                                                        cAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                     >  1:683791/1‑151 (MQ=255)                                                                        cAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                     >  1:574985/1‑151 (MQ=255)                                                                         acccAGTCCCGCCCCCGAACC‑ATGACGCGATGGTCGCGGGGAAA‑AGAATAATACTGGAAGCGAAGATGGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                        <  1:1030793/145‑1 (MQ=255)                                                                           ccAGTCCCCCCCCCGAACC‑AGGCCGCGAGGGTCGCGGGAAAA‑AAAATAATACGGGAAGGAAAGATTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                  <  1:254427/151‑1 (MQ=255)                                                                                 ccccccccAAACC‑ATGACGCGAGGGTCGCGGGAAAC‑AAAATAATACGGGAAGCAAAGATTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                               <  1:968609/145‑1 (MQ=255)                                                                                    cccccGAACC‑AGGCCGCGTGGGTCGCGGGAAAA‑AAAATAATCCTGAAAGCGAAGTTTGCCGGGTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                           <  1:963965/149‑1 (MQ=255)                                                                                       ccGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                       >  2:372682/1‑150 (MQ=255)                                                                                          accc‑AGGACGCGATGGTCGCGGGGAAC‑AAAATAATACTGAAGGGGAAAATTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                                      <  1:1009070/146‑1 (MQ=255)                                                                                          aaCC‑ATGACGCGATGGTCGCCGGAAAC‑AGAATAATACTGGAAGCGAAGTTTGCGGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                     <  1:1042775/149‑1 (MQ=255)                                                                                           acccgtGACGCGATGGTCCCCGGAAAA‑AAAATAATAGGGGAAGAAAAGATTGCGGGTTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                    <  2:257479/145‑1 (MQ=255)                                                                                           aCC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                  >  1:786119/1‑151 (MQ=255)                                                                                                   cgcgATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                           >  1:985634/1‑151 (MQ=255)                                                                                                    gcgAGGGTCGCGGGAAAA‑AAAATAATACTGGAAGGAAAGATTGCGGGAATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                            <  1:791845/149‑1 (MQ=255)                                                                                                      gttGGTCGCGGGGAAA‑AAAATAATACGGGAGGCAATGTGTGGCGGGTTCC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                            <  2:1033072/145‑1 (MQ=255)                                                                                                        tGGTCGACGGGTAC‑AGGATAATACTGGAAGCGCAGATTGCCGGGATTA‑CCCGCGCCATATTCACTTTCAGCGGGAAATGTGTGCTCGGTGCGGCATAGACACGACGACCTTGCTGACGTTTAGCGTAGCTTACCACAGTCCGGCGTTGatc                                      >  2:964388/1‑149 (MQ=255)                                                                                                         gttcaccGAGAAA‑AAAATAATACGGAAAGCGAGGATGGCTGGGTCTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                          <  2:370555/141‑1 (MQ=255)                                                                                                          ggccccGGGAAAAAGAAAAATATTGGAAGCGAGGATTGCCGGTTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                       <  2:193658/145‑1 (MQ=255)                                                                                                            cGCCGGGCAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACg                                   >  2:1013878/1‑150 (MQ=255)                                                                                                             gCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCt                                 >  1:43695/1‑151 (MQ=255)                                                                                                                 ggaaaaaaaataatACTGAAAGCGAAATTTCCCGGTATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaa                              <  1:866317/143‑1 (MQ=255)                                                                                                                   aaC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                            >  2:629587/1‑150 (MQ=255)                                                                                                                     c‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                         >  2:137/1‑151 (MQ=255)                                                                                                                     c‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                         >  1:568383/1‑151 (MQ=255)                                                                                                                     aaaaaaaaataCTGGACGCGAATATTG‑CGGGTTTACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                          <  2:96603/144‑1 (MQ=255)                                                                                                                              taCGGGAAGAGAAGATTGCGGGAATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                  <  1:420368/150‑1 (MQ=255)                                                                                                                                  ggaggaaaaatttcccGGGTTTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa               <  1:813527/135‑1 (MQ=255)                                                                                                                                     aGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc          >  2:870027/1‑151 (MQ=255)                                                                                                                                     aGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc          >  2:196997/1‑151 (MQ=255)                                                                                                                                       cGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc        >  2:502327/1‑151 (MQ=255)                                                                                                                                        gAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt       >  2:355472/1‑151 (MQ=255)                                                                                                                                         aaGATTGCCGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa         <  1:401612/148‑1 (MQ=255)                                                                                                                                             ttGCCGGTATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAACATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa   <  1:158447/150‑1 (MQ=255)                                                                                                                                             ttGCCGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa    <  1:408775/149‑1 (MQ=255)                                                                                                                                             ttGCCGGGATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa    <  1:968929/149‑1 (MQ=255)                                                                                                                                              tgcggtaatca‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGGCGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCCCGCTCAACAAATACAACAAAGAACGTCACTGCaaa  <  2:800186/140‑1 (MQ=255)                                                                                                                                               gcggggATTC‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaaa  <  1:291308/146‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                             CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACC‑ATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAA  >  NC_000913/3443128‑3443416

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