Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,362,693 T→C L751P (CTC→CCC)  yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,6930TC99.6% 1830.7 / ‑4.6 556L751P (CTC→CCC) yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (254/300);  total (254/302)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 565 total aligned reads displayed.

CCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCA  >  NC_000913/3362603‑3362781
                                                                                          |                                                                                        
  gACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAAt                                                                                      <  2:1326844/93‑1 (MQ=255)
     ggCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAAc                                                                                   >  2:6908061/1‑93 (MQ=255)
       cAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACAt                                                                                 <  2:9242858/93‑1 (MQ=255)
       cAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACAt                                                                                 <  2:2782443/93‑1 (MQ=255)
        aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                <  2:4468550/93‑1 (MQ=255)
        aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                <  1:517660/93‑1 (MQ=255)
         aaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCCACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCa                                                                               >  2:1865551/1‑93 (MQ=255)
         aaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCa                                                                               <  1:13328/93‑1 (MQ=255)
         aaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCa                                                                               <  2:8109156/93‑1 (MQ=255)
           aaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCaaa                                                                             <  1:1199792/93‑1 (MQ=255)
           aaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCaaa                                                                             <  2:5302958/93‑1 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                            <  2:1818214/93‑1 (MQ=255)
             aCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAcc                                                                           <  1:1392187/93‑1 (MQ=255)
             aCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTAGTATTCCCCAATAACATCAAAcc                                                                           >  2:5497534/1‑93 (MQ=255)
              cTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCg                                                                          <  2:161780/93‑1 (MQ=255)
              cTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCg                                                                          <  1:7123836/93‑1 (MQ=255)
               tGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCgg                                                                         <  2:6106089/93‑1 (MQ=255)
                gCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGa                                                                        >  1:2353169/1‑93 (MQ=255)
                gCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGa                                                                        <  1:9545666/93‑1 (MQ=255)
                 cccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGAc                                                                       <  2:397942/93‑1 (MQ=255)
                 cccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGAc                                                                       <  1:9981972/93‑1 (MQ=255)
                  ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                      >  1:7696060/1‑93 (MQ=255)
                   cAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACAt                                                                     >  1:8711735/1‑93 (MQ=255)
                    agagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATg                                                                    <  2:9672045/93‑1 (MQ=255)
                    agagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATg                                                                    >  1:2398136/1‑93 (MQ=255)
                    agagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATg                                                                    >  1:8497075/1‑93 (MQ=255)
                     gagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                   <  2:9879282/93‑1 (MQ=255)
                      agGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGt                                                                  <  2:3216854/93‑1 (MQ=255)
                       ggCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTg                                                                 >  1:6700741/1‑93 (MQ=255)
                        gCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGc                                                                <  2:9339307/93‑1 (MQ=255)
                        gCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGc                                                                <  2:6476585/93‑1 (MQ=255)
                           ttGCGATTGAAGATATTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAg                                                             >  1:2681640/1‑93 (MQ=255)
                           ttGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAg                                                             >  2:7629219/1‑93 (MQ=255)
                           ttGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAg                                                             >  1:455920/1‑93 (MQ=255)
                             gCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGat                                                           <  2:4999498/93‑1 (MQ=255)
                                atggAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                        <  1:677832/90‑1 (MQ=255)
                                aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                        >  2:2066682/1‑93 (MQ=255)
                                aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                        >  1:3233363/1‑93 (MQ=255)
                                 ttGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaa                                                       <  2:2508358/93‑1 (MQ=255)
                                 ttGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCAAGATATaa                                                       >  1:7786756/1‑93 (MQ=255)
                                    aaGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAg                                                    >  2:3512442/1‑93 (MQ=255)
                                     aGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGa                                                   <  1:1927101/93‑1 (MQ=255)
                                     aGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGa                                                   <  1:360172/93‑1 (MQ=255)
                                     aGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGa                                                   >  1:6250319/1‑93 (MQ=255)
                                      gATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGAt                                                  >  2:2211851/1‑93 (MQ=255)
                                      gATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGAt                                                  >  2:5047864/1‑93 (MQ=255)
                                        tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                <  2:575012/93‑1 (MQ=255)
                                        tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                >  2:586601/1‑93 (MQ=255)
                                        tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGAACATGGTGCTAGATATAAAAGATGa                                                <  2:668795/93‑1 (MQ=255)
                                         aGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGag                                               <  2:2595110/93‑1 (MQ=255)
                                          gTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagc                                              <  1:3326712/93‑1 (MQ=255)
                                            gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                            >  2:3920013/1‑93 (MQ=255)
                                            gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                            >  1:8634823/1‑93 (MQ=255)
                                            gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                            >  2:7763212/1‑93 (MQ=255)
                                            gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGGTGTATTCCCCAATAACATCAAACCGGACATGGGGCTAGATATAAAAGATGagcag                                            >  1:3857532/1‑93 (MQ=255)
                                             aaGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcagc                                           <  1:5270158/93‑1 (MQ=255)
                                             aaGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcagc                                           >  1:5348174/1‑93 (MQ=255)
                                               ggCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCaa                                         >  2:2354686/1‑93 (MQ=255)
                                               ggCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCaa                                         >  2:7106446/1‑93 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                      <  2:9463409/93‑1 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                      <  2:8541912/93‑1 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                      <  2:8289000/93‑1 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                      <  2:6245014/93‑1 (MQ=255)
                                                    ttATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATg                                    >  2:6162642/2‑93 (MQ=255)
                                                    cTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATg                                    >  1:8368588/1‑93 (MQ=255)
                                                     tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                   <  1:9693098/93‑1 (MQ=255)
                                                      aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                  <  2:8757665/93‑1 (MQ=255)
                                                      aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                  <  2:8499623/93‑1 (MQ=255)
                                                       tCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCg                                 >  1:6547365/1‑93 (MQ=255)
                                                         tGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTa                               <  1:3499872/93‑1 (MQ=255)
                                                            caacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTc                            <  1:9142541/93‑1 (MQ=255)
                                                             aacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCa                           >  1:8048836/1‑93 (MQ=255)
                                                              acaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAc                          >  1:4086566/1‑93 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                       <  2:7410964/93‑1 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                       >  1:3779127/1‑93 (MQ=255)
                                                                  cAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTa                      <  1:3006658/93‑1 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                    <  1:9345516/93‑1 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                    >  2:4811495/1‑93 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                    <  1:1456675/93‑1 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                    >  2:332891/1‑93 (MQ=255)
                                                                      ggATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACAtt                  <  2:7314900/93‑1 (MQ=255)
                                                                       gATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTc                 <  1:5574042/93‑1 (MQ=255)
                                                                            cGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCaga            <  2:41473/93‑1 (MQ=255)
                                                                            cGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCaga            <  2:3688451/93‑1 (MQ=255)
                                                                             ggTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaa           >  2:310742/1‑93 (MQ=255)
                                                                             ggTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaa           >  1:595973/1‑93 (MQ=255)
                                                                             cgtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaa           <  2:2960363/92‑1 (MQ=255)
                                                                              gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag          >  1:628250/1‑93 (MQ=255)
                                                                              gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag          >  2:9080066/1‑93 (MQ=255)
                                                                               tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga         <  2:1213920/93‑1 (MQ=255)
                                                                               tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga         >  1:572983/1‑93 (MQ=255)
                                                                               ggttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga         <  2:463183/92‑1 (MQ=255)
                                                                                gttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGAt        >  2:7962671/1‑93 (MQ=255)
                                                                                 ttgtATTGCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCGGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATg       <  1:1806430/93‑1 (MQ=255)
                                                                                 ttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATg       <  1:1016902/93‑1 (MQ=255)
                                                                                 ttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATg       <  1:4915507/93‑1 (MQ=255)
                                                                                    tATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCAc    <  1:7387019/93‑1 (MQ=255)
                                                                                    tATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCAc    >  1:2840147/1‑93 (MQ=255)
                                                                                    tATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCAc    <  1:2022524/93‑1 (MQ=255)
                                                                                      ttCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCa  >  1:894625/1‑93 (MQ=255)
                                                                                          |                                                                                        
CCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCA  >  NC_000913/3362603‑3362781

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