Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,957 (G)6→7 coding (71/1386 nt) yifK → putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9511.G97.2% 255.3 / ‑0.8 71coding (65/1386 nt)yifKputative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (35/34);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑GGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAA  >  NC_000913/3980818‑3981082
                                                                                                                                       |                                                                                                                                    
gCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑ggg‑gg                                                                                                                                >  1:1032305/1‑139 (MQ=255)
       tCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCa                                                                                                                          >  2:153759/1‑139 (MQ=255)
       tCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCa                                                                                                                          >  2:1065020/1‑139 (MQ=255)
                ccATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCg                                                                                                                 >  2:803886/1‑139 (MQ=255)
                  aTAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGc                                                                                                               <  1:1065020/139‑1 (MQ=255)
                  aTAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGc                                                                                                               <  2:253133/139‑1 (MQ=255)
                  aTAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGc                                                                                                               <  1:153759/139‑1 (MQ=255)
                   tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGcc                                                                                                              >  1:818704/1‑139 (MQ=255)
                   tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGcc                                                                                                              >  2:966476/1‑139 (MQ=255)
                   tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑ggg‑gggg                                                                                                                              <  2:672678/122‑2 (MQ=255)
                   tAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑ggg‑gggg                                                                                                                              >  1:672678/1‑121 (MQ=255)
                     aCGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTg                                                                                                            >  2:1008017/1‑139 (MQ=255)
                         tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTa                                                                                                        >  2:208182/1‑139 (MQ=255)
                              ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                   <  2:151422/139‑1 (MQ=255)
                              ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                   >  2:149709/1‑139 (MQ=255)
                              ggTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATgggg                                                                                                   <  1:169041/139‑1 (MQ=255)
                               gTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATggggg                                                                                                  >  2:877769/1‑139 (MQ=255)
                               gTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATggggg                                                                                                  <  2:818704/139‑1 (MQ=255)
                                   cacCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGCTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGgcctc                                                                                              <  1:730977/139‑3 (MQ=255)
                                        gAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                                                                         >  2:1084743/1‑139 (MQ=255)
                                        gAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑ggg‑gggg                                                                                                                              >  1:1115961/1‑100 (MQ=255)
                                        gAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑ggg‑gggg                                                                                                                              <  2:1115961/101‑2 (MQ=255)
                                          aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc                                                                                       >  2:284414/1‑139 (MQ=255)
                                          aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc                                                                                       <  1:233569/139‑1 (MQ=255)
                                               tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaa                                                                                  <  2:950225/139‑1 (MQ=255)
                                                 cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                                                                <  1:208182/139‑1 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTT‑GGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  1:873978/1‑139 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGt                                                                                                           >  2:15782/1‑103 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGt                                                                                                           <  1:15782/103‑1 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  2:1032305/139‑1 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  1:877769/139‑1 (MQ=255)
                                                          cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑GGGAGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  1:738479/1‑139 (MQ=255)
                                                               ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                                                                  <  1:966476/139‑1 (MQ=255)
                                                                gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                                                                           >  1:453295/1‑113 (MQ=255)
                                                                gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                                                                           <  2:453295/113‑1 (MQ=255)
                                                                gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                                                                 >  2:87033/1‑139 (MQ=255)
                                                                     aTGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                                                            <  2:873978/139‑1 (MQ=255)
                                                                     aTGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                                                            <  1:700066/139‑1 (MQ=255)
                                                                            aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                                                                 >  1:656160/1‑95 (MQ=255)
                                                                            aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                                                                 <  2:656160/95‑1 (MQ=255)
                                                                              aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtat                                                   >  2:85159/1‑139 (MQ=255)
                                                                                cAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTatat                                                 <  1:1108322/139‑1 (MQ=255)
                                                                                cAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTatat                                                 >  1:508984/1‑139 (MQ=255)
                                                                                   aCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat                                              >  2:808596/1‑139 (MQ=255)
                                                                                     cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCg                                            <  2:508984/139‑1 (MQ=255)
                                                                                        agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  2:872953/1‑109 (MQ=255)
                                                                                        agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       >  1:193775/1‑109 (MQ=255)
                                                                                        agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  2:193775/109‑1 (MQ=255)
                                                                                        agCTACAGCGTGGGCTGGAAGCTCGACATATCCAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                                                                       <  1:872953/109‑1 (MQ=255)
                                                                                             cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGcc                                                                                                              >  2:117015/1‑65 (MQ=255)
                                                                                             cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGcc                                                                                                              <  1:117015/65‑1 (MQ=255)
                                                                                             cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTg                                    >  2:234353/1‑139 (MQ=255)
                                                                                              aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGt                                   <  1:284414/139‑1 (MQ=255)
                                                                                                cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTc                                 >  1:84128/1‑139 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCAt                                                                                                                         >  2:438744/1‑49 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCAt                                                                                                                         <  1:438744/49‑1 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGt                               <  1:149709/139‑1 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGt                               >  2:650621/1‑139 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGCCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                                                                         >  2:605384/1‑97 (MQ=255)
                                                                                                  tGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGCCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                                                                         <  1:605384/97‑1 (MQ=255)
                                                                                                    ggCTGGAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCt                             >  2:145212/1‑139 (MQ=255)
                                                                                                        ggAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCttttt                         >  2:265085/1‑139 (MQ=255)
                                                                                                         gAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtc                        <  2:738479/139‑1 (MQ=255)
                                                                                                         gAAGCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTtc                        <  1:87033/139‑1 (MQ=255)
                                                                                                            gCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                            >  2:222866/1‑68 (MQ=255)
                                                                                                            gCTCGACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                                                            <  1:222866/68‑1 (MQ=255)
                                                                                                                gACATATCGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGc                 >  2:968368/1‑139 (MQ=255)
                                                                                                                       cGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                                                             >  1:1185019/1‑72 (MQ=255)
                                                                                                                       cGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                                                             <  2:1185019/72‑1 (MQ=255)
                                                                                                                       cGAACTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAAt          >  1:911346/1‑139 (MQ=255)
                                                                                                                          aCTCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATggg       >  1:427863/1‑139 (MQ=255)
                                                                                                                            tCATCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCg     <  1:234353/139‑1 (MQ=255)
                                                                                                                              aTCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                           <  1:459443/83‑1 (MQ=255)
                                                                                                                              aTCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                                                           >  2:459443/1‑83 (MQ=255)
                                                                                                                              aTCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGaa   <  1:650621/139‑1 (MQ=255)
                                                                                                                               tCGCCC‑TGGGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGaaa  >  2:830714/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCC‑T‑GGG‑GGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAA  >  NC_000913/3980818‑3981082

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