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% 834.3 / 3.6 265coding (1136/1386 nt)hofBHofB
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (141/120);  total (143/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 = 1.00e+00

Only 100 of 293 total aligned reads displayed.

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051
                                                                                              |                                                                                         
gACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGt                                                                                            <  1:2761612/93‑1 (MQ=255)
gACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGt                                                                                            >  2:4773175/1‑93 (MQ=255)
 aCGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTg                                                                                          <  2:4416760/93‑1 (MQ=255)
   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTggg                                                                                       <  2:1860733/93‑1 (MQ=255)
   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTggg                                                                                       >  1:2047272/1‑93 (MQ=255)
   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTggg                                                                                       >  2:2024418/1‑93 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggg                                                                                      >  1:2885833/1‑93 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggg                                                                                      >  2:2652182/1‑93 (MQ=255)
    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCACAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCATTgggg                                                                                      >  2:5596287/1‑93 (MQ=255)
       gACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggca                                                                                   <  1:141564/93‑4 (MQ=255)
       gACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggca                                                                                   >  1:4390002/1‑90 (MQ=255)
       gACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTgggggca                                                                                   <  2:3141041/93‑4 (MQ=255)
          cGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcg                                                                                >  1:5941239/1‑93 (MQ=255)
          cGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcg                                                                                <  1:6229611/93‑1 (MQ=255)
           ggCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcgg                                                                               <  1:1779513/93‑1 (MQ=255)
            gCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcggc                                                                              >  1:3770590/1‑93 (MQ=255)
             cGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAgcggcg                                                                             <  1:2869736/93‑1 (MQ=255)
              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGa                                                                            >  1:5864432/1‑93 (MQ=255)
               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGAt                                                                           <  2:6142991/93‑1 (MQ=255)
               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGAt                                                                           <  2:6862262/93‑1 (MQ=255)
                tATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATg                                                                          >  1:5540179/1‑93 (MQ=255)
                 aTGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAACCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATgg                                                                         <  1:4876324/93‑1 (MQ=255)
                  tGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGc                                                                        >  1:2894398/1‑93 (MQ=255)
                   gggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGcc                                                                       >  2:2091428/1‑93 (MQ=255)
                   gggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGcc                                                                       >  2:6227208/1‑93 (MQ=255)
                   gggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGcc                                                                       >  1:7276203/1‑93 (MQ=255)
                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCa                                                                      >  1:5531591/1‑93 (MQ=255)
                     gCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCAt                                                                     >  2:7207188/1‑93 (MQ=255)
                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATAc                                                                   >  1:407865/1‑93 (MQ=255)
                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATAc                                                                   >  2:3724198/1‑93 (MQ=255)
                        gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACACCCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACa                                                                  <  2:839933/93‑1 (MQ=255)
                        gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACa                                                                  <  1:2627640/93‑1 (MQ=255)
                         aaCTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACAt                                                                 <  1:1211483/93‑1 (MQ=255)
                          aCTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGACAGCGGCGATGGCCATACAtt                                                                <  2:7782953/93‑1 (MQ=255)
                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGt                                                              >  2:7835523/1‑93 (MQ=255)
                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGt                                                              >  2:1678816/1‑93 (MQ=255)
                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGt                                                              <  2:3403999/93‑1 (MQ=255)
                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGt                                                              <  2:4220853/93‑1 (MQ=255)
                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTg                                                           <  1:6544269/93‑1 (MQ=255)
                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTg                                                           >  2:161348/1‑93 (MQ=255)
                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTg                                                           <  2:3858182/93‑1 (MQ=255)
                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                          >  2:7694876/1‑93 (MQ=255)
                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                          <  1:1551043/93‑1 (MQ=255)
                                 aTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGa                                                         >  1:3330987/1‑93 (MQ=255)
                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGaa                                                        >  1:3581270/1‑93 (MQ=255)
                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGaa                                                        <  2:465554/93‑1 (MQ=255)
                                   aaGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATCCATTGTCTGGagg                                                       >  1:1910486/1‑91 (MQ=255)
                                    aGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAAtg                                                      <  2:6064662/93‑1 (MQ=255)
                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAAtgg                                                     >  2:5743890/1‑92 (MQ=255)
                                      gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAAtgtg                                                    <  1:6117060/93‑1 (MQ=255)
                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGa                                                  <  1:7160248/93‑1 (MQ=255)
                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  1:6805882/93‑1 (MQ=255)
                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 >  2:2263976/1‑93 (MQ=255)
                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                >  2:5971826/1‑93 (MQ=255)
                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                             <  2:6465183/93‑1 (MQ=255)
                                              aCCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCt                                            >  2:5837294/1‑93 (MQ=255)
                                              aCCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCt                                            >  2:6169449/1‑93 (MQ=255)
                                              aCCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCt                                            >  1:5235825/1‑93 (MQ=255)
                                               ccATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTc                                           >  1:1106144/1‑93 (MQ=255)
                                               ccATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTc                                           <  2:6323001/93‑1 (MQ=255)
                                                cATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 >  1:4574862/1‑86 (MQ=255)
                                                 aTAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTccc                                         >  2:2755344/1‑93 (MQ=255)
                                                 aTAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTccc                                         <  2:4624390/93‑1 (MQ=255)
                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        >  2:2442558/1‑93 (MQ=255)
                                                     aaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtg                                     <  2:1795885/93‑1 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  1:4538087/93‑1 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  1:4798262/93‑1 (MQ=255)
                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  2:1238889/93‑1 (MQ=255)
                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   <  2:1388723/93‑1 (MQ=255)
                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   >  2:7435782/1‑93 (MQ=255)
                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   >  2:1512054/1‑93 (MQ=255)
                                                         cGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgctgc                                 <  1:1806200/93‑1 (MQ=255)
                                                         cGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgctgc                                 <  1:4236241/93‑1 (MQ=255)
                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                >  1:5155218/1‑93 (MQ=255)
                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                <  2:3498277/93‑1 (MQ=255)
                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                >  1:7610383/1‑93 (MQ=255)
                                                            ggTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCgg                              <  1:5672711/93‑1 (MQ=255)
                                                              tAGCAATGTACACAACCGGGTGCCTGCCAGCGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATATGGATGGGCTCCCCTTGCTGCCGGCg                            >  1:937007/1‑93 (MQ=255)
                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAc                          >  1:637912/1‑93 (MQ=255)
                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACa                         >  1:2101184/1‑93 (MQ=255)
                                                                 cAATGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACa                         >  2:2864953/1‑93 (MQ=255)
                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  1:6479768/93‑1 (MQ=255)
                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  1:6776889/93‑1 (MQ=255)
                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  2:4713534/93‑1 (MQ=255)
                                                                   aTGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                       >  1:5079092/1‑93 (MQ=255)
                                                                    tGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      <  2:5086588/93‑1 (MQ=255)
                                                                    tGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      >  2:7594768/1‑93 (MQ=255)
                                                                    tGTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      >  2:6568005/1‑93 (MQ=255)
                                                                     gTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgt                     <  1:7433404/93‑1 (MQ=255)
                                                                     gTACACAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgt                     <  1:3054408/93‑1 (MQ=255)
                                                                      tACACAGCCGGGTGCCTGCCAGTGGAGTCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgtg                    <  1:6960108/93‑1 (MQ=255)
                                                                       acacAACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   <  1:5420775/93‑1 (MQ=255)
                                                                        cacaACCGGGTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTggg                  >  1:1680373/1‑93 (MQ=255)
                                                                         acaACCGGGTGCCTGCGAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGc                 >  2:5218746/1‑93 (MQ=255)
                                                                               gggTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           >  2:5017819/1‑93 (MQ=255)
                                                                                ggTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGttt          <  1:2262562/93‑1 (MQ=255)
                                                                                 gTGCCTGCCAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg         <  1:2531701/93‑1 (MQ=255)
                                                                                    cctgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGt      >  2:1011701/1‑93 (MQ=255)
                                                                                     ctgccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTa     >  1:7913916/1‑93 (MQ=255)
                                                                                       gccAGTGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCTGGCGACAATGTGGGCAAAGTTTGCGTAcc   <  2:6726709/93‑1 (MQ=255)
                                                                                              |                                                                                         
GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

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