Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,166,935 2 bp→CG coding (312‑313/483 nt) b3023 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,9350GC99.1% 783.3 / NA 221L105V (CTT→GTT) b3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (97/121);  minor base A (2/0);  total (99/122)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,166,9360CG99.6% 788.1 / ‑4.8 221Q104H (CAG→CACb3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (99/121);  total (99/122)
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.66e-01

Only 100 of 227 total aligned reads displayed.

TGCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGC  >  NC_000913/3166849‑3167021
                                                                                      ||                                                                                     
 gCTTTACGCATTTTCTCGCCTCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                 <  1:2042002/93‑1 (MQ=255)
 gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                 >  2:3827515/1‑93 (MQ=255)
 gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                 >  2:84911/1‑93 (MQ=255)
  cTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCggg                                                                                <  2:5726422/93‑1 (MQ=255)
   tttACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                                               >  2:1749153/1‑93 (MQ=255)
   tttACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                                               >  2:2534165/1‑93 (MQ=255)
   tttACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                                               >  2:221279/1‑93 (MQ=255)
     tACGCTTTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCa                                                                             <  2:4163586/93‑1 (MQ=255)
      aCGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGcccc                                                                            >  2:1664854/1‑91 (MQ=255)
        gCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                                          <  2:4057975/93‑1 (MQ=255)
          aTTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                                        >  1:5335372/1‑93 (MQ=255)
          aTTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                                        >  2:3021394/1‑93 (MQ=255)
           ttttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGGCACGGCAt                                                                       >  1:2622222/1‑93 (MQ=255)
            tttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATa                                                                      >  2:4557723/1‑93 (MQ=255)
            tttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGGGGCGGGCCACGGCATa                                                                      >  2:1463444/1‑93 (MQ=255)
              tctcGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAAc                                                                    <  2:4209910/93‑1 (MQ=255)
                 cgccgcTTGCAGGGAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                                 <  1:1566973/93‑1 (MQ=255)
                 cgccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                                 >  2:3394132/1‑93 (MQ=255)
                    cgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTTGCTAACGTGGCGGGCCACGGCATAACGTCCAc                                                              >  1:5643537/1‑93 (MQ=255)
                    cgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                                              <  2:264512/93‑1 (MQ=255)
                     gcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGGCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcc                                                             >  2:4285045/1‑93 (MQ=255)
                      cTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCg                                                            <  2:298262/93‑1 (MQ=255)
                        tGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACTGGt                                                          <  2:542350/93‑1 (MQ=255)
                         gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTa                                                         <  2:5634740/93‑1 (MQ=255)
                          caggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTaa                                                        <  1:1058811/93‑1 (MQ=255)
                                  cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                                >  2:3441428/1‑93 (MQ=255)
                                  cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                                <  2:3104223/93‑1 (MQ=255)
                                  cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                                <  2:3905702/93‑1 (MQ=255)
                                   cgTGGCGAATGATGCCCCATGCCGTGTGTGAAATCTCGTCCAGCTCGCCACCGTGGCGGGCCACGGCCTACCGTCCACCGGTAAGTTCACCAt                                               <  2:5495227/92‑1 (MQ=255)
                                   aGTGGCGAATGATGCGCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                                               >  2:5994658/1‑93 (MQ=255)
                                   aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                                               >  1:5508038/1‑93 (MQ=255)
                                   aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                                               <  2:646489/93‑1 (MQ=255)
                                   aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                                               >  2:5234897/1‑93 (MQ=255)
                                    gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCCCGGCATAACGTCCACCCGTAAGTTCCCCAtt                                              >  1:83908/1‑93 (MQ=255)
                                    gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAtt                                              >  1:2059763/1‑93 (MQ=255)
                                     tGGCGAATGACGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                                             >  1:2335820/1‑93 (MQ=255)
                                         gAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGATACGGCATAACGTCCACCGGTAAGTTCACCATTGCTca                                         >  2:4396012/1‑93 (MQ=255)
                                            tgatgCCCCATTCCTTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCCCCGGTAAGTTCACCATTGCTCacac                                      <  1:352172/93‑1 (MQ=255)
                                             gatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACAcc                                     >  2:6082912/1‑93 (MQ=255)
                                              atgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCa                                    >  2:4657229/1‑93 (MQ=255)
                                               tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                                   <  2:1545772/93‑1 (MQ=255)
                                               tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                                   <  1:3547960/93‑1 (MQ=255)
                                               tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                                   <  1:2210589/93‑1 (MQ=255)
                                               tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                                   <  2:1006223/93‑1 (MQ=255)
                                                 ccccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATaa                                 <  2:4698567/93‑1 (MQ=255)
                                                  gcaaTACCGTGTGAGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                                <  2:4871583/90‑1 (MQ=255)
                                                  cccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                                <  2:2950566/93‑1 (MQ=255)
                                                  cccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                                >  2:2736967/1‑93 (MQ=255)
                                                   ccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACg                               >  2:1364154/1‑93 (MQ=255)
                                                   ccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACg                               >  1:543398/1‑93 (MQ=255)
                                                    cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa                              >  1:1363177/1‑93 (MQ=255)
                                                    cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa                              >  1:4659828/1‑93 (MQ=255)
                                                    cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa                              <  2:2566426/93‑1 (MQ=255)
                                                    cATACCGTGTGTGAAAGATCGTCCAGCTCGCCAACGTGGCGGGACACGGCATAACGTCCAGCGGTAAGTTCACCATTGCTCACACCATAACGa                              <  2:2492350/93‑1 (MQ=255)
                                                     aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAt                             >  1:1193796/1‑93 (MQ=255)
                                                     aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACAGGTAAGTTCACCATTGCTCACACCATAACGAt                             >  1:2231041/1‑93 (MQ=255)
                                                      tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt                            <  2:5875756/93‑1 (MQ=255)
                                                       aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa                           <  1:2836812/93‑1 (MQ=255)
                                                       aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa                           <  1:4480027/93‑1 (MQ=255)
                                                        ccGTGGTTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTAt                          <  2:2086772/93‑1 (MQ=255)
                                                          gtgtgtGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCg                        >  1:3026452/1‑93 (MQ=255)
                                                          gtgtgtGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCg                        <  1:4693578/93‑1 (MQ=255)
                                                           tgtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg                       <  1:4015777/93‑1 (MQ=255)
                                                           tgtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACAGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg                       <  2:2330606/93‑1 (MQ=255)
                                                            gtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCggg                      <  2:2364135/93‑1 (MQ=255)
                                                            gtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCggg                      >  2:817753/1‑93 (MQ=255)
                                                              gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                    >  2:3680413/1‑93 (MQ=255)
                                                              gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                    <  2:936725/93‑1 (MQ=255)
                                                               tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt                   <  2:2006448/93‑1 (MQ=255)
                                                               tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt                   <  2:1855850/93‑1 (MQ=255)
                                                                gAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATc                  >  1:3458568/1‑93 (MQ=255)
                                                                gAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATc                  <  1:2221374/93‑1 (MQ=255)
                                                                gAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATc                  >  1:3674606/1‑93 (MQ=255)
                                                                 aaaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCg                 <  1:2104777/93‑1 (MQ=255)
                                                                 aaaTATCGTCCAGCTCGCCAACGTGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCgg                <  1:76068/93‑1 (MQ=255)
                                                                   atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGt               >  2:1071575/1‑93 (MQ=255)
                                                                    tatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGtt              <  2:1221098/93‑1 (MQ=255)
                                                                    tatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGtt              <  1:3728140/93‑1 (MQ=255)
                                                                    tatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGtt              <  1:4247416/93‑1 (MQ=255)
                                                                     atCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTc             >  1:1087592/1‑93 (MQ=255)
                                                                     atCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTc             >  1:4052199/1‑93 (MQ=255)
                                                                     atCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTc             <  1:5930571/93‑1 (MQ=255)
                                                                      tCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCg            >  1:3116331/1‑93 (MQ=255)
                                                                       tgTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGc           <  1:5677879/92‑1 (MQ=255)
                                                                       cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGc           <  2:5465344/93‑1 (MQ=255)
                                                                       cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGc           >  1:4870063/1‑93 (MQ=255)
                                                                        ctCCAGCGCCCCAACGGGGCGGGCCACGGCATAACGTCCACCGGTAAATTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCt          <  1:4779488/92‑1 (MQ=255)
                                                                         tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTa         >  1:4942198/1‑93 (MQ=255)
                                                                          ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa        <  1:336593/93‑1 (MQ=255)
                                                                          ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa        <  2:2847709/93‑1 (MQ=255)
                                                                          ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa        <  1:6117109/93‑1 (MQ=255)
                                                                           cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAAc       >  1:4115714/1‑93 (MQ=255)
                                                                           cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAAc       <  1:721706/93‑1 (MQ=255)
                                                                           cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAAc       <  1:2574603/93‑1 (MQ=255)
                                                                            aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACg      <  1:3884127/93‑1 (MQ=255)
                                                                             gCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGc     >  2:2794574/1‑93 (MQ=255)
                                                                             gCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGc     >  1:311355/1‑93 (MQ=255)
                                                                             gCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGc     <  1:319358/93‑1 (MQ=255)
                                                                              cTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCt    <  2:4776631/93‑1 (MQ=255)
                                                                               tCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTg   >  2:5375205/1‑93 (MQ=255)
                                                                                      ||                                                                                     
TGCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGC  >  NC_000913/3166849‑3167021

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