Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 368,911 G→A V359V (GTG→GTA)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913368,9110GA99.7% 1065.7 / ‑5.8 315V359V (GTG→GTA) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (154/160);  total (154/161)
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.34e-01

Only 100 of 339 total aligned reads displayed.

AGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCT  >  NC_000913/368821‑369000
                                                                                          |                                                                                         
aGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACt                                                                                         <  1:6777854/93‑1 (MQ=255)
aGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACt                                                                                         >  2:3243582/1‑93 (MQ=255)
 gCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTg                                                                                        >  1:4275627/1‑93 (MQ=255)
   ccGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGc                                                                                      <  2:7750740/93‑1 (MQ=255)
   ccGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGc                                                                                      <  2:8822841/93‑1 (MQ=255)
     gcgATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTc                                                                                    <  1:1034418/93‑1 (MQ=255)
       gATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccg                                                                                  <  2:5725712/93‑1 (MQ=255)
         tGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgcc                                                                                >  2:5087526/1‑93 (MQ=255)
          gcgcTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgccg                                                                               <  2:1293998/93‑1 (MQ=255)
          gcgcTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgccg                                                                               <  1:2334926/93‑1 (MQ=255)
             ctgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGcaa                                                                            >  2:7990951/1‑93 (MQ=255)
              tgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTTCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAAc                                                                           >  1:134098/1‑93 (MQ=255)
              tgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAAc                                                                           <  2:3398891/93‑1 (MQ=255)
               gctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGCCTCCGCCGAAAcc                                                                          >  1:2604829/1‑92 (MQ=255)
               gctCGATACCTATCAACAAGAACGTCGCGATCACGACAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTAATGGGTGCGCCGaaaag                                                                          >  2:3782070/1‑91 (MQ=255)
                 tCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCt                                                                        <  2:2860066/93‑1 (MQ=255)
                 tCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAAAGCt                                                                        >  2:8203338/1‑93 (MQ=255)
                  cGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTg                                                                       >  2:711487/1‑93 (MQ=255)
                    aTACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGc                                                                     >  1:8169411/1‑93 (MQ=255)
                    aTACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGc                                                                     <  2:3434923/93‑1 (MQ=255)
                     tACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCa                                                                    >  2:3066310/1‑93 (MQ=255)
                     tACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCa                                                                    >  1:5379631/1‑93 (MQ=255)
                     tACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGCACTGGCTCCGCCGAAACCCTGGCa                                                                    >  2:3125076/1‑93 (MQ=255)
                      aCCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAg                                                                   <  1:8060665/93‑1 (MQ=255)
                       ccTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAgg                                                                  >  1:7272089/1‑93 (MQ=255)
                        cTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAggg                                                                 <  2:5325106/93‑1 (MQ=255)
                         tATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGt                                                                >  1:688400/1‑93 (MQ=255)
                          aTCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGta                                                               <  1:3481511/93‑1 (MQ=255)
                          aTCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGta                                                               <  2:2593512/93‑1 (MQ=255)
                          aTCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGta                                                               <  2:651355/93‑1 (MQ=255)
                           tCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtac                                                              <  1:4095280/93‑1 (MQ=255)
                             aacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgt                                                            <  1:5794123/93‑1 (MQ=255)
                             aacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgt                                                            <  1:3006446/93‑1 (MQ=255)
                              acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGCCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGgacgtt                                                           >  1:4485723/1‑93 (MQ=255)
                              acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgtt                                                           >  1:5842762/1‑93 (MQ=255)
                              acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgtt                                                           <  1:5786101/93‑1 (MQ=255)
                               caaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGgacgtta                                                          >  1:4770195/1‑93 (MQ=255)
                                aaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgttac                                                         >  1:7833912/1‑93 (MQ=255)
                                aaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGATCCGCCGAAACGCTGGCAGGGgacgttac                                                         >  1:7914841/1‑93 (MQ=255)
                                 aGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGgacattacg                                                        >  1:8395932/1‑93 (MQ=255)
                                  gAACGTCGCGGTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCCGCCTAATCGTTCCGCCGAAACGCTGGCAGGGgacgttacgg                                                       >  2:6065076/1‑92 (MQ=255)
                                  gAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgttacgt                                                       >  2:5464917/1‑93 (MQ=255)
                                  gAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgttacgt                                                       <  1:8604729/93‑1 (MQ=255)
                                    aCGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGa                                                     >  1:7452463/1‑93 (MQ=255)
                                      gTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCACGGGCAGTTATGTGACg                                                   >  1:6751715/1‑93 (MQ=255)
                                       tCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACgg                                                  <  1:4703486/93‑1 (MQ=255)
                                         gcgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGCAACGCTGGCAGGGTAGGTTACGTGACAGCg                                                >  1:5715877/1‑93 (MQ=255)
                                          cgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGt                                               <  1:5089521/93‑1 (MQ=255)
                                          cgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGt                                               >  1:6461528/1‑93 (MQ=255)
                                          cgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGGTACGTGACGGCGt                                               >  2:1668443/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGGAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  2:4624161/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  1:7379210/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             <  2:137612/93‑1 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  1:2534080/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTCGCACGGGACGCTACGTGACGGCGttt                                             >  1:8026334/1‑93 (MQ=255)
                                             tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTc                                            >  2:4743151/1‑93 (MQ=255)
                                              cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTcc                                           <  2:8407595/93‑1 (MQ=255)
                                              cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTcc                                           <  1:3800747/93‑1 (MQ=255)
                                               aCGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCt                                          <  2:7526664/93‑1 (MQ=255)
                                               aCGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGAAGTTACGTGACGGCGTTTCCt                                          >  2:8752398/1‑93 (MQ=255)
                                                cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTg                                         >  1:3368587/1‑93 (MQ=255)
                                                  ccAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGGTTCCTGGc                                       >  2:8642993/1‑93 (MQ=255)
                                                   cAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCt                                      <  2:4264154/93‑1 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTg                                     >  1:5968573/1‑93 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTg                                     <  1:1726781/93‑1 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGGTTCCTGGCTg                                     >  2:1221156/1‑93 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTg                                     >  2:59367/1‑93 (MQ=255)
                                                     aaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGGGTTTCCTGGCTGt                                    >  1:2410776/1‑93 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTTTTGaa                                >  2:7316426/1‑93 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                >  1:6380503/1‑93 (MQ=255)
                                                          atgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAAt                               >  2:7650021/1‑93 (MQ=255)
                                                            gatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTa                             <  2:1963955/93‑1 (MQ=255)
                                                             atTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTAt                            <  1:1147985/93‑1 (MQ=255)
                                                             atTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTAt                            >  1:4238602/1‑93 (MQ=255)
                                                               tGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                          <  1:7056830/93‑1 (MQ=255)
                                                               tGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                          <  1:5998379/93‑1 (MQ=255)
                                                                  ccTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgcc                       <  1:7142046/93‑1 (MQ=255)
                                                                  ccTGTCCGTGACGGCGGGCAACGTACTGGCACCGCCAAAGCGATGGCAGGGTACGTTACGTGACGGCGTTTCGTGGCTGTTGAATTATCTgcc                       <  2:5983556/93‑1 (MQ=255)
                                                                    tGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgccgc                     <  2:2767382/93‑1 (MQ=255)
                                                                      tCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCGGGGTACGTTACGTGGCGGCGTTTCCCGGCTGGTGAATTATCTTCCGACa                   >  2:1641545/1‑93 (MQ=255)
                                                                       ccGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                  <  1:1912012/93‑1 (MQ=255)
                                                                       ccGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                  <  2:2610830/93‑1 (MQ=255)
                                                                        cGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGGGTTTCCTGGCTGTTGAATTATCTGCCGCCAGt                 >  2:4475656/1‑93 (MQ=255)
                                                                        cGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGt                 <  1:6160321/93‑1 (MQ=255)
                                                                            aCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa             >  2:1161009/1‑93 (MQ=255)
                                                                            aCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa             <  2:8714113/93‑1 (MQ=255)
                                                                             cggcggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAAc            <  1:5612361/93‑1 (MQ=255)
                                                                             cggcggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAAc            <  2:403309/93‑1 (MQ=255)
                                                                               gcggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGc          <  2:8551516/93‑1 (MQ=255)
                                                                                cggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCt         >  2:4163412/1‑93 (MQ=255)
                                                                                 gggCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTa        >  2:1039712/1‑93 (MQ=255)
                                                                                   gCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACt      >  2:7869833/1‑93 (MQ=255)
                                                                                   gCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACactact      >  2:2195081/1‑93 (MQ=255)
                                                                                    cAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACtt     <  2:6938385/93‑1 (MQ=255)
                                                                                     aaCGTACTGGGTCCGCCGAAACGCTGGCAGGGGACGGTACGTGACGGCGGTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTCCTTc    >  2:6229836/1‑93 (MQ=255)
                                                                                     aaCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc    >  1:5458966/1‑93 (MQ=255)
                                                                                     aaCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc    >  2:3157171/1‑93 (MQ=255)
                                                                                      aCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc   >  1:4283662/1‑93 (MQ=255)
                                                                                      aCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc   >  2:2225854/1‑93 (MQ=255)
                                                                                       cGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAcctcct  <  1:2033726/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCT  >  NC_000913/368821‑369000

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