Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,240 (G)7→8 intergenic (+18/‑165) cstA → / → ybdD carbon starvation protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,2331.G83.4% 618.7 / 114.3 199intergenic (+11/‑172)cstA/ybdDcarbon starvation protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (22/11);  new base G (100/66);  total (122/77)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 239 total aligned reads displayed.

GCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGC  >  NC_000913/631141‑631322
                                                                                             |                                                                                           
gCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGa                                                                                              <  2:5771917/93‑1 (MQ=255)
   aaCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggg                                                                                        >  2:6099402/1‑93 (MQ=255)
   aaCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggg                                                                                        <  1:1801126/93‑1 (MQ=255)
    aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggg                                                                                       <  2:1422539/93‑1 (MQ=255)
    aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggg                                                                                       <  1:4444797/93‑1 (MQ=255)
    aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggg                                                                                       <  2:1926892/93‑1 (MQ=255)
    aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggg                                                                                       >  1:1011523/1‑93 (MQ=255)
     cGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggggg                                                                                      <  1:2413390/93‑1 (MQ=255)
     cGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggggg                                                                                      >  2:3464078/1‑93 (MQ=255)
      ggCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggggg                                                                                     <  2:6579166/93‑1 (MQ=255)
      ggCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑gggggg                                                                                     <  1:1271375/93‑1 (MQ=255)
         gaaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggggggggg                                                                                  >  2:4643945/1‑91 (MQ=255)
          aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑ggggggggat                                                                                 <  1:2287282/93‑4 (MQ=255)
             gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGTAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGtt                                                                               <  2:6394066/93‑1 (MQ=255)
             gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTg                                                                              >  2:5790406/1‑93 (MQ=255)
              aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTg                                                                              >  2:1849577/1‑93 (MQ=255)
                aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTggg                                                                            >  2:5448557/1‑92 (MQ=255)
                aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGc                                                                            >  1:3908443/1‑93 (MQ=255)
                aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCg                                                                           >  2:4211571/1‑93 (MQ=255)
                 cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCg                                                                           >  2:2228434/1‑93 (MQ=255)
                  gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGGGa                                                                          >  2:1527649/1‑93 (MQ=255)
                    cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGTGTTGGCGAAt                                                                        >  1:1434638/1‑93 (MQ=255)
                    cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAAt                                                                        >  1:2870598/1‑93 (MQ=255)
                      tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGt                                                                      <  1:3092142/93‑1 (MQ=255)
                         gAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGATGTTGGCGAATGTTg                                                                    <  2:5887697/93‑1 (MQ=255)
                          aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGc                                                                  >  1:5891760/1‑93 (MQ=255)
                            ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCtt                                                                >  1:2197119/1‑93 (MQ=255)
                             cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTa                                                               <  2:2892862/93‑1 (MQ=255)
                             cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTa                                                               >  1:414036/1‑93 (MQ=255)
                              aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAg                                                              >  2:2515420/1‑93 (MQ=255)
                              aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAg                                                              <  2:1285003/93‑1 (MQ=255)
                              aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGTTGTTGGCGAATGTTGGCTTAGt                                                             >  1:4502165/1‑93 (MQ=255)
                              aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGGGAGGGGGGGGGTGTTGGCGAATGTTGGCTTAg                                                              >  1:2486526/1‑93 (MQ=255)
                               aTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGt                                                             <  1:3695971/93‑1 (MQ=255)
                                tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTg                                                            >  1:3325328/1‑93 (MQ=255)
                                 gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGc                                                           <  2:1654046/93‑1 (MQ=255)
                                  ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGcc                                                          >  1:1395622/1‑93 (MQ=255)
                                   cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGccc                                                         <  2:5545706/93‑1 (MQ=255)
                                     gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                      >  1:1643342/1‑93 (MQ=255)
                                      aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                      <  2:2067378/93‑1 (MQ=255)
                                       aaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                     <  1:5185344/93‑1 (MQ=255)
                                       aaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                     >  2:2212023/1‑93 (MQ=255)
                                         aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGtt                                                   <  2:5250954/93‑1 (MQ=255)
                                            tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                <  1:4297621/93‑1 (MQ=255)
                                            tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                <  1:2404652/93‑1 (MQ=255)
                                             cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCct                                               >  1:205657/1‑93 (MQ=255)
                                              gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                              >  1:1386255/1‑93 (MQ=255)
                                              gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                              >  2:2734466/1‑93 (MQ=255)
                                               aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct                                             <  1:5321589/93‑1 (MQ=255)
                                                 gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctcc                                           >  2:819645/1‑92 (MQ=255)
                                                 gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                           >  2:1088009/1‑93 (MQ=255)
                                                  agaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAc                                          >  2:4936432/1‑93 (MQ=255)
                                                  agaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAc                                          >  2:787067/1‑93 (MQ=255)
                                                   gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                         >  1:63240/1‑93 (MQ=255)
                                                   gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                         <  2:3249662/93‑1 (MQ=255)
                                                    aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                        >  2:5518178/1‑93 (MQ=255)
                                                    aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                        >  2:3672938/1‑93 (MQ=255)
                                                    aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                        >  1:556215/1‑93 (MQ=255)
                                                     tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCt                                       >  1:4924846/1‑93 (MQ=255)
                                                     tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                      >  1:3115436/1‑93 (MQ=255)
                                                         ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAcc                                   <  2:2842071/93‑1 (MQ=255)
                                                           cgcAGGCAAAAGGCGCACCCTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                 <  1:5390909/93‑1 (MQ=255)
                                                             cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctct                               <  1:3842356/93‑1 (MQ=255)
                                                             cAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCtctc                              >  2:5999269/1‑93 (MQ=255)
                                                              aGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                             >  1:6440030/1‑93 (MQ=255)
                                                               ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                             >  1:1621193/1‑93 (MQ=255)
                                                                gCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                            <  2:5315850/93‑1 (MQ=255)
                                                                 cAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcccc                           <  1:1849577/93‑1 (MQ=255)
                                                                 cAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GCGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCGCCCTAACCCTCTcccc                           <  1:5435013/93‑1 (MQ=255)
                                                                  aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCg                          >  1:5950933/1‑93 (MQ=255)
                                                                     aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                       <  1:2529024/93‑1 (MQ=255)
                                                                      ggCGCACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                      >  1:5039275/1‑93 (MQ=255)
                                                                        cgcACACTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                    >  1:3766114/1‑93 (MQ=255)
                                                                         gcACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                  <  1:3170408/93‑1 (MQ=255)
                                                                          cacaCTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                  >  2:4296139/1‑93 (MQ=255)
                                                                           acacTAAAGTCAGAGTGAGGGGGGGGATGTTAGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                 >  2:5127188/1‑93 (MQ=255)
                                                                            cacTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                <  1:4370323/93‑1 (MQ=255)
                                                                            cacTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               >  1:960397/1‑93 (MQ=255)
                                                                             acTGAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               >  1:2230580/1‑93 (MQ=255)
                                                                             acTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTTCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               >  2:2359603/1‑93 (MQ=255)
                                                                             acTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               >  1:2778360/1‑93 (MQ=255)
                                                                             acTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               <  2:5865146/93‑1 (MQ=255)
                                                                             acTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg               >  1:6074234/1‑93 (MQ=255)
                                                                              cTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg              <  2:4424278/93‑1 (MQ=255)
                                                                              cTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg              >  2:1667553/1‑93 (MQ=255)
                                                                              cTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg              <  2:6372206/93‑1 (MQ=255)
                                                                              cTAAAGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg              >  1:348736/1‑93 (MQ=255)
                                                                                aaaGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc           <  2:2924147/93‑1 (MQ=255)
                                                                                 aaGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc           <  1:6032352/93‑1 (MQ=255)
                                                                                 aaGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc           <  2:3063577/93‑1 (MQ=255)
                                                                                 aaGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc           <  2:2469477/93‑1 (MQ=255)
                                                                                  aGTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt          >  1:2647556/1‑93 (MQ=255)
                                                                                  aGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg         >  1:4058654/1‑93 (MQ=255)
                                                                                   gTCAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg         >  1:768447/1‑93 (MQ=255)
                                                                                     cAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc       <  2:353521/93‑1 (MQ=255)
                                                                                     cAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc       <  1:5544182/93‑1 (MQ=255)
                                                                                     cAGAGTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc       <  1:1200330/93‑1 (MQ=255)
                                                                                       gagTGAGGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg     >  2:6254088/1‑93 (MQ=255)
                                                                                       gagTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa    <  2:3810854/93‑1 (MQ=255)
                                                                                        agTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg   >  1:5737395/1‑93 (MQ=255)
                                                                                             |                                                                                           
GCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑‑GGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGC  >  NC_000913/631141‑631322

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