Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,627,142 +G intergenic (‑24/+23) yhiH ← / ← yhiI putative ATP‑binding component of a transport system, fragment 1/putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,627,1401.G99.6% 745.9 / ‑4.3 233intergenic (‑22/+25)yhiH/yhiIputative ATP‑binding component of a transport system, fragment 1/putative membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (134/98);  total (135/98)
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.99e-01

Only 100 of 277 total aligned reads displayed.

ACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCA  >  NC_000913/3627048‑3627227
                                                                                             |                                                                                       
aCCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                          <  2:3776689/93‑1 (MQ=255)
 ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCg                                                                                        >  2:6083111/1‑93 (MQ=255)
 ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCg                                                                                        <  2:1898925/93‑1 (MQ=255)
 ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCg                                                                                        <  2:2444146/93‑1 (MQ=255)
 ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCg                                                                                        <  1:3467865/93‑1 (MQ=255)
 ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTCCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCg                                                                                        <  2:2714517/93‑1 (MQ=255)
  cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCgg                                                                                       >  1:987636/1‑93 (MQ=255)
  cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGATGCGCGAAAGGCgg                                                                                       >  2:3186105/1‑93 (MQ=255)
   aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCggg                                                                                      <  1:2214603/93‑2 (MQ=255)
   aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCggg                                                                                      >  1:4342352/1‑92 (MQ=255)
    tACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCggga                                                                                     >  1:233825/1‑91 (MQ=255)
      cAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCggggcg                                                                                   >  1:5040681/1‑89 (MQ=255)
      cAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCgggacg                                                                                   <  1:5576349/93‑5 (MQ=255)
      cAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCgggacg                                                                                   <  2:4466478/93‑5 (MQ=255)
        gCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACggg                                                                                  >  2:549831/1‑93 (MQ=255)
        gCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACggg                                                                                  <  2:2949362/93‑1 (MQ=255)
        gCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACggg                                                                                  <  1:3129365/93‑1 (MQ=255)
                  atatCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                        >  1:773058/1‑93 (MQ=255)
                   tatCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                       <  1:4914715/93‑1 (MQ=255)
                    atCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAg                                                                      >  1:2788960/1‑93 (MQ=255)
                     tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                     <  1:4064684/93‑1 (MQ=255)
                     tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                     >  1:2607756/1‑93 (MQ=255)
                     tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                     <  1:4498134/93‑1 (MQ=255)
                      cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                    <  1:6044456/93‑1 (MQ=255)
                      cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                    <  1:3211936/93‑1 (MQ=255)
                        agagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGc                                                                  >  2:3892948/1‑93 (MQ=255)
                             gATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCAt                                                             <  2:6252708/93‑1 (MQ=255)
                              atatTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGGACCAGGTCCAGATGCGTCAtt                                                            >  2:1310811/1‑93 (MQ=255)
                                 ttGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCg                                                         >  1:5768559/1‑93 (MQ=255)
                                 ttGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCg                                                         >  2:193529/1‑93 (MQ=255)
                                    ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCa                                                      >  1:4979504/1‑93 (MQ=255)
                                     tCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                     <  1:2525138/93‑1 (MQ=255)
                                     tCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                     >  2:4720263/1‑93 (MQ=255)
                                     tCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                     >  2:692706/1‑93 (MQ=255)
                                      cAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACCGGAACCAGTTCCCGATGCGTCATTGCCGCCAc                                                    >  2:2546020/1‑93 (MQ=255)
                                       aGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                   <  1:4622190/93‑1 (MQ=255)
                                       aGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                   >  2:5085269/1‑93 (MQ=255)
                                        gcgcAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCt                                                  >  1:1790195/1‑93 (MQ=255)
                                          gcAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTca                                                >  1:2410825/1‑93 (MQ=255)
                                             aCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                             <  2:4902634/93‑1 (MQ=255)
                                             aCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                             <  2:4548958/93‑1 (MQ=255)
                                              cGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcaccac                                            >  1:6071581/1‑93 (MQ=255)
                                               ggTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACg                                           >  2:4255137/1‑93 (MQ=255)
                                                gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGa                                          <  1:2911538/93‑1 (MQ=255)
                                                gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGa                                          <  2:4017329/93‑1 (MQ=255)
                                                 tttttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGTCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                         >  1:4449638/1‑93 (MQ=255)
                                                 tttttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                         <  2:568522/93‑1 (MQ=255)
                                                  ttttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                        >  1:3593376/1‑93 (MQ=255)
                                                  ttttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                        <  2:1247780/93‑1 (MQ=255)
                                                   tttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                       <  1:6180266/93‑1 (MQ=255)
                                                   tttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                       >  2:3927336/1‑93 (MQ=255)
                                                    ttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                      >  1:5698503/1‑93 (MQ=255)
                                                     tCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcg                                     >  1:1311231/1‑93 (MQ=255)
                                                      ccATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                    >  2:1474172/1‑93 (MQ=255)
                                                      ccATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                    >  1:4000036/1‑93 (MQ=255)
                                                      ccATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCAACACGAGgtcgt                                    >  1:3296101/1‑93 (MQ=255)
                                                       cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtc                                   >  2:2311443/1‑93 (MQ=255)
                                                        aTAATGCTGTCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCa                                  <  2:5841838/93‑1 (MQ=255)
                                                        aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCCCCACGAGGTCGTCa                                  <  1:5268781/93‑1 (MQ=255)
                                                        aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtcc                                  >  2:1600396/1‑92 (MQ=255)
                                                        aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCa                                  <  2:4098155/93‑1 (MQ=255)
                                                          aaTGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAgg                                >  1:856873/1‑93 (MQ=255)
                                                           aTGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCCACCTCACCACGAGGTCGTCCGGc                               >  2:3601859/1‑93 (MQ=255)
                                                           aTGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                               >  2:3105425/1‑93 (MQ=255)
                                                            tGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGcc                              >  2:1084536/1‑93 (MQ=255)
                                                            tGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGcc                              >  2:3373060/1‑93 (MQ=255)
                                                             gCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCa                             >  1:1446469/1‑93 (MQ=255)
                                                             gCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCa                             >  1:4215803/1‑93 (MQ=255)
                                                              cTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                            >  2:2483410/1‑93 (MQ=255)
                                                              cTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                            <  1:4597381/93‑1 (MQ=255)
                                                               tGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACg                           >  2:1197771/1‑93 (MQ=255)
                                                                ggCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGGACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGGCGTCAGGCCACgg                          >  1:3243067/1‑93 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGTGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                       <  1:1958991/93‑1 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGacg                       >  2:6071009/1‑91 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                       >  2:479020/1‑93 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                       >  1:2494280/1‑93 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                       <  2:1455185/93‑1 (MQ=255)
                                                                   tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                       >  2:5768053/1‑93 (MQ=255)
                                                                    cACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGt                      >  2:5126804/1‑93 (MQ=255)
                                                                     aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGtt                     >  2:2780032/1‑93 (MQ=255)
                                                                     aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGtt                     >  1:5108768/1‑93 (MQ=255)
                                                                      cgcccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                    <  2:95092/93‑1 (MQ=255)
                                                                          cgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttca                >  2:1765422/1‑93 (MQ=255)
                                                                           gccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcat               >  1:1242273/1‑93 (MQ=255)
                                                                           gccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcat               >  2:5380004/1‑93 (MQ=255)
                                                                            ccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcatt              >  1:1142111/1‑93 (MQ=255)
                                                                            ccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcatt              <  2:2800193/93‑1 (MQ=255)
                                                                            ccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcatt              >  1:5977336/1‑93 (MQ=255)
                                                                             cAGTTGCGCGACAGGCGGGCCGGGAACCAGTTCCCGATGCGTCATTGGGGCAACCTCACCACGAGGACGTCAGGGCACGGGAGATattcattc             >  2:524946/1‑93 (MQ=255)
                                                                             cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc             >  2:3738102/1‑93 (MQ=255)
                                                                             cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc             >  1:2805195/1‑93 (MQ=255)
                                                                             cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc             >  2:2605625/1‑93 (MQ=255)
                                                                              aGTTGCGCGACAGTCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattca            >  1:3384066/1‑93 (MQ=255)
                                                                              aGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattca            <  1:3720829/93‑1 (MQ=255)
                                                                              aGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattca            <  2:4170693/93‑1 (MQ=255)
                                                                                 tGCGCGAGAGGCGGGACGGGAACCAGTGCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAccc         <  2:561944/93‑1 (MQ=255)
                                                                                 tGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAccc         <  2:2073291/93‑1 (MQ=255)
                                                                                      gACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc    <  1:5413083/93‑1 (MQ=255)
                                                                                      gACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc    >  2:50134/1‑93 (MQ=255)
                                                                                       aCAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc   >  1:2232115/1‑93 (MQ=255)
                                                                                             |                                                                                       
ACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCA  >  NC_000913/3627048‑3627227

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