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.7% 1677.6 / ‑3.0 373intergenic (+59/‑25)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (184/187);  total (184/189)
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.74e-01
*NC_0009134,019,8610C.99.5% 1683.5 / ‑1.0 373intergenic (+60/‑24)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (184/187);  minor base G (0/2);  total (184/189)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 380 total aligned reads displayed.

AAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGG  >  NC_000913/4019775‑4019942
                                                                                     ||                                                                                 
aaGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTgggg                                                                           >  1:4355588/1‑93 (MQ=255)
 aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTTATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          >  1:2072372/1‑93 (MQ=255)
 aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          >  1:7480689/1‑93 (MQ=255)
 aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          <  2:5187567/93‑1 (MQ=255)
 aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          <  1:2236875/93‑1 (MQ=255)
    gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       <  2:2262267/93‑1 (MQ=255)
    gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       >  2:6742447/1‑93 (MQ=255)
     ccACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCa                                                                      >  1:5582519/1‑93 (MQ=255)
     ccACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCa                                                                      >  1:1896695/1‑93 (MQ=255)
      cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     <  1:2401275/93‑1 (MQ=255)
      cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     >  2:1425235/1‑93 (MQ=255)
      cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     <  1:9509634/93‑1 (MQ=255)
       aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    >  1:2812270/1‑93 (MQ=255)
       aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    >  1:7934246/1‑93 (MQ=255)
       aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    <  1:9587753/93‑1 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  2:3584980/1‑93 (MQ=255)
         gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  2:8699202/1‑93 (MQ=255)
          aTAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                                                 >  1:7041429/1‑93 (MQ=255)
           tAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                                                >  2:2670353/1‑93 (MQ=255)
            aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                                               <  1:9883092/93‑1 (MQ=255)
            aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                                               >  1:7793980/1‑93 (MQ=255)
            aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                                               >  1:6249348/1‑93 (MQ=255)
            aaaGAGCAGGTGTAATCCGTGCTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCGTCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                                               >  1:8168493/1‑93 (MQ=255)
             aaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                                              <  2:8619452/93‑1 (MQ=255)
              agagCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGcc                                                             <  1:3940403/93‑1 (MQ=255)
               gagCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                                            <  1:1779874/93‑1 (MQ=255)
               aagCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                                            <  1:6338241/92‑1 (MQ=255)
                 gCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                                          >  2:9877941/1‑93 (MQ=255)
                 gCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                                          >  1:6027048/1‑93 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  1:5657519/1‑93 (MQ=255)
                   aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  1:3437408/1‑93 (MQ=255)
                    ggTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                                       <  2:8123311/93‑1 (MQ=255)
                    ggTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                                       <  1:3669317/93‑1 (MQ=255)
                    ggTCTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                                       <  1:738204/93‑1 (MQ=255)
                     gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                                                      >  2:6618922/1‑93 (MQ=255)
                      tgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                                     <  2:3866228/93‑1 (MQ=255)
                        tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                                   >  2:1330254/1‑93 (MQ=255)
                          aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaa                                                 <  2:3480987/93‑1 (MQ=255)
                           tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  2:9413041/1‑93 (MQ=255)
                            ccGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaaa                                               <  2:4035988/93‑1 (MQ=255)
                             cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              >  2:1029831/1‑93 (MQ=255)
                             cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              <  2:5881373/93‑1 (MQ=255)
                              gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGTCGGTAAAAACg                                             >  1:1667836/1‑93 (MQ=255)
                              gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                                             >  1:4785104/1‑93 (MQ=255)
                                gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGt                                           >  2:3431137/1‑93 (MQ=255)
                                 tttGATATCGTTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          <  2:9175843/93‑1 (MQ=255)
                                 tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          <  2:1579561/93‑1 (MQ=255)
                                 tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          >  2:6830856/1‑93 (MQ=255)
                                   tGATATCGGTTTTAGCGAACTGCTATTGTTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                                        >  2:6083534/1‑93 (MQ=255)
                                   tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                                        <  1:3032664/93‑1 (MQ=255)
                                    gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                                       <  1:8387936/93‑1 (MQ=255)
                                    gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                                       >  1:46011/1‑93 (MQ=255)
                                     atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCgg                                      <  1:6498747/93‑1 (MQ=255)
                                     atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCgg                                      <  1:5025507/93‑1 (MQ=255)
                                      tatCGGTTTTAGCGAACTGCTAGTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCggg                                     >  2:6409402/1‑93 (MQ=255)
                                       atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGc                                    <  1:1513784/93‑1 (MQ=255)
                                        tCGGTTTTAGCGAACTGCGATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt                                   <  1:5219578/93‑1 (MQ=255)
                                             tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt                              <  2:4551399/93‑1 (MQ=255)
                                             tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt                              >  2:7160156/1‑93 (MQ=255)
                                               tAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcg                            <  1:3849428/93‑1 (MQ=255)
                                                tgccaaCTGCTATTGGTTTTCCTCATCCGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc                           <  1:3184453/89‑1 (MQ=255)
                                                  cGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgc                         <  2:8370992/93‑1 (MQ=255)
                                                   gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        <  2:6242780/93‑1 (MQ=255)
                                                   gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        >  1:9971365/1‑93 (MQ=255)
                                                    aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgt                       >  2:7632782/1‑93 (MQ=255)
                                                    aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGCAGCGGGCTGGATTCGCgcgt                       >  2:2040989/1‑93 (MQ=255)
                                                     aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      >  1:3494947/1‑93 (MQ=255)
                                                     aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      <  2:7725175/93‑1 (MQ=255)
                                                      cTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttg                     <  2:9665536/93‑1 (MQ=255)
                                                      cTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttg                     >  2:5782839/1‑93 (MQ=255)
                                                       tGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgc                    <  2:178319/93‑1 (MQ=255)
                                                          tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgtt                 >  1:7108329/1‑93 (MQ=255)
                                                          tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgtt                 >  1:86781/1‑93 (MQ=255)
                                                          tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgtt                 >  1:9039449/1‑93 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgttg                >  2:3310727/1‑92 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                >  2:3751276/1‑93 (MQ=255)
                                                           aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                <  2:3112538/93‑1 (MQ=255)
                                                            ttGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCa               <  1:2507683/93‑1 (MQ=255)
                                                            ttGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCa               >  2:8243116/1‑93 (MQ=255)
                                                             tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              <  1:5452918/93‑1 (MQ=255)
                                                             tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              <  2:7793979/93‑1 (MQ=255)
                                                             tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              >  1:5906299/1‑93 (MQ=255)
                                                               gtgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTg            >  2:2910987/1‑93 (MQ=255)
                                                                tgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTgg           >  1:1035261/1‑93 (MQ=255)
                                                                tgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTgg           <  2:9318123/93‑1 (MQ=255)
                                                                 gtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGc          >  2:5511559/1‑93 (MQ=255)
                                                                  ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTgggg         >  1:9710669/1‑91 (MQ=255)
                                                                  ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGca         <  1:4198378/93‑2 (MQ=255)
                                                                  ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCg         >  2:4854297/1‑93 (MQ=255)
                                                                     atcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACa      >  2:7074902/1‑93 (MQ=255)
                                                                     atcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACa      <  1:1364871/93‑1 (MQ=255)
                                                                      tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTTGATTCGCGCGTTGCGTTCACTGGCGACaa     <  2:9064651/93‑1 (MQ=255)
                                                                      tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     >  1:9300520/1‑93 (MQ=255)
                                                                      tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     <  1:1009997/93‑1 (MQ=255)
                                                                       catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    <  2:4775775/93‑1 (MQ=255)
                                                                       catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    >  2:4368267/1‑93 (MQ=255)
                                                                       catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    >  1:8800739/1‑93 (MQ=255)
                                                                        atcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGCAAAAACGGTAGCGGGCTGGATTGGCGCGTTGCGTTCACTGGCGAAAACg   >  1:3311110/1‑93 (MQ=255)
                                                                         tcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACgg  <  2:5307926/93‑1 (MQ=255)
                                                                         tcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACgg  <  2:8283180/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: