Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,205 A→G F345S (TTC→TCC)  yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,2050AG92.5% 1567.3 / 84.5 564F345S (TTC→TCC) yieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base G (262/260);  total (285/279)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 571 total aligned reads displayed.

ATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGAT  >  NC_000913/3937115‑3937294
                                                                                          |                                                                                         
aTCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGAc                                                                                         >  2:2172383/1‑93 (MQ=255)
   gcgGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGTAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGa                                                                                      <  2:5255800/93‑1 (MQ=255)
   gcgGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGa                                                                                      <  1:87547/93‑1 (MQ=255)
   gcgGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGa                                                                                      <  1:8094347/93‑1 (MQ=255)
     ggTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGc                                                                                    <  1:4957053/93‑1 (MQ=255)
      gTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCg                                                                                   >  1:8934911/1‑93 (MQ=255)
       tAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGa                                                                                  >  2:5994782/1‑93 (MQ=255)
         aaCTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATc                                                                                <  1:9857835/93‑1 (MQ=255)
           cTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCAt                                                                              <  1:7117124/93‑1 (MQ=255)
            tGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCAtt                                                                             <  2:1418013/93‑1 (MQ=255)
             gCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTa                                                                            >  2:6799/1‑93 (MQ=255)
             gCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTa                                                                            >  1:8433374/1‑93 (MQ=255)
               gTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGc                                                                          >  2:5562408/1‑93 (MQ=255)
                tCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGcc                                                                         >  2:3207620/1‑93 (MQ=255)
                 cGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGccc                                                                        <  1:2049077/93‑1 (MQ=255)
                  gACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCa                                                                       >  2:5886679/1‑93 (MQ=255)
                   aCATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCaa                                                                      <  2:5698421/93‑1 (MQ=255)
                    cATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCaat                                                                     <  1:1219087/93‑1 (MQ=255)
                     aTAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCaata                                                                    <  2:6557527/93‑1 (MQ=255)
                      tAGCCATCCCTAATATAAACAACGGTAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCaataa                                                                   >  1:4342547/1‑93 (MQ=255)
                      tAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCaataa                                                                   >  2:1722157/1‑93 (MQ=255)
                       aGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCaataat                                                                  <  2:7611884/93‑1 (MQ=255)
                        gCCATCCCTAATATAAACAACCGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATc                                                                 >  1:1844572/1‑93 (MQ=255)
                          cATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCAc                                                               <  2:7306989/93‑1 (MQ=255)
                           aTCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCAcc                                                              <  2:6072803/93‑1 (MQ=255)
                            tCCCTAATATAAATAACGGCAAGATCTGCATCCATATAGCCATTGCTGGTGATAGCAATGAGGACTGAGCGATCATTAGCCCAATAATCACCg                                                             >  2:8979767/1‑93 (MQ=255)
                            tCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCg                                                             <  1:5402200/93‑1 (MQ=255)
                            tCCATAATATAAACAACGGCAAGATCAGCATCAATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCg                                                             <  1:7255418/93‑1 (MQ=255)
                             cccTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGt                                                            >  1:5257348/1‑93 (MQ=255)
                             cccTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGt                                                            <  1:8583358/93‑1 (MQ=255)
                             cccTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGt                                                            >  2:8034271/1‑93 (MQ=255)
                              ccTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTg                                                           >  1:4846835/1‑93 (MQ=255)
                                   tataAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATccc                                                      >  2:2918744/1‑93 (MQ=255)
                                   tataAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATccc                                                      >  2:9889106/1‑93 (MQ=255)
                                   tataAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATccc                                                      >  1:5891959/1‑93 (MQ=255)
                                   aataAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATccc                                                      <  1:7874117/92‑1 (MQ=255)
                                     taAACAACGGCAAGATCAGCATCCATATAGCTATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCa                                                    <  2:1586392/93‑1 (MQ=255)
                                     taAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCa                                                    <  1:6613315/93‑1 (MQ=255)
                                     taAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCa                                                    >  2:10013328/1‑93 (MQ=255)
                                     taAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCa                                                    >  1:9444707/1‑93 (MQ=255)
                                      aaaCAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCAc                                                   <  2:6409268/93‑1 (MQ=255)
                                        acaacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTa                                                 <  1:1103317/93‑1 (MQ=255)
                                         caacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTaa                                                <  1:2644285/93‑1 (MQ=255)
                                         caacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTaa                                                <  1:4419358/93‑1 (MQ=255)
                                         caacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTaa                                                >  2:2536445/1‑93 (MQ=255)
                                          aacGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAAc                                               >  2:3814114/1‑93 (MQ=255)
                                             ggCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTa                                            <  1:9269149/93‑1 (MQ=255)
                                             ggCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGCGCGATCATTAGCCCAATAATCACCAAAATGCCCACTAACGTa                                            >  1:5747935/1‑93 (MQ=255)
                                              gCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTAt                                           >  2:8482892/1‑93 (MQ=255)
                                              gCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTAt                                           <  1:3738927/93‑1 (MQ=255)
                                              gCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTAt                                           <  1:4060231/93‑1 (MQ=255)
                                              gCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTAt                                           <  2:1624236/93‑1 (MQ=255)
                                               taaGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATg                                          <  2:1988854/92‑1 (MQ=255)
                                                   aTCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGa                                      >  1:602734/1‑93 (MQ=255)
                                                   aTCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGa                                      >  2:2365105/1‑93 (MQ=255)
                                                    tCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGAt                                     <  1:2781771/93‑1 (MQ=255)
                                                    tCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGAt                                     >  2:4664768/1‑93 (MQ=255)
                                                    tCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGAt                                     <  1:2987400/93‑1 (MQ=255)
                                                     cAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATa                                    <  1:6509245/93‑1 (MQ=255)
                                                       gCATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGc                                  >  2:397555/1‑93 (MQ=255)
                                                        cATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGcc                                 >  2:9570848/1‑93 (MQ=255)
                                                        cATCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGcc                                 <  2:3628374/93‑1 (MQ=255)
                                                          tCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCa                               >  1:8215665/1‑93 (MQ=255)
                                                          tCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCa                               <  1:3322550/93‑1 (MQ=255)
                                                          tCCATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCa                               >  1:7685116/1‑93 (MQ=255)
                                                            cATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAga                             <  1:3567300/93‑1 (MQ=255)
                                                            cATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAga                             <  2:1206976/93‑1 (MQ=255)
                                                            cATATAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAga                             <  2:4492535/93‑1 (MQ=255)
                                                             atatAGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgac                            <  1:4223315/93‑1 (MQ=255)
                                                              tataGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacg                           >  1:3100168/1‑93 (MQ=255)
                                                              tataGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacg                           >  1:5621678/1‑93 (MQ=255)
                                                               ataGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacga                          <  2:8742196/93‑1 (MQ=255)
                                                                taGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacgac                         >  2:5066454/1‑93 (MQ=255)
                                                                taGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacgac                         >  2:5459895/1‑93 (MQ=255)
                                                                taGCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAgacgac                         >  2:2698082/1‑93 (MQ=255)
                                                                  gCCATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGt                       >  1:4661170/1‑93 (MQ=255)
                                                                   ccATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTa                      >  2:422465/1‑93 (MQ=255)
                                                                    cATTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTaa                     >  2:4989253/1‑93 (MQ=255)
                                                                     aTTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAg                    <  1:1720006/93‑1 (MQ=255)
                                                                     aTTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAg                    <  2:7765245/93‑1 (MQ=255)
                                                                     aTTGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAg                    <  1:1497991/93‑1 (MQ=255)
                                                                      ttGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGa                   <  2:8879469/93‑1 (MQ=255)
                                                                      ttGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGa                   <  2:7098550/93‑1 (MQ=255)
                                                                      ttGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGa                   >  2:9763674/1‑93 (MQ=255)
                                                                      ttGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGa                   >  1:6921520/1‑93 (MQ=255)
                                                                       tGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGAc                  <  2:5764278/93‑1 (MQ=255)
                                                                       tGCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGAc                  >  1:8995952/1‑93 (MQ=255)
                                                                        gCTGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACt                 >  1:2496405/1‑93 (MQ=255)
                                                                        gCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACt                 <  1:6254246/93‑1 (MQ=255)
                                                                          tGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTg               >  1:2158235/1‑93 (MQ=255)
                                                                          tGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTg               <  2:8342465/93‑1 (MQ=255)
                                                                          tGGTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTg               <  2:9238612/93‑1 (MQ=255)
                                                                           ggTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTgg              >  2:3469212/1‑93 (MQ=255)
                                                                            gTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTggg             <  1:4549943/93‑1 (MQ=255)
                                                                            gTGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTggg             >  1:5221335/1‑93 (MQ=255)
                                                                             tGATTGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTggga            >  1:6765643/1‑92 (MQ=255)
                                                                                ttGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAAc         >  1:7507374/1‑93 (MQ=255)
                                                                                ttGCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAAc         >  2:625122/1‑93 (MQ=255)
                                                                                  gCAAAGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGAGCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCa       >  1:3969855/1‑93 (MQ=255)
                                                                                     aaGAGGACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCg    >  1:9201979/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTAACGTATGGCGATAGCCCAGACGACGTAAGACTTGGGTAACCATCGAT  >  NC_000913/3937115‑3937294

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