Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,525,135 (A)6→5 coding (9/765 nt) yfeN → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,1300A.97.1% 744.7 / 12.0 173coding (4/765 nt)yfeNputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base . (83/85);  total (84/89)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 183 total aligned reads displayed.

GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTG  >  NC_000913/2525031‑2525225
                                                                                                   |                                                                                               
ggCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGa                                                                                                 >  2:948738/1‑100 (MQ=255)
  cGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaa                                                                                               <  1:134336/100‑1 (MQ=255)
  cGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaa                                                                                               <  2:1529261/100‑1 (MQ=255)
     tatTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaac                                                                                            <  1:2468046/100‑2 (MQ=255)
         aaTTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCt                                                                                       <  1:1305684/100‑1 (MQ=255)
         aaTTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCt                                                                                       >  1:1704614/1‑100 (MQ=255)
            tCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCtttt                                                                                    >  2:373236/1‑100 (MQ=255)
             ccTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCtt                                                                                      <  1:1577399/97‑1 (MQ=255)
              cTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTaa                                                                                  >  1:2671493/1‑100 (MQ=255)
              cTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTaa                                                                                  >  1:3225312/1‑100 (MQ=255)
               tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                                 >  1:2818507/1‑100 (MQ=255)
               tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                                 >  1:3006449/1‑100 (MQ=255)
               tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                                 <  2:2754140/100‑1 (MQ=255)
               tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                                 <  2:1949681/100‑1 (MQ=255)
               tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                                 <  2:1174753/100‑1 (MQ=255)
                  gCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActct                                                                              <  2:2259378/100‑1 (MQ=255)
                  gCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActct                                                                              >  2:2985488/1‑100 (MQ=255)
                   cATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTg                                                                             <  2:3434298/100‑1 (MQ=255)
                   cATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTg                                                                             >  2:1555915/1‑100 (MQ=255)
                   cATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTg                                                                             >  1:710112/1‑100 (MQ=255)
                    aTTAACAATATGTCGTTTTTAGCAAACCCCTTTATTAAGCAATCTACTAATAACGATCGGCATCGGGATTTATGCTATG‑AAAAACATCTTTTAACTCTg                                                                             >  2:3313951/1‑99 (MQ=255)
                     ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCACTCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGa                                                                            >  2:80529/1‑99 (MQ=255)
                     ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGac                                                                           >  2:1983170/1‑100 (MQ=255)
                     ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGac                                                                           <  2:2676811/100‑1 (MQ=255)
                      tAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGaca                                                                          <  1:2192558/100‑1 (MQ=255)
                      tAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGaca                                                                          >  2:1100302/1‑100 (MQ=255)
                       aaCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGacac                                                                         >  1:436457/1‑100 (MQ=255)
                       aaCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGacac                                                                         <  2:1854038/100‑1 (MQ=255)
                        aCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACt                                                                        >  1:2296572/1‑100 (MQ=255)
                        aCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACt                                                                        >  2:1369290/1‑100 (MQ=255)
                          aaTATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACttt                                                                      >  2:124170/1‑100 (MQ=255)
                            tatGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTcc                                                                    <  2:260082/100‑1 (MQ=255)
                             ntGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCct                                                                   >  1:3183641/2‑100 (MQ=255)
                              tGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCctc                                                                  >  2:2186848/1‑100 (MQ=255)
                              tGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAAAAAAAATA‑AAAAAAACATCTTTTAACTCTGACACTTTCctc                                                                  <  1:3387633/100‑1 (MQ=255)
                               gTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAAGTATGCTATG‑AAAAACATCTGTTAACTCGGACACTTTCctct                                                                 >  2:1895824/1‑100 (MQ=255)
                                tCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCta                                                                <  1:952643/100‑1 (MQ=255)
                                tCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCta                                                                <  2:3369685/100‑1 (MQ=255)
                                 ngTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtat                                                               >  1:3473521/2‑100 (MQ=255)
                                 cGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTccagatc                                                               >  1:1970506/1‑95 (MQ=255)
                                 cGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtat                                                               >  1:2608260/1‑100 (MQ=255)
                                 cGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtat                                                               >  1:2531169/1‑100 (MQ=255)
                                 cGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtat                                                               <  1:1916660/100‑1 (MQ=255)
                                  gTTTTTAGCAAGCCATTTTATTCAGCAATCTGCTATTAACGATCGGCATAGGGAATTATGCTATG‑ACAAACATCTTTTATCTCTGACACTTTCCTCtata                                                              <  2:1378195/100‑1 (MQ=255)
                                  gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtatn                                                              <  2:3473577/100‑2 (MQ=255)
                                  gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtata                                                              <  2:2496994/100‑1 (MQ=255)
                                     tttAGCAAACCACATTATTAAGCAATCTACTAATAACGATCGGCATAGGGAA‑AAAAAAAAAAAAAAACATCTTTTAACTCTGACACTTTCCTCTATAtt                                                            <  1:1397784/99‑1 (MQ=255)
                                       tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         >  1:2004963/1‑100 (MQ=255)
                                       tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         <  2:1054757/100‑1 (MQ=255)
                                       tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         <  2:620004/100‑1 (MQ=255)
                                       tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         <  1:1100288/100‑1 (MQ=255)
                                       ccgatctccACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         <  2:91223/93‑1 (MQ=255)
                                        aGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCg                                                        >  2:2174845/1‑100 (MQ=255)
                                         gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGTAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGa                                                       >  1:3236169/1‑100 (MQ=255)
                                         gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGa                                                       >  2:365955/1‑100 (MQ=255)
                                         gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGa                                                       >  1:947689/1‑100 (MQ=255)
                                         gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGa                                                       >  1:280101/1‑100 (MQ=255)
                                             aCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATAcc                                                   <  1:717729/100‑1 (MQ=255)
                                             aCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATAcc                                                   <  1:1216200/100‑1 (MQ=255)
                                              ccACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGcagatcgg                                                  >  1:91222/1‑93 (MQ=255)
                                              ccACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCt                                                  >  2:324881/1‑100 (MQ=255)
                                               cACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtg                                                 >  2:857046/1‑100 (MQ=255)
                                                aCTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtg                                                 >  1:1795858/1‑100 (MQ=255)
                                                 cTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgtt                                               >  2:2610131/1‑100 (MQ=255)
                                                   ttattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgttgt                                             >  1:2401917/1‑100 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            <  1:3377503/100‑1 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            >  1:2586109/1‑100 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            <  1:2862487/100‑1 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            >  1:2703448/1‑100 (MQ=255)
                                                      ttaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTTTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtc                                          <  1:1861707/100‑1 (MQ=255)
                                                        aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                        >  1:597301/1‑100 (MQ=255)
                                                        aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                        >  1:1260066/1‑100 (MQ=255)
                                                         aGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCa                                       <  2:784193/100‑1 (MQ=255)
                                                          gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCAt                                      <  1:96179/100‑1 (MQ=255)
                                                          gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCAt                                      <  1:1369271/100‑1 (MQ=255)
                                                          gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCAt                                      <  1:3006561/100‑1 (MQ=255)
                                                           cAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATg                                     <  2:472075/100‑1 (MQ=255)
                                                            aaTCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGc                                    <  1:1485146/100‑1 (MQ=255)
                                                                tactaATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGaa                                >  2:2836893/1‑100 (MQ=255)
                                                                tactaATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGaa                                >  1:3135509/1‑100 (MQ=255)
                                                                   taataaCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAttt                             <  2:3254170/100‑1 (MQ=255)
                                                                        aCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAgg                        <  1:2918764/100‑1 (MQ=255)
                                                                        aCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAgg                        <  1:2308247/100‑1 (MQ=255)
                                                                          gATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCg                      >  2:2055729/1‑100 (MQ=255)
                                                                           aTCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCgg                     <  2:593790/100‑1 (MQ=255)
                                                                            tCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGt                    <  1:104095/100‑1 (MQ=255)
                                                                                cATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTg                >  1:1532174/1‑100 (MQ=255)
                                                                                 aTAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc               <  1:399975/100‑1 (MQ=255)
                                                                                 aTAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc               <  1:2138488/100‑1 (MQ=255)
                                                                                 aTAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc               <  2:1704637/100‑1 (MQ=255)
                                                                                  tAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCa              <  2:2304209/100‑1 (MQ=255)
                                                                                    gggAATTATGCTATG‑AAATACATCTTTTAACTCTGACTCTTTCCTCTATATTAGCGATACCTGTTGTATCTTATGCAGAATTTAAAGGCGGTTTTGCata            <  1:1017350/100‑3 (MQ=255)
                                                                                          tATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCgg      <  1:1653079/100‑1 (MQ=255)
                                                                                          tATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCgg      <  2:3127705/100‑1 (MQ=255)
                                                                                           aTGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGc     >  2:1754057/1‑100 (MQ=255)
                                                                                            tGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCg    <  2:577752/100‑1 (MQ=255)
                                                                                            tGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCg    <  2:701595/100‑1 (MQ=255)
                                                                                            tGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCg    <  2:3006498/100‑1 (MQ=255)
                                                                                             gCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGt   >  1:124713/1‑100 (MQ=255)
                                                                                              cTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTg  >  1:834021/1‑100 (MQ=255)
                                                                                                   |                                                                                               
GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTG  >  NC_000913/2525031‑2525225

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