Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 115,968 +G coding (1132/1386 nt) hofB ← HofB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913115,9641.G98.5% 653.8 / 1.3 204coding (1136/1386 nt)hofBHofB
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (92/109);  total (93/111)
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.96e-01

Only 100 of 230 total aligned reads displayed.

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051
                                                                                             |                                                                                        
gACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGt                                                                                           >  1:45466/1‑93 (MQ=255)
 aCGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTg                                                                                         <  2:4040789/93‑1 (MQ=255)
  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgg                                                                                        <  1:1236884/93‑1 (MQ=255)
   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTggg                                                                                       >  2:3540336/1‑93 (MQ=255)
   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTggg                                                                                       >  1:2773527/1‑93 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggg                                                                                     <  2:324743/93‑1 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggg                                                                                     >  1:4695334/1‑93 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggg                                                                                     >  1:3112896/1‑93 (MQ=255)
      tGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggc                                                                                   <  2:102061/93‑3 (MQ=255)
       gACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggca                                                                                  >  1:5412338/1‑90 (MQ=255)
        aCCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggcag                                                                                 <  2:2275851/93‑5 (MQ=255)
           ggCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcgg                                                                               >  2:5131252/1‑93 (MQ=255)
           ggCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcgg                                                                               <  2:4105551/93‑1 (MQ=255)
           ggCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcgg                                                                               >  1:63312/1‑93 (MQ=255)
             cGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcggcg                                                                             >  1:1149988/1‑93 (MQ=255)
             cGTTATGGGCAGAACTTCAAATAAGGCCGTAAGACCATAAAAACCGTGGTAGCAATGTACACAACCGAGTGCCTGCCAGTCGGGGCAgcggcg                                                                             >  2:4949935/1‑93 (MQ=255)
              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGa                                                                            <  1:6102459/93‑1 (MQ=255)
               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGAt                                                                           <  2:2231207/93‑1 (MQ=255)
               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGAt                                                                           <  2:2573650/93‑1 (MQ=255)
                  tGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGc                                                                        <  2:5679847/93‑1 (MQ=255)
                  tGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGc                                                                        <  1:3923547/93‑1 (MQ=255)
                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCa                                                                      <  2:5733637/93‑1 (MQ=255)
                     gCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCAt                                                                     >  1:1840112/1‑93 (MQ=255)
                     gCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCAt                                                                     <  1:2740518/93‑1 (MQ=255)
                      cAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATa                                                                    <  1:1186611/93‑1 (MQ=255)
                      cAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATa                                                                    <  1:5750903/93‑1 (MQ=255)
                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATAc                                                                   <  2:4579135/93‑1 (MQ=255)
                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATAc                                                                   >  1:3773721/1‑93 (MQ=255)
                        gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACa                                                                  >  1:4615304/1‑93 (MQ=255)
                        gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACa                                                                  <  1:2865513/93‑1 (MQ=255)
                          aCTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACAtt                                                                >  2:361602/1‑93 (MQ=255)
                           cTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTg                                                               <  1:3850513/93‑1 (MQ=255)
                              cAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCt                                                            >  1:921454/1‑93 (MQ=255)
                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTg                                                           <  2:4458406/93‑1 (MQ=255)
                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCAGCGATGGCCATACATTGTCTg                                                           <  1:4995031/93‑1 (MQ=255)
                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                          >  1:405782/1‑93 (MQ=255)
                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                          >  1:5577004/1‑93 (MQ=255)
                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                          >  1:2223973/1‑93 (MQ=255)
                                 aTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGa                                                         <  1:6058683/93‑1 (MQ=255)
                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAAtgt                                                     <  2:3832258/93‑1 (MQ=255)
                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGa                                                  <  2:1496032/93‑1 (MQ=255)
                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGa                                                  <  2:3408481/93‑1 (MQ=255)
                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGAGGCGATGGCCATACATTGTCTGGAATGTGGa                                                  >  2:5975290/1‑93 (MQ=255)
                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 >  1:2576963/1‑93 (MQ=255)
                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  2:2687158/93‑1 (MQ=255)
                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  1:79289/93‑1 (MQ=255)
                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  1:352483/93‑1 (MQ=255)
                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  2:5020019/93‑1 (MQ=255)
                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  1:5932327/93‑1 (MQ=255)
                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                >  1:5219014/1‑93 (MQ=255)
                                           aCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATgg                                               <  1:4169680/93‑1 (MQ=255)
                                            cGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATggg                                              <  1:1005718/93‑1 (MQ=255)
                                            cGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATggg                                              >  2:1320812/1‑93 (MQ=255)
                                            cGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATATGGATggg                                              <  2:6144829/93‑1 (MQ=255)
                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                             <  2:3397317/93‑1 (MQ=255)
                                              aCCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCt                                            >  1:245610/1‑93 (MQ=255)
                                                 aTAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTccc                                         >  2:1682912/1‑93 (MQ=255)
                                                 aTAAAAACCGTGGTAGCAATGTACACAACCGGGAGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTccc                                         >  1:253913/1‑93 (MQ=255)
                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        <  2:3287393/93‑1 (MQ=255)
                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        <  1:1679421/93‑1 (MQ=255)
                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCt                                       <  2:5443270/93‑1 (MQ=255)
                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCt                                       <  2:826893/93‑1 (MQ=255)
                                                    aaaaCCGTGGTAGCAATGTACACAAACGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      >  2:1254401/1‑93 (MQ=255)
                                                     aaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtg                                     >  1:2521734/1‑93 (MQ=255)
                                                      agccGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  1:5651230/91‑1 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    >  1:4793656/1‑93 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  2:1556711/93‑1 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  1:3723257/93‑1 (MQ=255)
                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   >  2:2829931/1‑93 (MQ=255)
                                                         cGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgctgc                                 >  1:5318013/1‑93 (MQ=255)
                                                              tAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCg                            >  2:4809379/1‑93 (MQ=255)
                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAc                          >  2:761148/1‑93 (MQ=255)
                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAc                          >  1:1607703/1‑93 (MQ=255)
                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAc                          >  2:5918478/1‑93 (MQ=255)
                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  2:5221289/93‑1 (MQ=255)
                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCAATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        >  2:3844499/1‑93 (MQ=255)
                                                                   aTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                       >  1:5375232/1‑93 (MQ=255)
                                                                    tGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      <  2:301763/93‑1 (MQ=255)
                                                                    tGTACACAACCGGGTGCCTGCCAGTGGGCGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      >  2:5516021/1‑93 (MQ=255)
                                                                     gTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgt                     >  2:6099739/1‑93 (MQ=255)
                                                                       acacAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   >  1:5182028/1‑93 (MQ=255)
                                                                       acacAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   <  2:65452/93‑1 (MQ=255)
                                                                       acacAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   <  2:1937095/93‑1 (MQ=255)
                                                                           aaCCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaa               >  1:6179044/1‑93 (MQ=255)
                                                                           aaCCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaa               <  1:1989510/93‑1 (MQ=255)
                                                                            aCCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaaa              >  1:1905867/1‑93 (MQ=255)
                                                                              cGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGt            >  2:4708266/1‑93 (MQ=255)
                                                                               gggTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           <  2:2402857/93‑1 (MQ=255)
                                                                                ggTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGttt          <  1:2310616/93‑1 (MQ=255)
                                                                                 gTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg         <  2:3210626/93‑1 (MQ=255)
                                                                                   gcctgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCg       >  2:175082/1‑93 (MQ=255)
                                                                                   gcctgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCg       <  2:1837066/93‑1 (MQ=255)
                                                                                      tgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    >  2:4177840/1‑93 (MQ=255)
                                                                                      tgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    <  1:1481646/93‑1 (MQ=255)
                                                                                      tgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    >  2:4005521/1‑93 (MQ=255)
                                                                                      tgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    >  2:3876208/1‑93 (MQ=255)
                                                                                      tgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    >  2:141189/1‑93 (MQ=255)
                                                                                       gccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAcc   >  1:3278977/1‑93 (MQ=255)
                                                                                        ccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCa  <  1:2620320/93‑1 (MQ=255)
                                                                                        ccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCa  <  1:2415429/93‑1 (MQ=255)
                                                                                             |                                                                                        
GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

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