Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,860 Δ2 bp intergenic (+59/‑24) b3836 → / → b3838 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,019,8600C.99.8% 2251.0 / ‑3.8 485intergenic (+59/‑25)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (272/212);  total (273/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01
*NC_0009134,019,8610C.99.8% 2255.0 / ‑4.0 485intergenic (+60/‑24)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (272/212);  minor base G (1/0);  total (273/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 500 total aligned reads displayed.

AAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGG  >  NC_000913/4019775‑4019942
                                                                                     ||                                                                                 
aaGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTgggg                                                                           >  1:7416322/1‑93 (MQ=255)
 aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          >  2:806365/1‑93 (MQ=255)
  gcgcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                                         >  1:7094295/1‑93 (MQ=255)
    gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       >  1:265754/1‑93 (MQ=255)
    gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       <  1:4233438/93‑1 (MQ=255)
    gacACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       <  1:2059374/91‑1 (MQ=255)
     ccACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCa                                                                      >  2:8806692/1‑93 (MQ=255)
     ccACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCa                                                                      >  1:6660881/1‑93 (MQ=255)
      cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     <  1:6895898/93‑1 (MQ=255)
       aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    >  2:8399604/1‑93 (MQ=255)
       aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    >  1:8513305/1‑93 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  <  2:4980963/93‑1 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  1:8648153/1‑93 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  1:5262630/1‑93 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCAGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  <  2:1260493/93‑1 (MQ=255)
          aTAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                                                 <  1:4979152/93‑1 (MQ=255)
           tAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                                                <  1:3143017/93‑1 (MQ=255)
           tAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                                                >  2:4804774/1‑93 (MQ=255)
            aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACtc                                                               >  1:6440641/1‑92 (MQ=255)
             aaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                                              >  2:1725826/1‑93 (MQ=255)
             aaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                                              >  1:4504147/1‑93 (MQ=255)
              agagCAGGTGTAATCCGTGTTTGATATCGGTTGTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGcc                                                             <  1:5495250/93‑1 (MQ=255)
                 gCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                                          >  1:2826032/1‑93 (MQ=255)
                 gCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                                          >  2:8186837/1‑93 (MQ=255)
                 gCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                                          <  2:1639141/93‑1 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        <  2:7556009/93‑1 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  1:4174304/1‑93 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  2:5219746/1‑93 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  2:1472714/1‑93 (MQ=255)
                     gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTGTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                                                      <  2:2366980/93‑1 (MQ=255)
                     gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                                                      >  1:3024929/1‑93 (MQ=255)
                     gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                                                      >  1:3592202/1‑93 (MQ=255)
                      tgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCATGTGGCgg                                                     >  1:7618798/1‑93 (MQ=255)
                       gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGt                                                    <  1:5804051/93‑1 (MQ=255)
                       gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGt                                                    >  1:387780/1‑93 (MQ=255)
                        tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                                   >  2:3619067/1‑93 (MQ=255)
                         aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  >  2:3168757/1‑93 (MQ=255)
                         aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  >  1:7889511/1‑93 (MQ=255)
                          aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaa                                                 >  1:895070/1‑93 (MQ=255)
                          aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCAGCCT‑‑CTGCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGCaaa                                                 >  2:5495352/1‑93 (MQ=255)
                           tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                <  2:6013807/93‑1 (MQ=255)
                           tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  2:5367156/1‑93 (MQ=255)
                           tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  1:4399954/1‑93 (MQ=255)
                            ctgtgtgtGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaaa                                               <  2:8654535/91‑1 (MQ=255)
                             cGTGTTTGATGTCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              <  1:226464/93‑1 (MQ=255)
                             cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              <  1:6999996/93‑1 (MQ=255)
                             cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              <  2:4614254/93‑1 (MQ=255)
                              gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                                             <  2:4359505/93‑1 (MQ=255)
                              gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGCCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                                             <  1:1856588/93‑1 (MQ=255)
                               tgtTTGATATCGGTTTTAGCGAACCGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGCCTGCCTGTGGCGGTAAAAACgg                                            >  1:7523087/1‑93 (MQ=255)
                                gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGt                                           <  2:6406448/93‑1 (MQ=255)
                                 tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          <  2:136172/93‑1 (MQ=255)
                                 tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          >  1:3283554/1‑93 (MQ=255)
                                 gttGATATCGGTTTTAGTTAACTGCTATTGCTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCGGTGTCGGTAAAAACGGTa                                          <  1:4844833/92‑1 (MQ=255)
                                  ttGATATCGGTTTTAGCTAGCTGTTGTTGGTGTTCATCATCGCGCTCGTAG‑‑TGCTGGGGCCGCAACGACTGCCTGTGTCGGTAAAAACGGTAg                                         <  2:4732303/93‑1 (MQ=255)
                                  ttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAg                                         >  1:6148258/1‑93 (MQ=255)
                                   tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAgg                                        >  2:2247618/1‑92 (MQ=255)
                                   tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                                        <  1:6057633/93‑1 (MQ=255)
                                     atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCgg                                      >  2:4038436/1‑93 (MQ=255)
                                      tatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCggg                                     >  2:3721447/1‑93 (MQ=255)
                                       atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAACGGTAGCGGGc                                    >  1:1186948/1‑93 (MQ=255)
                                         cGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTg                                  <  2:2184319/93‑1 (MQ=255)
                                         cGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTg                                  <  1:5675405/93‑1 (MQ=255)
                                           gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTggt                                >  1:6829241/1‑92 (MQ=255)
                                           gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGa                                >  2:2508853/1‑93 (MQ=255)
                                           gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGa                                <  2:2856212/93‑1 (MQ=255)
                                           gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGa                                >  1:730388/1‑93 (MQ=255)
                                           gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAAAGGTAGCGGGCTGGa                                >  2:5885382/1‑93 (MQ=255)
                                             tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt                              >  1:1017259/1‑93 (MQ=255)
                                             tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt                              >  2:2983524/1‑93 (MQ=255)
                                                aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGGTTcgc                           >  2:6724821/1‑93 (MQ=255)
                                                aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc                           >  1:6348370/1‑93 (MQ=255)
                                                 gCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcg                          >  1:240521/1‑93 (MQ=255)
                                                   gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        >  2:6604417/1‑93 (MQ=255)
                                                   gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        <  1:2700355/93‑1 (MQ=255)
                                                    aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgt                       >  1:3774073/1‑93 (MQ=255)
                                                     aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      <  1:8001672/93‑1 (MQ=255)
                                                     aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      >  1:8356279/1‑93 (MQ=255)
                                                     aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      <  1:7834186/93‑1 (MQ=255)
                                                      cTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGAAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttg                     <  2:4131793/93‑1 (MQ=255)
                                                         cTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgt                  >  1:4755887/1‑93 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                >  1:6527696/1‑93 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                >  1:3398780/1‑93 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                <  2:7346455/93‑1 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                <  2:3838354/93‑1 (MQ=255)
                                                             tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              <  1:1296177/93‑1 (MQ=255)
                                                             tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              <  1:1659730/93‑1 (MQ=255)
                                                              ggTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACt             >  2:3366157/1‑93 (MQ=255)
                                                               gtgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTg            >  1:818873/1‑93 (MQ=255)
                                                                tgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAACGGAGGCGGGCTGGATTGGCGCGTTGCGTTCACCgg           >  2:6849508/1‑93 (MQ=255)
                                                                  ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCg         >  1:1559115/1‑93 (MQ=255)
                                                                  ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCg         >  2:2738608/1‑93 (MQ=255)
                                                                   tcatcaTCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGa        <  2:6123664/93‑1 (MQ=255)
                                                                    catcatCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGAc       >  2:1991762/1‑93 (MQ=255)
                                                                    catcatCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGAc       >  2:1172932/1‑93 (MQ=255)
                                                                      tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGACTCGCGCGTTGCGTTCACTGGCGACaa     >  2:5616432/1‑93 (MQ=255)
                                                                       catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    <  2:7931271/93‑1 (MQ=255)
                                                                       catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    <  1:1996794/93‑1 (MQ=255)
                                                                        atcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACg   >  2:7160596/1‑93 (MQ=255)
                                                                         tcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACgg  <  1:5363406/93‑1 (MQ=255)
                                                                                     ||                                                                                 
AAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGG  >  NC_000913/4019775‑4019942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: