Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 144,939 G→A A121A (GCG→GCA yadI → putative PTS enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,9390GA96.3% 620.2 / 4.3 192A121A (GCG→GCAyadIputative PTS enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base A (89/95);  total (89/103)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

Only 100 of 201 total aligned reads displayed.

TGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTC  >  NC_000913/144849‑145027
                                                                                          |                                                                                         
tgtCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCAcc                                                                                         >  1:2767601/1‑93 (MQ=255)
tgtCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGCATTGTCGAACTGGGTGCAcc                                                                                         >  1:1885897/1‑93 (MQ=255)
  tCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCgg                                                                                       >  1:3690515/1‑93 (MQ=255)
    aCGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGCTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAg                                                                                     >  1:592896/1‑93 (MQ=255)
      gTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAggg                                                                                   >  2:468064/1‑92 (MQ=255)
      gTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAggg                                                                                   >  2:3198840/1‑92 (MQ=255)
      gTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGt                                                                                   <  2:4269187/93‑1 (MQ=255)
      gTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGt                                                                                   <  2:1311108/93‑1 (MQ=255)
        tACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGa                                                                                 <  1:810277/93‑1 (MQ=255)
        tACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGa                                                                                 <  2:3231216/93‑1 (MQ=255)
          ccGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagt                                                                               >  1:3388711/1‑93 (MQ=255)
           cGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagta                                                                              <  2:3334872/93‑1 (MQ=255)
            gTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGCGGTGagtag                                                                             >  2:3662256/1‑93 (MQ=255)
            gTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtag                                                                             >  2:374793/1‑93 (MQ=255)
            gTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtag                                                                             <  1:3220379/93‑1 (MQ=255)
            gTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtag                                                                             >  2:2500085/1‑93 (MQ=255)
            gTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtag                                                                             >  2:2688485/1‑93 (MQ=255)
             ttAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtagt                                                                            <  2:492314/93‑1 (MQ=255)
             ttAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTACAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGagtagt                                                                            >  2:3615272/1‑93 (MQ=255)
              tAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTc                                                                           >  2:3879632/1‑93 (MQ=255)
                aTTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCtt                                                                         >  1:139199/1‑93 (MQ=255)
                aTTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCtt                                                                         <  2:1453710/93‑1 (MQ=255)
                 ttGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGGGCACCGGAGGGGAGTAGTCttt                                                                        >  2:4393017/1‑93 (MQ=255)
                  tGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTg                                                                       <  2:441319/93‑1 (MQ=255)
                    aaCAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGc                                                                     <  2:3008682/93‑1 (MQ=255)
                     aCAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTggac                                                                    >  1:1425008/1‑91 (MQ=255)
                      cAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCAc                                                                   <  2:2078203/93‑1 (MQ=255)
                      cAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCAc                                                                   <  2:1336861/93‑1 (MQ=255)
                      cAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCAc                                                                   <  1:3326124/93‑1 (MQ=255)
                       aGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCAcc                                                                  >  2:1878541/1‑93 (MQ=255)
                       aGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCAca                                                                  >  2:3898217/1‑92 (MQ=255)
                        gatgatGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACca                                                                 >  2:230676/1‑93 (MQ=255)
                        gatgatGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACca                                                                 <  2:2861694/93‑1 (MQ=255)
                         atgatgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaa                                                                <  2:446563/93‑1 (MQ=255)
                           gatgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaaca                                                              >  2:245230/1‑93 (MQ=255)
                           gatgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaaca                                                              >  1:1621593/1‑93 (MQ=255)
                            atgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaa                                                             <  1:2362999/93‑1 (MQ=255)
                             tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACcaacaac                                                            <  1:3528168/93‑1 (MQ=255)
                             tgGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGGAGTCTTTGGCACcaacaaa                                                            >  1:2459183/1‑92 (MQ=255)
                              ggCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaaca                                                          >  2:2474089/1‑93 (MQ=255)
                               gCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaacaa                                                         <  2:2234523/93‑1 (MQ=255)
                               gCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACcaacaacaa                                                         <  2:3434178/93‑1 (MQ=255)
                                 ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                       >  2:204282/1‑93 (MQ=255)
                                 ttGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaa                                                       >  2:4159797/1‑93 (MQ=255)
                                    ccGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaa                                                    <  1:3797645/93‑1 (MQ=255)
                                    ccGTGAAACCATAACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaa                                                    <  1:110078/93‑1 (MQ=255)
                                     cGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaa                                                   <  1:2817314/93‑2 (MQ=255)
                                     cGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaa                                                   >  2:597602/1‑92 (MQ=255)
                                      gTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACaaaaaaaat                                                  >  1:2261021/1‑91 (MQ=255)
                                       tGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGGGAGTAGTCTTTGGCACCAACAACaaaaaaaatc                                                 >  2:1031683/1‑90 (MQ=255)
                                          aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                               >  2:794891/1‑93 (MQ=255)
                                          aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                               >  2:2196237/1‑93 (MQ=255)
                                          aaCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGc                                               <  1:2124815/93‑1 (MQ=255)
                                           aCCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGcc                                              >  1:3342104/1‑93 (MQ=255)
                                            ccATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCt                                             <  2:3420935/93‑1 (MQ=255)
                                             cATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCtt                                            >  2:2938394/1‑93 (MQ=255)
                                             cATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCtt                                            >  1:1678219/1‑93 (MQ=255)
                                               tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                          <  2:1175307/93‑1 (MQ=255)
                                               tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                          >  2:2287865/1‑93 (MQ=255)
                                               tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                          >  2:2648645/1‑93 (MQ=255)
                                               tGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTc                                          <  1:894206/93‑1 (MQ=255)
                                                gACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCg                                         <  2:3464464/93‑1 (MQ=255)
                                                   cAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTcc                                      <  2:73550/93‑1 (MQ=255)
                                                   cAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTcc                                      >  1:1695498/1‑93 (MQ=255)
                                                      ttCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCa                                   >  2:2535801/1‑93 (MQ=255)
                                                       tCAGAGTTTCGCTAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaa                                  <  2:1372557/93‑1 (MQ=255)
                                                       tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaa                                  >  1:376626/1‑93 (MQ=255)
                                                       tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaa                                  >  1:1951693/1‑93 (MQ=255)
                                                       tCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAAGCCGCCTTTCGTCCTCaa                                  <  1:3241789/93‑1 (MQ=255)
                                                        cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                 <  1:2651659/93‑1 (MQ=255)
                                                        cAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCaaa                                 <  2:325787/93‑1 (MQ=255)
                                                         agagTTTTGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAAc                                <  1:4654880/93‑1 (MQ=255)
                                                         agagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAAc                                <  1:689722/93‑1 (MQ=255)
                                                         agagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAAc                                >  1:1491419/1‑93 (MQ=255)
                                                          gagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACa                               <  2:4347613/93‑1 (MQ=255)
                                                          gagTTTCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACa                               >  2:1673337/1‑93 (MQ=255)
                                                             tttCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATaa                            <  1:410884/93‑1 (MQ=255)
                                                             tttCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATaa                            <  2:1575061/93‑1 (MQ=255)
                                                               tCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                          <  1:2975539/93‑1 (MQ=255)
                                                               tCGCGAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAAtt                          <  1:4315427/93‑1 (MQ=255)
                                                                 gcgAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTg                        <  1:420456/93‑1 (MQ=255)
                                                                 gcgAGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTg                        >  1:4099727/1‑93 (MQ=255)
                                                                   gAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTa                      <  1:3571713/93‑1 (MQ=255)
                                                                    aGCGTATTGTTGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTAt                     <  1:228541/93‑1 (MQ=255)
                                                                    aGCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTAt                     >  1:565029/1‑93 (MQ=255)
                                                                     gCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                    <  1:3860634/93‑1 (MQ=255)
                                                                     gCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                    <  1:2944223/93‑1 (MQ=255)
                                                                     gCGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATg                    <  2:903635/93‑1 (MQ=255)
                                                                      cGTATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGa                   <  2:2393828/93‑1 (MQ=255)
                                                                        tATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGt                 >  2:3142710/1‑93 (MQ=255)
                                                                        tATTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGt                 <  1:1119131/93‑1 (MQ=255)
                                                                         aTTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                <  2:4662690/93‑1 (MQ=255)
                                                                         aTTGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTa                <  2:67923/93‑1 (MQ=255)
                                                                           tGTCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTAtt              <  1:4393407/93‑1 (MQ=255)
                                                                             tCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTaa            <  1:4215097/93‑1 (MQ=255)
                                                                             tCGAACTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTaa            <  1:2070698/93‑1 (MQ=255)
                                                                                aaCTGGGTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAAccc         >  2:1204199/1‑93 (MQ=255)
                                                                                  cTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATGTGTATGAGTATTAACCcgc       <  1:1306307/93‑1 (MQ=255)
                                                                                     ggTGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGAt    <  1:4081438/93‑1 (MQ=255)
                                                                                       tGCACCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTc  <  1:1982018/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTTGGCACCAACAACAAAAAAATCCGCCTTTCGTCCTCAAACATAATTTGTATGAGTATTAACCCGCGATTC  >  NC_000913/144849‑145027

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