Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,119 G→A 17.9% intergenic (+40/‑49) dnaK → / → dnaJ chaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,1190GA17.9% 444.2 / 77.3 195intergenic (+40/‑49)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base G (88/72);  new base A (18/17);  total (106/89)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

Only 100 of 205 total aligned reads displayed.

AACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAA  >  NC_000913/14020‑14218
                                                                                                   |                                                                                                   
aacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAgg                                                                                                     <  1:2447266/100‑1 (MQ=255)
    cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGaa                                                                                                 >  2:4231165/1‑100 (MQ=255)
    cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGaa                                                                                                 >  1:4361274/1‑100 (MQ=255)
       aaGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAttt                                                                                              >  2:94868/1‑100 (MQ=255)
         gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                            <  2:3676554/100‑1 (MQ=255)
         gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                            <  1:3609364/100‑1 (MQ=255)
         gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                            <  1:3160624/100‑1 (MQ=255)
           tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                            <  1:2712355/98‑1 (MQ=255)
           tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctc                                                                                          >  2:1887657/1‑100 (MQ=255)
             aCGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctct                                                                                         >  2:4056101/1‑99 (MQ=255)
             aCGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCtctc                                                                                        <  2:1017233/100‑1 (MQ=255)
              cGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTcc                                                                                       >  2:3795059/1‑100 (MQ=255)
                aTGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGc                                                                                     <  1:4384025/100‑1 (MQ=38)
                  gttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCACTCCGcc                                                                                    >  1:3991414/1‑99 (MQ=255)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                 <  2:3746913/100‑1 (MQ=37)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                 <  2:2896255/100‑1 (MQ=37)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                 <  2:1393417/100‑1 (MQ=37)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                 <  1:1074200/100‑1 (MQ=37)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGt                                                                                 >  1:1600179/1‑100 (MQ=255)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGt                                                                                 <  1:3583813/100‑1 (MQ=255)
                    tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACAGGCGAAGGGGAATTTCCTCTCAGCCCGt                                                                                 >  1:1883833/1‑100 (MQ=255)
                    ngtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGt                                                                                 >  1:1181412/2‑100 (MQ=255)
                     gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                  <  1:1532254/98‑1 (MQ=37)
                     gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                <  2:3257721/100‑1 (MQ=255)
                     gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                <  1:325368/100‑1 (MQ=255)
                     gtCGACGCTGAATTTAAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                 <  1:497730/99‑1 (MQ=25)
                      tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGc                                                                               >  2:229707/1‑100 (MQ=255)
                        gACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAGACGGGTAATTATACTGACACTGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                <  1:1662125/97‑1 (MQ=255)
                           gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                          <  2:4153548/100‑1 (MQ=255)
                           gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                          <  1:2991950/100‑1 (MQ=255)
                             tGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATc                                                                        <  2:678520/100‑1 (MQ=255)
                               aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                      >  2:3781312/1‑100 (MQ=255)
                                aTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAg                                                                     >  1:3101922/1‑100 (MQ=255)
                                  ttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                                                                   >  1:1935466/1‑100 (MQ=255)
                                    gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                                 >  2:2720107/1‑100 (MQ=255)
                                        aagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAtt                                                             <  1:1022951/100‑1 (MQ=255)
                                         agTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAttt                                                            >  1:2071457/1‑100 (MQ=255)
                                          gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa                                                           >  2:2438033/1‑100 (MQ=255)
                                          gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa                                                           >  2:1425498/1‑100 (MQ=255)
                                           tCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaa                                                          <  2:1843552/100‑1 (MQ=255)
                                            cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                         <  2:4003244/100‑1 (MQ=255)
                                            cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                         <  2:3334414/100‑1 (MQ=255)
                                            cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                         <  1:4108651/100‑1 (MQ=255)
                                            cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGAGAAATTCCTCTCCGCCCGTGCATTCCTCTAGGGGCAATTTaaa                                                         >  1:2820867/1‑100 (MQ=255)
                                              aaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaa                                                       <  1:615911/100‑1 (MQ=255)
                                                 ataaaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTACGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                     <  1:1436478/97‑1 (MQ=255)
                                                 aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                    >  2:2409712/1‑100 (MQ=255)
                                                 aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                    <  1:2091000/100‑1 (MQ=255)
                                                 aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                    <  2:2996696/100‑1 (MQ=255)
                                                  cAAAAACTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTACCACTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                   >  2:1225406/1‑100 (MQ=255)
                                                  cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGCAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                   >  1:265309/1‑100 (MQ=255)
                                                  cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                   >  1:3216666/1‑100 (MQ=255)
                                                  cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                   >  1:1346657/1‑100 (MQ=255)
                                                  cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                   >  2:312886/1‑100 (MQ=255)
                                                         taatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAg                                            <  1:3936202/100‑1 (MQ=255)
                                                            tCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCaa                                         <  2:3514365/100‑1 (MQ=255)
                                                             cGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                        >  1:2359327/1‑100 (MQ=255)
                                                             cGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                        >  1:384464/1‑100 (MQ=255)
                                                               cccTATAAACGGGTAATTATACTGACACGGACGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                      >  1:2891708/1‑100 (MQ=255)
                                                                ccTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatt                                     <  1:917901/100‑1 (MQ=255)
                                                                ccTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGa                                       >  1:1696070/1‑98 (MQ=255)
                                                                 cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCACTTTAAAAAAGATGGCTAAGCAAGatta                                    >  1:3791489/1‑100 (MQ=255)
                                                                 cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatta                                    <  2:3889684/100‑1 (MQ=255)
                                                                  tataAACGGGTAATTATACTGACAAGGGAAAAGGGGAATTTCGTGTCCGCCCGAGCAGTCATAAAGGTGAAAATTGAAAAAGCTAGCTaa                                             >  1:769127/1‑90 (MQ=255)
                                                                   ataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCACGattatt                                  >  2:876224/1‑100 (MQ=255)
                                                                    taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                 >  1:2188856/1‑100 (MQ=255)
                                                                      aaCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACg                               <  1:2206827/100‑1 (MQ=255)
                                                                       aCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGAAATTTAAAAAAGATGGCTAAGCGAGATTATTACga                              >  1:3192736/1‑100 (MQ=255)
                                                                          ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                           >  2:144403/1‑100 (MQ=255)
                                                                             aagtataCTGAGAGGGGCGAAGGGGAATGTCCTCTCCGGCGGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtttt                        <  1:3689837/97‑1 (MQ=255)
                                                                              aTTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGGTTATTACGAGATTTTa                       >  2:132021/1‑100 (MQ=255)
                                                                              aTTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                       <  1:48655/100‑1 (MQ=255)
                                                                                tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                     >  1:338582/1‑100 (MQ=255)
                                                                                tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                     <  2:1782988/100‑1 (MQ=255)
                                                                                   aCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGt                  >  2:35792/1‑100 (MQ=255)
                                                                                    cTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt                 >  2:536361/1‑100 (MQ=255)
                                                                                    cTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt                 >  2:3696958/1‑100 (MQ=255)
                                                                                    cTGACACGGGCGAAGGGGAATTTCCCCTCCGCCCGTGCATGCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCg                   >  1:1235785/1‑98 (MQ=255)
                                                                                        cacGGGCGGAGAGGAAGTTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCCAGATTATTACGAGATTTTAGGCGTTTCCa             >  1:1157379/1‑100 (MQ=38)
                                                                                        cacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa             >  1:1844400/1‑100 (MQ=37)
                                                                                        cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa             <  1:3719100/100‑1 (MQ=255)
                                                                                        cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa             <  2:4361305/100‑1 (MQ=255)
                                                                                        cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa             <  2:2671453/100‑1 (MQ=255)
                                                                                          cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGGTGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa           >  2:1168816/1‑100 (MQ=255)
                                                                                          cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa           <  2:75647/100‑1 (MQ=255)
                                                                                          cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa           <  1:2409701/100‑1 (MQ=255)
                                                                                           gggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa          >  2:2926349/1‑100 (MQ=255)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  1:3476891/1‑100 (MQ=38)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  1:4190745/1‑100 (MQ=38)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  1:2075629/1‑100 (MQ=38)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  2:1625029/1‑100 (MQ=38)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  2:2142612/1‑100 (MQ=38)
                                                                                            ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGGTTCCAAAAc         >  1:646776/1‑100 (MQ=38)
                                                                                            ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         <  1:1288958/100‑1 (MQ=255)
                                                                                            ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc         >  2:3536120/1‑100 (MQ=255)
                                                                                            ggCGAAGGAGAATTTCCTCTCCGCCCGTGCATTCATCTAGAGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa           <  1:3551193/98‑1 (MQ=255)
                                                                                             gCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACa        <  1:1193409/100‑1 (MQ=255)
                                                                                               gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGc      <  2:3967594/100‑1 (MQ=255)
                                                                                                 aGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg    <  2:2790764/100‑1 (MQ=255)
                                                                                                   gggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGaa  <  2:247790/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAA  >  NC_000913/14020‑14218

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: