Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,305,259 T→C N187D (AAT→GAT)  mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,305,2590TC98.2% 362.0 / ‑2.4 109N187D (AAT→GAT) mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (59/48);  total (60/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 109 total aligned reads displayed.

CGTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATTGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGCACAA  >  NC_000913/3305163‑3305354
                                                                                                |                                                                                               
cGTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCt                                                                                              >  2:1735870/1‑100 (MQ=255)
cGTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCt                                                                                              >  2:2485156/1‑100 (MQ=255)
 gTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCtt                                                                                             >  1:256518/1‑100 (MQ=255)
 gTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCtt                                                                                             >  1:509513/1‑100 (MQ=255)
    aTACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATGAACAGATACGGTGCATAAGACGCAGCGCTTTCg                                                                                          >  2:404177/1‑100 (MQ=255)
    aTACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCg                                                                                          <  2:1404128/100‑1 (MQ=255)
     tACCTCATCAGGCTCGGCACGTTACCCTGATAACCAAACGACGCCAGACAGAACGCCACGATCATCAACAGATACGGGGCATAAGACGCATCACTTTCgg                                                                                         <  2:664801/100‑1 (MQ=255)
      aCTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGc                                                                                        <  2:334712/100‑1 (MQ=255)
      aCTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGc                                                                                        >  1:590555/1‑100 (MQ=255)
      aCTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGc                                                                                        <  1:1521363/100‑1 (MQ=255)
       cTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCg                                                                                       >  1:63537/1‑100 (MQ=255)
         tcatcaGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGAc                                                                                     <  1:64317/100‑1 (MQ=255)
          catcaGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACg                                                                                    >  1:2471608/1‑100 (MQ=255)
          catcaGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACg                                                                                    >  2:1857305/1‑100 (MQ=255)
           atcaGGCTTGGGACGTTCCCGTGATAACCAAACGATGCCAGACAGAACGTCAGGGTCATGAACAGATACGGTGCATACGGCGCATCGCTTTCGGCGACgc                                                                                   <  2:756423/100‑2 (MQ=255)
           agcaGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACATATACGGTGCATAAGACGCATCGCTTTCGGCGACGt                                                                                   <  1:2184637/98‑1 (MQ=255)
            tcaGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGtt                                                                                  >  1:2566201/1‑100 (MQ=255)
              aGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTgn                                                                                <  1:2093421/100‑2 (MQ=255)
               ggCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGaa                                                                               >  2:1215707/1‑100 (MQ=255)
               ggCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGaa                                                                               >  1:903382/1‑100 (MQ=255)
                gCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAAc                                                                              <  1:548802/100‑1 (MQ=255)
                gCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAAc                                                                              <  2:1899790/100‑1 (MQ=255)
                   tGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAAt                                                                           >  1:1425489/1‑100 (MQ=255)
                    ggCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATg                                                                          >  1:2224921/1‑100 (MQ=255)
                    ggCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATg                                                                          <  2:1415110/100‑1 (MQ=255)
                    ggCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATg                                                                          >  2:121988/1‑100 (MQ=255)
                     gCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACACATCGCTTTCGGCGACGTTGAACAATGt                                                                         >  1:845125/1‑100 (MQ=255)
                      cACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTc                                                                        >  1:474024/1‑100 (MQ=255)
                       aCGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCg                                                                       >  2:154254/1‑100 (MQ=255)
                           tACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACATTGTTCAACGTCGccga                                                                   <  2:685122/100‑4 (MQ=255)
                             ccGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTGGCAGGCt                                                                 >  1:2604918/1‑100 (MQ=255)
                             ccGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCt                                                                 >  1:2591737/1‑100 (MQ=255)
                              cGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTg                                                                >  1:2595477/1‑100 (MQ=255)
                               gTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGc                                                               <  2:407307/100‑1 (MQ=255)
                               gTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCCTTCGGCGACGTTGAACAATGTCGCAGGCTGc                                                               >  2:1647326/1‑100 (MQ=255)
                                tGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGca                                                              <  2:82710/100‑1 (MQ=255)
                                tGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGca                                                              <  1:173453/100‑1 (MQ=255)
                                 gATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGcac                                                             <  1:1801700/100‑1 (MQ=255)
                                 gATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGcac                                                             <  2:2123944/100‑1 (MQ=255)
                                  aTAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGcaca                                                            <  1:2496977/100‑1 (MQ=255)
                                  aTAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGcaca                                                            <  1:2230739/100‑1 (MQ=255)
                                  aTAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGcaca                                                            <  2:1279272/100‑1 (MQ=255)
                                   tAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATTGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACAt                                                           <  1:1512768/100‑1 (MQ=255)
                                   tAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCGGCACAt                                                           >  2:1167482/1‑100 (MQ=255)
                                     aCCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGc                                                         >  1:1654920/1‑100 (MQ=255)
                                      ccAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCGGGCTGCACATGcc                                                        >  2:701034/1‑100 (MQ=255)
                                      ccAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGcc                                                        >  1:521079/1‑100 (MQ=255)
                                       cAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGccc                                                       >  1:2420746/1‑100 (MQ=255)
                                       cAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCATCCGGCGACGTTGAACAATGTCGCCGGCTGCACATGcac                                                       >  2:1403762/1‑98 (MQ=255)
                                           cGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCcagc                                                   >  2:1651458/1‑100 (MQ=255)
                                            gATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCcagca                                                  >  1:888296/1‑100 (MQ=255)
                                             aTGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCcagcag                                                 <  1:1010759/100‑1 (MQ=255)
                                              tgacagacagaACGGCAGGGTCATCAACAGATTCGTTTCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAgg                                                >  2:2140151/4‑100 (MQ=255)
                                              tGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAgg                                                <  1:825033/100‑1 (MQ=255)
                                               gCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGc                                               <  1:1127414/100‑1 (MQ=255)
                                               gCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGc                                               >  1:1431331/1‑100 (MQ=255)
                                                ccAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCt                                              >  2:137474/1‑100 (MQ=255)
                                                 cagacagaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTa                                             >  2:2301849/1‑100 (MQ=255)
                                                   gacagaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTAcc                                           >  1:1855477/1‑100 (MQ=255)
                                                     cagaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCaa                                         >  1:1355883/1‑100 (MQ=255)
                                                       gaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATTGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAg                                       >  1:1685319/1‑100 (MQ=255)
                                                       gaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGGCGTTGAACAATGGCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAg                                       >  1:1466532/1‑100 (MQ=255)
                                                       gaACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAg                                       >  1:1402311/1‑100 (MQ=255)
                                                        aaCGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAgg                                      <  1:2366725/100‑1 (MQ=255)
                                                         aCGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGt                                     <  2:940830/100‑1 (MQ=255)
                                                         aCGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGt                                     <  1:100272/100‑1 (MQ=255)
                                                         aCGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGt                                     >  2:2162353/1‑100 (MQ=255)
                                                          cGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTg                                    <  2:231/100‑1 (MQ=255)
                                                          cGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTg                                    >  2:1372192/1‑100 (MQ=255)
                                                              aGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGga                                <  1:1892951/100‑1 (MQ=255)
                                                              aGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGga                                <  1:1454016/100‑1 (MQ=255)
                                                              aGGGTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGga                                >  1:1907359/1‑100 (MQ=255)
                                                                 gTCATCAACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGgaaga                             >  2:317540/1‑100 (MQ=255)
                                                                  tcatcaACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGgaagaa                            <  1:2540292/100‑1 (MQ=255)
                                                                     tcaACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGt                         >  1:1842005/1‑100 (MQ=255)
                                                                      caACAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTa                        >  1:1215454/1‑100 (MQ=255)
                                                                       aaCAGATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTaa                       <  1:2020097/100‑1 (MQ=255)
                                                                         cAGCTACGGTGCATAAGACGCATCGCGTTCGGCGACGGTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATg                     <  2:2482892/100‑1 (MQ=255)
                                                                          aGATACAGTGCATAAGACGCATCGCTTTCGTGGACGTTGAACAATATCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGCTCATGa                    <  2:2441996/100‑1 (MQ=255)
                                                                           gATACGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGAc                   <  2:675397/100‑1 (MQ=255)
                                                                              aCGGCGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGATGCACATGCCCCAGCAGGCCACCAAAGGTGAGGAAGAAGGTAAGGACttt                >  2:2444000/1‑100 (MQ=255)
                                                                               cGGTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTc               >  2:491534/1‑100 (MQ=255)
                                                                                ggTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAGGGTAATGACTTTCg              <  1:1200565/100‑1 (MQ=255)
                                                                                ggTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCg              >  1:1784322/1‑100 (MQ=255)
                                                                                ggTGCATAAGACACATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCg              <  2:830450/100‑1 (MQ=255)
                                                                                 gTGCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGc             >  1:2024058/1‑100 (MQ=255)
                                                                                   gCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGccc           <  1:2531919/100‑1 (MQ=255)
                                                                                   gCATAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGccc           <  2:1275431/100‑1 (MQ=255)
                                                                                     aTAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGccccc         >  2:1862346/1‑100 (MQ=255)
                                                                                     aTAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGccccc         >  1:2151386/1‑100 (MQ=255)
                                                                                      tAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCa        >  2:1434206/1‑100 (MQ=255)
                                                                                      tAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCa        <  1:1481525/100‑1 (MQ=255)
                                                                                      tAAGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCa        <  1:1481404/100‑1 (MQ=255)
                                                                                       aaGACGCATCGCTTTCGGCGACGTTGACCAATGTCGCAGGCTGCACATGCCCCAGCAGGGTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAg       <  1:2290957/100‑1 (MQ=255)
                                                                                       aaGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAg       >  2:2366842/1‑100 (MQ=255)
                                                                                        aGACGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGc      >  2:884525/1‑100 (MQ=255)
                                                                                           cGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGcaca   <  2:67459/100‑1 (MQ=255)
                                                                                           cGCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGcaca   >  1:2004618/1‑100 (MQ=255)
                                                                                            gCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGCACaa  >  1:1192569/1‑100 (MQ=255)
                                                                                            gCATCGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGCACaa  <  2:2052861/100‑1 (MQ=255)
                                                                                                |                                                                                               
CGTAATACTTCATCAGGCTTGGCACGTTACCGTGATAACCAAACGATGCCAGACAGAACGGCAGGGTCATCAACAGATACGGTGCATAAGACGCATTGCTTTCGGCGACGTTGAACAATGTCGCAGGCTGCACATGCCCCAGCAGGCTACCAAAGGTGAGGAAGAAGGTAATGACTTTCGCCCCCAGCACAA  >  NC_000913/3305163‑3305354

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