Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,907 +GG intergenic (+109/‑11) pldB → / → b4555 lysophospholipase L(2)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,9071.G92.4% 1490.0 / 97.5 411intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (11/20);  new base G (194/185);  total (206/205)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_0009134,007,9072.G92.5% 1476.4 / 97.5 412intergenic (+109/‑11)pldB/b4555lysophospholipase L(2)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (11/20);  new base G (195/186);  total (206/206)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 466 total aligned reads displayed.

CAGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCTCCA  >  NC_000913/4007812‑4007999
                                                                                                    ||                                                                                            
cAGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgc                                                                                                <  1:13936788/96‑1 (MQ=255)
 aGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgcg                                                                                               >  1:9170417/1‑95 (MQ=255)
 aGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgcg                                                                                               >  2:13454191/1‑95 (MQ=255)
 aGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgcg                                                                                               <  2:7650713/96‑2 (MQ=255)
 aGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgcg                                                                                               >  2:6319483/1‑95 (MQ=255)
    ttgttgCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCg‑‑cg‑‑cgcggca                                                                                            >  1:10647694/1‑92 (MQ=255)
               gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTtgtg                                                                                 <  2:6522727/96‑1 (MQ=255)
               gATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTtgtg                                                                                 >  2:3812326/1‑96 (MQ=255)
                aTTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTgtgt                                                                                >  1:5757387/1‑96 (MQ=255)
                aTTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTgtgt                                                                                >  1:677048/1‑96 (MQ=255)
                aTTTAGATGGCACGTTACTTTCTCCCCACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CGCACGG‑‑CATCAACTTTgtgt                                                                                >  2:11067590/1‑96 (MQ=255)
                 tttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGtt                                                                               <  1:6906175/96‑1 (MQ=255)
                  ttAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGttt                                                                              <  2:510673/96‑1 (MQ=255)
                   tAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTg                                                                             >  2:2012610/1‑96 (MQ=255)
                     gATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCg                                                                           <  1:7435694/96‑1 (MQ=255)
                      aTGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCGAAGAAACTCTGAAGCTGCTCACCGCACG‑‑CGG‑‑CATCAACTTTTTGTTTGCgt                                                                          >  2:9550185/1‑95 (MQ=255)
                      aTGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGa                                                                          <  1:3291766/96‑1 (MQ=255)
                       tGGCACGTTACTTTCTCCCGACCCTACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGAc                                                                         <  2:10146201/96‑1 (MQ=255)
                        ggCTCGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGAcc                                                                        <  1:3019515/96‑1 (MQ=255)
                         gCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGAGACCg                                                                       >  1:10035831/1‑96 (MQ=255)
                          cACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGGTTGCGACCgg                                                                      >  2:13214107/1‑96 (MQ=255)
                           aCGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTCCGACCGgt                                                                     <  1:10490531/96‑1 (MQ=255)
                                aCTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTca                                                                <  1:13477494/96‑1 (MQ=255)
                                 cTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                               <  2:1173964/96‑1 (MQ=255)
                                 cTTCCTCCCGCCCACCCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcac                                                               <  1:1587248/96‑1 (MQ=255)
                                  tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacc                                                              <  2:11961259/96‑1 (MQ=255)
                                  tttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACAGGTCGTcacc                                                              >  2:10737264/1‑96 (MQ=255)
                                   ttCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacca                                                             >  1:3614804/1‑96 (MQ=255)
                                    tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcacctc                                                            <  1:3418973/96‑3 (MQ=255)
                                    tctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTcaccac                                                            <  2:8040438/96‑1 (MQ=255)
                                     ctcCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACg                                                           >  1:11978648/1‑96 (MQ=255)
                                       cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                                         <  2:362301/96‑1 (MQ=255)
                                       cccGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGtt                                                         <  1:6431098/96‑1 (MQ=255)
                                         cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGa                                                       <  2:12841142/96‑1 (MQ=255)
                                         cGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGa                                                       >  2:9696010/1‑96 (MQ=255)
                                          gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTgtt                                                      <  1:11080359/96‑3 (MQ=255)
                                          gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                      >  2:7468423/1‑96 (MQ=255)
                                          gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                      >  1:12558090/1‑96 (MQ=255)
                                          gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                      <  1:13698000/96‑1 (MQ=255)
                                          gACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAt                                                      >  1:2880772/1‑96 (MQ=255)
                                           aCCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtg                                                     <  1:4555227/96‑1 (MQ=255)
                                            ccATACGTTATCCCCTTACGCCAAAGAAACTCTTAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTCTGTGTATGCGACCGGTAGTCACCACGTTGAtgt                                                    >  1:14059528/1‑96 (MQ=255)
                                            ccATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCGGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgt                                                    >  2:12633350/1‑96 (MQ=255)
                                              acaCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgttg                                                  <  1:7640426/94‑3 (MQ=255)
                                               tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGAtgtgcg                                                 >  1:948533/1‑94 (MQ=255)
                                               tACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTggg                                                 >  2:8916379/1‑96 (MQ=255)
                                                aCGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTgggg                                                <  2:6999684/96‑1 (MQ=255)
                                                  gTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCa                                              >  2:4855583/1‑96 (MQ=255)
                                                   ttATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaa                                             >  1:423269/1‑96 (MQ=255)
                                                    tATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCaaa                                            >  2:11373351/1‑96 (MQ=255)
                                                     aTCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAt                                           >  2:14020539/1‑96 (MQ=255)
                                                       ccccTTTCGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                         <  1:11901876/96‑1 (MQ=255)
                                                       ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAttt                                         <  1:11332976/96‑2 (MQ=255)
                                                       ccccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTc                                         <  1:3965248/96‑1 (MQ=255)
                                                       ccccTGTCGCCAACGCCACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAAtcc                                         <  1:14337827/96‑3 (MQ=255)
                                                        cccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcg                                        <  1:707952/96‑1 (MQ=255)
                                                         ccTTACGCCAAAGAAACTCTGAATCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                       >  2:5879437/1‑96 (MQ=255)
                                                         ccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTGGTCACCACGGTGATGGGGGGCAAATTcgc                                       >  1:6654952/1‑96 (MQ=255)
                                                         ccTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgc                                       >  1:8559763/1‑96 (MQ=255)
                                                          cTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgcg                                      <  2:2619490/96‑1 (MQ=255)
                                                          cTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTcgcg                                      >  1:12322243/1‑96 (MQ=255)
                                                            tACGCCAAAGAAACTCTGAAGCTGCTCCCCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGAt                                    <  1:2450079/96‑1 (MQ=255)
                                                             aCGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATa                                   >  2:3371697/1‑96 (MQ=255)
                                                              cGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATaa                                  <  2:10456782/96‑1 (MQ=255)
                                                               gCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAAt                                 >  2:4425513/1‑96 (MQ=255)
                                                                ccAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGGCCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGatcatc                                >  1:1303894/1‑96 (MQ=255)
                                                                 cAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                               >  1:11187212/1‑96 (MQ=255)
                                                                 cAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCt                               <  1:5656897/96‑1 (MQ=255)
                                                                  aaaGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTg                              >  1:2255414/1‑96 (MQ=255)
                                                                    aGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                            >  2:561310/1‑96 (MQ=255)
                                                                    aGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                            <  1:8042831/96‑1 (MQ=255)
                                                                    aGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGga                            >  2:3097664/1‑96 (MQ=255)
                                                                      aaaCTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGgaga                          >  1:1710096/1‑96 (MQ=255)
                                                                       aaCTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAt                         <  1:4470007/96‑1 (MQ=255)
                                                                        aCTCTGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGAtt                        <  1:7575301/96‑1 (MQ=255)
                                                                         ctctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                       <  2:232194/96‑1 (MQ=255)
                                                                         ctctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTa                       <  2:596874/96‑1 (MQ=255)
                                                                          tctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                      >  2:1498281/1‑96 (MQ=255)
                                                                          tctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                      >  2:14250057/1‑96 (MQ=255)
                                                                          tctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                      >  2:12176886/1‑96 (MQ=255)
                                                                          tctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTaa                      >  2:7255294/1‑96 (MQ=255)
                                                                           ctGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTacg                     >  1:12164478/1‑94 (MQ=255)
                                                                            tGAAGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGt                    >  2:9107312/1‑96 (MQ=255)
                                                                               aGCTGCTCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCtt                 <  2:4081057/96‑1 (MQ=255)
                                                                                gctgctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTa                <  2:774157/96‑1 (MQ=255)
                                                                                 ctgctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTAc               >  2:8166636/1‑96 (MQ=255)
                                                                                  tgctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACa              >  1:12676745/1‑96 (MQ=255)
                                                                                  tgctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTATTACa              >  2:10259756/1‑96 (MQ=255)
                                                                                   gctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACAt             <  1:4831350/96‑1 (MQ=255)
                                                                                   gctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACAt             >  2:8772318/1‑96 (MQ=255)
                                                                                   gctCACCG‑‑CG‑‑CGCGGCATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACAt             >  1:6192822/1‑96 (MQ=255)
                                                                                      ccgcg‑‑cg‑‑cgG‑‑CATCACCTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTa        <  1:9597363/93‑1 (MQ=255)
                                                                                      ccgcg‑‑cg‑‑cgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTa        >  2:8606432/4‑96 (MQ=255)
                                                                                       cgcg‑‑cg‑‑cgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc       <  2:13835113/94‑1 (MQ=255)
                                                                                       cgcg‑‑cg‑‑cgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc       <  2:13641310/94‑1 (MQ=255)
                                                                                       cgcg‑‑cg‑‑cgG‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc       <  1:13350197/94‑1 (MQ=255)
                                                                                       cgag‑‑cg‑‑cgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAc       <  1:3880492/93‑1 (MQ=255)
                                                                                          g‑‑cg‑‑cgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTAccac    <  1:11505142/96‑3 (MQ=255)
                                                                                          g‑‑cg‑‑cgg‑‑CATCAACTTTATGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCTc    <  2:8318609/96‑1 (MQ=255)
                                                                                              g‑‑cgg‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCTCCa  >  1:8717326/1‑96 (MQ=255)
                                                                                                    ||                                                                                            
CAGGTTGTTGCGTCTGATTTAGATGGCACGTTACTTTCTCCCGACCATACGTTATCCCCTTACGCCAAAGAAACTCTGAAGCTGCTCACCG‑‑CG‑‑CGC‑‑CATCAACTTTGTGTTTGCGACCGGTCGTCACCACGTTGATGTGGGGCAAATTCGCGATAATCTGGAGATTAAGTCTTACATGATTACCTCCA  >  NC_000913/4007812‑4007999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: