Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 106,554 A→G 55.3% intergenic (+98/‑3) ftsZ → / → lpxC GTP‑binding tubulin‑like cell division protein/UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913106,5540AG55.3% 74.7 / 196.5 163intergenic (+98/‑3)ftsZ/lpxCGTP‑binding tubulin‑like cell division protein/UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (46/44);  total (80/83)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

AAGAATTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTGCGCCGG  >  NC_000913/106455‑106653
                                                                                                   |                                                                                                   
 aGAATTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATAc                                                                                                    >  2:2140740/1‑100 (MQ=255)
 aGAATTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATAc                                                                                                    <  2:1480431/100‑1 (MQ=255)
   aaTTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCga                                                                                                  >  2:1950416/1‑100 (MQ=255)
    aTTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCgat                                                                                                 <  1:321289/100‑1 (MQ=255)
    aTTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGGACACTTCGGCTGGATAGGTAATTTGGAGAGATCATGCgat                                                                                                 >  2:3424852/1‑100 (MQ=255)
        accggcaTTTGGGTTTCGAGGCTCTTTGTGCTACACTGGCCCGCCGAATGTATAGTACACGTCTGTTGGATAGGTAATTTGGCGAGCTAACGCGATGATc                                                                                             >  2:207885/7‑100 (MQ=255)
        aCTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATc                                                                                             >  2:1986977/1‑100 (MQ=255)
         cTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATca                                                                                            >  1:2931822/1‑100 (MQ=255)
          tGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATcaa                                                                                           <  1:3168904/100‑1 (MQ=255)
              aTTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATcaaacaa                                                                                       <  2:442369/100‑1 (MQ=255)
                ttGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAg                                                                                     <  1:3797795/100‑1 (MQ=255)
                ttGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAg                                                                                     <  1:2860088/100‑1 (MQ=255)
                 tGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAgg                                                                                    <  2:1916874/100‑1 (MQ=255)
                 tGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAgg                                                                                    >  2:71570/1‑100 (MQ=255)
                 tGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAgg                                                                                    <  1:2079718/100‑1 (MQ=255)
                   ggTTTCGAGGCTGTTTGTGCTAAACTGGCACGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGac                                                                                  <  2:3551297/100‑1 (MQ=255)
                     tttCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGacac                                                                                <  1:2103071/100‑1 (MQ=255)
                      ttCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACt                                                                               >  1:4086239/1‑100 (MQ=255)
                      ttCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACt                                                                               <  1:1336526/100‑1 (MQ=255)
                           ggCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAAc                                                                          <  2:3316587/100‑1 (MQ=255)
                           ggCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAAc                                                                          >  1:963371/1‑100 (MQ=255)
                           ggCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAAc                                                                          <  1:1565253/100‑1 (MQ=255)
                           ggCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAAc                                                                          <  2:3735394/100‑1 (MQ=255)
                            gCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACg                                                                         <  2:2331546/100‑1 (MQ=255)
                            gCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACg                                                                         <  2:2305082/100‑1 (MQ=255)
                            gCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACg                                                                         >  2:4003568/1‑100 (MQ=255)
                                 ttGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCg                                                                    <  1:1462006/100‑1 (MQ=255)
                                    tgCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTc                                                                 >  1:1541407/1‑100 (MQ=255)
                                     gCTAAACTGGCCCGCCGAATGTATAGTACACTTCTGTTGGATAGGTAATTTGGCGAGATAAGGCGATGATCAAACAAAGGTCTCTTAAAAGTATCTTTCa                                                                >  2:3502913/1‑100 (MQ=255)
                                     gCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATGGTTCa                                                                >  2:727480/1‑100 (MQ=255)
                                     gCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCa                                                                >  2:858697/1‑100 (MQ=255)
                                      cTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAg                                                               <  1:2644378/100‑1 (MQ=255)
                                          aCTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAACCAAAGGACACTTAAACGTATCGTTCAGGCGa                                                           >  2:2097883/1‑100 (MQ=255)
                                            ttgtccgccgAATGTATAGTATACTTCGTTTGAATAGGTAATTTGTCAAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGcttc                                                          <  2:1839747/95‑4 (MQ=255)
                                             ggCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACgg                                                        <  2:1826126/100‑1 (MQ=255)
                                             ggCCCGCCGAATGTATAGTACACTTCGGATGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGAc                                                          >  2:1354859/1‑98 (MQ=255)
                                              gCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACggg                                                       >  1:330242/1‑100 (MQ=255)
                                              gCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACggg                                                       <  1:1577134/100‑1 (MQ=255)
                                              gCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACggg                                                       >  2:1616723/1‑100 (MQ=255)
                                              gCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACggg                                                       <  1:556409/100‑1 (MQ=255)
                                              gCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACgcg                                                       <  1:2792348/100‑3 (MQ=255)
                                                ccgccgAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGgtg                                                     >  1:1046178/1‑100 (MQ=255)
                                                ccgccgAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGgtg                                                     >  1:831283/1‑100 (MQ=255)
                                                    cgAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCgg                                                 <  2:1270850/100‑1 (MQ=255)
                                                      aaTGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGtt                                               <  1:1365005/100‑1 (MQ=255)
                                                       aTGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGttt                                              <  1:2126098/100‑1 (MQ=255)
                                                        ttccgatctcacTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTa                                             <  2:2417545/91‑1 (MQ=255)
                                                        tGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATGGGGCAGGCGACGGGTGTCGGTTTa                                             <  1:52751/100‑1 (MQ=255)
                                                        tGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAAAGTATCGTTCAGGCGACGGGTGTAGGTTTa                                             >  2:1026373/1‑100 (MQ=255)
                                                           ataGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTAc                                            <  2:3985282/98‑1 (MQ=255)
                                                           ataGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACAt                                          >  2:3954936/1‑100 (MQ=255)
                                                           ataGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACAt                                          <  1:1811408/100‑1 (MQ=255)
                                                             aGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATAc                                        <  1:3035723/100‑1 (MQ=255)
                                                              gTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATAcc                                       >  1:104259/1‑100 (MQ=255)
                                                              gTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATAcc                                       <  2:929293/100‑1 (MQ=255)
                                                                acacTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCgn                                     <  1:3594614/100‑2 (MQ=255)
                                                                acacTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCgg                                     >  2:2101511/1‑100 (MQ=255)
                                                                acacTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCg                                      >  2:1729202/1‑99 (MQ=255)
                                                                 cacTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTaagatcggaa                                    >  1:2417559/1‑91 (MQ=255)
                                                                  acTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCa                                   <  1:1213621/100‑1 (MQ=255)
                                                                   cTTCGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCaa                                  >  1:1681688/1‑100 (MQ=255)
                                                                      cGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGaa                               <  1:4040363/100‑1 (MQ=255)
                                                                      cGGTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGaa                               >  1:3048829/1‑100 (MQ=255)
                                                                       ggTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGaaa                              >  2:3893266/1‑100 (MQ=255)
                                                                        gTTGGATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAg                             <  1:3085384/100‑1 (MQ=255)
                                                                        gTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGaaat                             >  1:1924556/1‑99 (MQ=255)
                                                                          ttgATAGGTAAATTTGCGAGATTATACGATTATCAAACAAAGGATACTTAAACTTATCGTTCAGGCGACGTGTGTCGGTTTACATACCGGCAAGa                                <  2:3278926/93‑1 (MQ=255)
                                                                          tGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTc                           <  1:2558054/100‑1 (MQ=255)
                                                                           ggATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCa                          >  2:3618804/1‑100 (MQ=255)
                                                                            gATAGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCAc                         >  1:631410/1‑100 (MQ=255)
                                                                            gATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCAc                         <  1:1895011/100‑1 (MQ=255)
                                                                               aGGTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCt                      <  2:2968145/100‑1 (MQ=255)
                                                                                ggTAATTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTg                     >  1:4125306/1‑100 (MQ=255)
                                                                                ggTAATTTGGCGAGATAATACGCTGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTg                     >  2:1541410/1‑100 (MQ=255)
                                                                                ggTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTg                     >  1:2475069/1‑100 (MQ=255)
                                                                                ggTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGGCACCCTg                     >  2:407594/1‑100 (MQ=255)
                                                                                 gTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACTGGCAAGAAAGTCACCCTGa                    >  2:4133442/1‑100 (MQ=255)
                                                                                  tAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGAc                   <  2:49010/100‑1 (MQ=255)
                                                                                   aaTTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACg                  >  1:341201/1‑100 (MQ=255)
                                                                                   aaTTTGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACg                  <  1:1190153/100‑1 (MQ=255)
                                                                                     tttGGCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGtt                >  2:1637649/1‑100 (MQ=255)
                                                                                       tGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTAc              >  2:109619/1‑100 (MQ=255)
                                                                                        ggCGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACg             >  2:2211064/1‑100 (MQ=255)
                                                                                        ggCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACg             <  1:3365936/100‑1 (MQ=255)
                                                                                        ggCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACg             >  2:643165/1‑100 (MQ=255)
                                                                                        ggCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACg             >  2:2779508/1‑100 (MQ=255)
                                                                                         gCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTAc              >  2:202114/1‑98 (MQ=255)
                                                                                         gCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGc            <  1:3203179/100‑1 (MQ=255)
                                                                                          cGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGcc           >  2:3774704/1‑100 (MQ=255)
                                                                                          cGAGATAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGcc           >  1:3984524/1‑100 (MQ=255)
                                                                                          cGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGcc           >  2:1879590/1‑100 (MQ=255)
                                                                                           gagaTCATACGATGCTCAAACAAAGGACACTCAACCGTATCCTTCAGGCGACGCGCGCCGCTTTCCACACCGGCAAGACCCCCCCCCTGACCCTACGccc          >  2:3203716/1‑100 (MQ=255)
                                                                                           gagaTAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGccc          <  2:143676/100‑1 (MQ=255)
                                                                                           gagaTAATACGACGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGATTTACATACCGGCAAGAATGTCACCCTGACGTTACGccc          <  2:2530969/100‑1 (MQ=255)
                                                                                            agaTAATGCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCt         >  1:3122632/1‑100 (MQ=255)
                                                                                             gaTAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTg        >  2:3382725/1‑100 (MQ=255)
                                                                                               taatgCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTgcg      >  1:1906149/1‑100 (MQ=255)
                                                                                               taatgCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTgcg      <  1:2321579/100‑1 (MQ=255)
                                                                                               taataCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTgcg      >  2:3738813/1‑100 (MQ=255)
                                                                                                 ataCGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTGCGcc    >  2:3075254/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAGAATTGACTGGAATTTGGGTTTCGAGGCTCTTTGTGCTAAACTGGCCCGCCGAATGTATAGTACACTTCGGTTGGATAGGTAATTTGGCGAGATAATACGATGATCAAACAAAGGACACTTAAACGTATCGTTCAGGCGACGGGTGTCGGTTTACATACCGGCAAGAAAGTCACCCTGACGTTACGCCCTGCGCCGG  >  NC_000913/106455‑106653

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