Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,008 A→G 16.3% 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,0080AG16.3% 235.4 / 38.7 112T22A (ACC→GCC) tatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (51/41);  new base G (6/12);  total (58/54)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 113 total aligned reads displayed.

TGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGT  >  NC_000913/4021912‑4022106
                                                                                                |                                                                                                  
tGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCa                                                                                                 <  1:511289/100‑1 (MQ=255)
 gTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGTCACCaa                                                                                                >  2:1720122/1‑100 (MQ=255)
 gTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaa                                                                                                >  2:811697/1‑100 (MQ=255)
  tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaa                                                                                               >  2:2353384/1‑100 (MQ=255)
   ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaa                                                                                              >  1:1797275/1‑100 (MQ=255)
   ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaa                                                                                              >  2:2462362/1‑100 (MQ=255)
      atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                           <  2:892799/100‑1 (MQ=255)
      atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGc                                                                                           >  1:1301/1‑100 (MQ=255)
       tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:2660757/100‑2 (MQ=255)
        catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                         <  1:2585548/100‑3 (MQ=255)
         atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                        <  2:1500658/100‑1 (MQ=255)
         atcatcTACCACAGAGGAACATGTATGGGTGGTAGCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                        <  2:170460/100‑1 (MQ=255)
          tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  2:435082/100‑5 (MQ=255)
          tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  2:2543069/100‑5 (MQ=255)
           catcTACATCAGAGGAACATGTATGGGTGGCATCAGGATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtac                                                                                      <  2:1829466/100‑6 (MQ=255)
            atcTACCACAGAGGAACATGTATGGGCGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCt                                                                                     >  2:1772578/1‑100 (MQ=255)
              cttccACAGAGGAAGATGTATGGCTGGTATTAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTATACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  1:2017463/93‑5 (MQ=255)
              cTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtacaac                                                                                   <  1:2536122/100‑9 (MQ=38)
              cTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                   >  2:627043/1‑100 (MQ=255)
               tACCACAAAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                   <  1:627035/99‑1 (MQ=255)
                aCCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCAt                                                                                 >  1:175634/1‑100 (MQ=255)
                 ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  2:900875/1‑100 (MQ=255)
                 ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  2:1314666/1‑100 (MQ=255)
                 ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                <  2:1243678/100‑1 (MQ=255)
                  cacaGAGGAACATGTATGGGTGGTATCCGTATTTGGCAGTTATTGAGTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCGCGGGTCCATCg                                                                               <  1:2714369/100‑1 (MQ=255)
                  cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCg                                                                               >  1:763666/1‑100 (MQ=255)
                     agagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCTTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGtt                                                                            >  2:2272815/1‑100 (MQ=255)
                     agagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGtt                                                                            <  2:473876/100‑1 (MQ=255)
                      gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTc                                                                           <  2:1522677/100‑1 (MQ=255)
                       agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:992067/100‑1 (MQ=255)
                       agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:435849/100‑1 (MQ=255)
                       agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:457659/100‑1 (MQ=255)
                           aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                      <  2:2786381/100‑1 (MQ=255)
                            cATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCTGTTCCGATc                                                                      <  1:277022/99‑1 (MQ=255)
                              tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                   <  2:1037994/100‑1 (MQ=255)
                               gcaTGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  2:1501814/98‑1 (MQ=255)
                               gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  1:2107064/100‑1 (MQ=255)
                               gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  1:2106573/100‑1 (MQ=255)
                                tATGGGAGGTTTCAGTATTTGGCAGTTATTGATTATGGCCGTCATCGTTGTACTGCTTTATAAAAAAAAAAAGCTCGGCTCCATCGGTTCCGAt                                                                       <  1:263781/94‑1 (MQ=255)
                                  tGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGc                                                               >  1:812901/1‑100 (MQ=255)
                                   gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCg                                                              <  1:2848110/100‑1 (MQ=255)
                                    ggtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                             >  2:2073906/1‑100 (MQ=255)
                                     gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTc                                                            <  2:1840232/100‑1 (MQ=255)
                                      tgttATCAGTATTTGGTAGTGATTGATTATTGCTGTCATCGTTGTACTTCTTTTTGGCACCAAAAAGCTCGGCTCCATCCGTTCCGATCTTGGTGCGTCg                                                           <  1:540439/97‑1 (MQ=255)
                                      tggtATCAGTATTTGGCAGTTATTGATTATTTCTGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCg                                                           <  1:2296380/100‑1 (MQ=255)
                                        gtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                          >  2:1315555/1‑99 (MQ=255)
                                         tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        <  2:456920/100‑1 (MQ=255)
                                         tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        >  2:2140189/1‑100 (MQ=255)
                                         tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        >  2:2740911/1‑100 (MQ=255)
                                          aTCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCa                                                       >  2:275544/1‑100 (MQ=255)
                                           tgaGTATTTGGCATTTATTGAGTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        <  1:2285073/96‑1 (MQ=255)
                                           tCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaa                                                      >  1:863272/1‑100 (MQ=255)
                                            cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                     <  2:40651/100‑1 (MQ=255)
                                            cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                     >  2:821905/1‑100 (MQ=255)
                                                agttGGCAGTTATGGATGATGGCCGTCATCGTTGTACTGCTTGTTGGCCCCAAAGAGCTCGGCTCCATCGGGTCCGATCTTGGTGCGTCGATCAAAGGc                                                  <  1:724383/97‑1 (MQ=255)
                                                    ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  1:1542960/1‑100 (MQ=255)
                                                    ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  2:805919/1‑100 (MQ=255)
                                                    ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  2:780871/1‑100 (MQ=255)
                                                    ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  2:1864428/1‑100 (MQ=255)
                                                         ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                        <  2:2622917/100‑1 (MQ=255)
                                                            ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCaa                                     >  1:790850/1‑100 (MQ=255)
                                                             tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                    >  1:2508921/1‑100 (MQ=255)
                                                             tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                    >  2:2900942/1‑100 (MQ=255)
                                                             tGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACAAAAAAGCTCGGCTCAATCGGATCCGATCTTGGTGCGTCGATCAAAGGCTGTAAAAAAGCaa                                     >  1:2614562/1‑99 (MQ=255)
                                                             tGATTATGGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                    <  2:2703260/100‑1 (MQ=255)
                                                              gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                   >  1:2405240/1‑100 (MQ=255)
                                                                 tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGc                                <  1:275543/100‑1 (MQ=255)
                                                                     gCCGTCATCGTTGTACTGCTTTTTTTTTTTTTTTAGCTCGGCTCCATCGGTTCCGATCTTTGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatg                            >  1:2629263/1‑100 (MQ=255)
                                                                      ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                           >  1:2269831/1‑100 (MQ=255)
                                                                      ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                           >  1:2292750/1‑100 (MQ=255)
                                                                      ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatg                            <  2:2729411/99‑1 (MQ=255)
                                                                        gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGatgatg                         <  2:2187099/100‑1 (MQ=255)
                                                                          cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATTTTGGTGCGTAGATCAAAGGCCATAAAAAAGCAATGAACGATGATGaa                       >  1:2273745/1‑100 (MQ=255)
                                                                          cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTCAAAATAGCAATGAGCGAtgatga                        >  1:774155/1‑99 (MQ=255)
                                                                           ctCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAc                      <  2:2823177/99‑1 (MQ=255)
                                                                            tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAcc                     >  2:1835165/1‑100 (MQ=255)
                                                                            tCGTTGTACTGCTTTTTGGCACCAAAAACCCCGGCTCCATCGGGTCCGATCTTGGTGCATCGACCAAAGACTACCAAAAAACAATGAGCcaggataacc                      >  1:283065/1‑95 (MQ=255)
                                                                             cGTTGTACTGCTTTTTGGGACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAcc                     <  1:775228/99‑1 (MQ=255)
                                                                             cGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaaaag                    >  1:2823159/1‑97 (MQ=255)
                                                                             cGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCa                    >  2:2513506/1‑100 (MQ=255)
                                                                               ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCATCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCaa                   >  1:727214/1‑99 (MQ=255)
                                                                                 gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                >  1:2643800/1‑100 (MQ=255)
                                                                                  ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  1:2329875/5‑100 (MQ=255)
                                                                                  ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  1:2413738/5‑100 (MQ=255)
                                                                                  ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  2:815431/5‑100 (MQ=255)
                                                                                   cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTTGTTCAAAGGCTTTCAAACAGCAATGAGCGATGATCAACCAAAGCa               >  1:1201999/4‑99 (MQ=255)
                                                                                    gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:2724556/3‑100 (MQ=255)
                                                                                    cttctgtttGGCACTAATAGGCGCGGCTACATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             <  1:2456143/94‑1 (MQ=255)
                                                                                    cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:1837111/1‑100 (MQ=255)
                                                                                    cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGACGAACCAAAGCAgg             >  2:1019290/1‑100 (MQ=255)
                                                                                      gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt           >  2:1415586/1‑100 (MQ=255)
                                                                                        tttttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaa         >  1:2584023/1‑100 (MQ=255)
                                                                                           ttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAc      <  2:210402/100‑1 (MQ=255)
                                                                                            tGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc     <  2:322343/100‑1 (MQ=255)
                                                                                            tGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAcc     <  2:2677789/100‑1 (MQ=255)
                                                                                             ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa    >  2:2225300/1‑100 (MQ=255)
                                                                                             ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa    >  1:420598/1‑100 (MQ=255)
                                                                                              gCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAg   >  2:1391943/1‑100 (MQ=255)
                                                                                              gCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAg   >  1:1718709/1‑100 (MQ=255)
                                                                                               cACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt  <  2:2457548/100‑1 (MQ=255)
                                                                                                |                                                                                                  
TGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGT  >  NC_000913/4021912‑4022106

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