Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,008 A→G 15.1% T22A (ACC→GCC)  tatA → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,022,0080AG15.1% 316.3 / 39.9 135T22A (ACC→GCC) tatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (61/53);  new base G (3/18);  total (64/71)
Fisher's exact test for biased strand distribution p-value = 8.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 138 total aligned reads displayed.

CTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTC  >  NC_000913/4021909‑4022107
                                                                                                   |                                                                                                   
cTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                     <  2:2408526/100‑1 (MQ=255)
  ttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAcc                                                                                                   <  1:2734015/100‑1 (MQ=255)
  ttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAcc                                                                                                   <  2:2736333/100‑1 (MQ=255)
     tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTTCTGCTTTTTGGCACCaaa                                                                                                >  1:684193/1‑100 (MQ=255)
         atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGc                                                                                            <  1:1128729/100‑1 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:2274728/100‑2 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:2593230/100‑2 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCt                                                                                           >  1:1407515/1‑100 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCt                                                                                           >  2:513245/1‑100 (MQ=255)
           catcatcTCCCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  1:1862295/100‑3 (MQ=255)
           catcatcTACCACAGAGGATCATGTATGTCTCGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTc                                                                                          <  1:592964/100‑1 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  2:2771939/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  1:550177/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  1:1800663/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  1:1322140/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTc                                                                                          <  2:779778/100‑1 (MQ=255)
            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagt                                                                                         <  2:2252154/100‑4 (MQ=255)
            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagt                                                                                         <  2:975453/100‑4 (MQ=255)
            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagt                                                                                         <  1:2126087/100‑4 (MQ=255)
            atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagt                                                                                         <  2:154166/100‑4 (MQ=255)
             tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                        <  2:1396415/100‑5 (MQ=255)
             tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCgg                                                                                        >  1:835945/1‑100 (MQ=255)
             tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCgg                                                                                        <  2:1938245/100‑1 (MQ=255)
                tcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTc                                                                                     <  1:1397959/100‑1 (MQ=255)
                 cTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtacaac                                                                                    <  2:362746/100‑9 (MQ=38)
                  tACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCa                                                                                   >  2:1286762/1‑100 (MQ=255)
                    ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                 >  1:1315833/1‑100 (MQ=255)
                    ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                 >  2:1933058/1‑100 (MQ=255)
                      acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                               >  2:2501747/1‑100 (MQ=255)
                       catagtaaCATGTATGCGTGGTATCAGTATTTGGCAGTGATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGt                                                                              <  2:2091537/94‑1 (MQ=255)
                          agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                           <  2:1880688/100‑1 (MQ=255)
                            gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGa                                                                         <  1:1552136/100‑1 (MQ=255)
                              aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                       >  2:402097/1‑100 (MQ=255)
                              aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                       <  1:2035157/100‑1 (MQ=255)
                              aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                       <  2:1407521/100‑1 (MQ=255)
                              aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                       >  1:2082188/1‑100 (MQ=255)
                               cATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                      >  1:1050679/1‑100 (MQ=255)
                               cATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                      >  1:2703687/1‑100 (MQ=255)
                                aTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                     <  2:2084033/100‑1 (MQ=255)
                                aTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                      >  2:2440730/1‑99 (MQ=255)
                                 tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                    >  1:532181/1‑100 (MQ=255)
                                  gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                   <  2:2536980/100‑1 (MQ=255)
                                  gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                   <  2:1950086/100‑1 (MQ=255)
                                     tGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGc                                                                >  1:698407/1‑100 (MQ=255)
                                      gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTGTGGCGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCg                                                               <  1:1219490/100‑1 (MQ=255)
                                       ggtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCg                                                               >  1:2878946/1‑99 (MQ=255)
                                       ggtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                              >  1:1232385/1‑100 (MQ=255)
                                        gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTc                                                             <  1:2507424/100‑1 (MQ=255)
                                        gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTc                                                             >  1:1592338/1‑100 (MQ=255)
                                          ggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGTCACCAAAAAGCTCGTCTCCATCGGTTCCGATCTTTGttc                                                                >  1:626918/1‑93 (MQ=255)
                                            tATCAGTTTTTGGTAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCGCCAAAAAGCTCGGCTCCCTCGGTTCCGATCTTGGTGCGTCGATc                                                         <  1:2899982/100‑1 (MQ=255)
                                            tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                         <  2:2472717/100‑1 (MQ=255)
                                            tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAACTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                         <  2:2150307/100‑1 (MQ=255)
                                             atgaGTATTTGGCTGTTGTTAATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGGGCGTCGATCa                                                        <  1:981070/97‑1 (MQ=255)
                                             aTCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCa                                                        >  1:2811021/1‑100 (MQ=255)
                                             aTCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCa                                                        <  2:2163125/100‑1 (MQ=255)
                                               cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                      >  2:1763791/1‑100 (MQ=255)
                                                aGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAg                                                     >  2:451542/1‑100 (MQ=255)
                                                 gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAgg                                                    >  2:119600/1‑100 (MQ=255)
                                                    tttGGCAGTTAGTGATTATTGGCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCtt                                                 <  1:1491442/100‑1 (MQ=255)
                                                      tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGTTGCGTCGATCAAAGGCTTTa                                               <  2:2347457/100‑1 (MQ=255)
                                                      tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTa                                               >  2:1135528/1‑100 (MQ=255)
                                                           gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                          <  1:2243845/100‑1 (MQ=255)
                                                            ttATTGATTATTGCCGTCATCGTTGTACTGTTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                         <  1:1968482/100‑1 (MQ=255)
                                                            ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                         <  2:2682487/100‑1 (MQ=255)
                                                             tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGc                                        >  2:731314/1‑100 (MQ=255)
                                                               ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCaa                                      >  1:2185360/1‑100 (MQ=255)
                                                                tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                     <  2:2862730/100‑1 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    >  1:80606/1‑100 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    >  2:828834/1‑100 (MQ=255)
                                                                  attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                   >  1:2839105/1‑100 (MQ=255)
                                                                    tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGc                                 <  2:2545082/100‑1 (MQ=255)
                                                                         ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                            <  2:3037015/100‑1 (MQ=255)
                                                                          cGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatgat                           <  2:574491/100‑1 (MQ=255)
                                                                              aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAc                       <  2:698412/100‑1 (MQ=255)
                                                                                cGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCa                     <  2:2292509/100‑1 (MQ=255)
                                                                                cGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCa                     >  2:1476788/1‑100 (MQ=255)
                                                                                      cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               >  2:3025976/4‑100 (MQ=255)
                                                                                      aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               >  1:2393057/1‑100 (MQ=255)
                                                                                      aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               <  1:1984012/100‑1 (MQ=255)
                                                                                      aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               >  2:469781/1‑100 (MQ=255)
                                                                                      aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               >  2:1866994/1‑100 (MQ=255)
                                                                                       gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  1:943268/3‑100 (MQ=255)
                                                                                       cTGCTTTTTGGCACCGAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  1:1699940/1‑100 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              <  1:329556/100‑1 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              <  2:1050685/100‑1 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  2:193550/1‑100 (MQ=255)
                                                                                        tGCTTTTTGGCACCAAAAAGCTCGGCTGCATCGGTTCCGATCTTGGTGCGTCGAACAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               <  1:1086898/98‑1 (MQ=255)
                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt            >  2:608154/1‑100 (MQ=255)
                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt            <  1:2245968/100‑1 (MQ=255)
                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCAGTACCGATCTTGGTGAGTCGCTCAAAGGAAATAACAAAGAAATGAACAATGATGAACCAAAg                  >  1:2775969/1‑94 (MQ=255)
                                                                                          ttcgtaggCGCCAAAAAGCTAGGCTCCATCAGTTCCAATCTTGTTGCGCCGACCAAAGGCTTTAAAAAAGCAATGAGCGAGGAGGAACCTAAGCAGGac            <  1:335099/93‑2 (MQ=255)
                                                                                          cTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATa           >  2:902692/1‑100 (MQ=255)
                                                                                               tGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc      <  2:551771/100‑1 (MQ=255)
                                                                                               tGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc      <  1:1001457/100‑1 (MQ=255)
                                                                                                ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa     >  2:2417385/1‑100 (MQ=255)
                                                                                                ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa     >  2:2600182/1‑100 (MQ=255)
                                                                                                  cACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt   >  1:2939751/1‑100 (MQ=255)
                                                                                                   aCCAAAAAGCTCGGCTCGATCGGTTCCGATCTTGGTGCGTTGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc  <  1:2453711/100‑1 (MQ=255)
                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc  <  2:351155/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTC  >  NC_000913/4021909‑4022107

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: