Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,007 C→T 17.9% G21G (GGC→GGT 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,0070CT17.9% 233.7 / 39.2 111G21G (GGC→GGTtatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base C (51/38);  new base T (7/12);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 116 total aligned reads displayed.

GCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGT  >  NC_000913/4021908‑4022106
                                                                                                   |                                                                                                   
gCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGc                                                                                                     <  2:1885026/100‑1 (MQ=255)
gCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCGTTTTGGc                                                                                                     >  2:2622621/1‑100 (MQ=255)
gCTGTGTTGAGTCATCATCTGCCAGAGAGGAACATGTATGGGTGGTATCAGTATTTGGGAGTTAGTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGc                                                                                                     <  1:34654/100‑1 (MQ=255)
    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:2224155/100‑1 (MQ=255)
          atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                           <  2:892799/100‑1 (MQ=255)
          atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGc                                                                                           >  1:1301/1‑100 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:255392/100‑2 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:2660757/100‑2 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCt                                                                                          >  2:177844/1‑100 (MQ=255)
            catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                         <  1:2585548/100‑3 (MQ=255)
            aatcatcTACCACAGAGGAACATGGATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGGTGTACTGCTCTTTGGTGCCAAAAAGcag                                                                                         <  1:576964/99‑3 (MQ=255)
             atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                        <  2:1500658/100‑1 (MQ=255)
             atcatcTACCACAGAGGAACATGTATGGGTGGTAGCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                        <  2:170460/100‑1 (MQ=255)
              tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  2:2543069/100‑5 (MQ=255)
              tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  2:435082/100‑5 (MQ=255)
               catcTACATCAGAGGAACATGTATGGGTGGCATCAGGATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtac                                                                                      <  2:1829466/100‑6 (MQ=255)
                atcTACCACAGAGGAACATGTATGGGCGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCt                                                                                     >  2:1772578/1‑100 (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:1243678/100‑1 (MQ=255)
                     ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  2:1314666/1‑100 (MQ=255)
                     ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  2:900875/1‑100 (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)
                          gagGAACATTTATGGGTCGTAGCAGTATTTGTACGTTATTGATTAGTGCCGTCATCGTTGTACTGCTTTTTTGCACTAAAAAGCTCGGCTCCATCg                                                                               <  1:2066085/96‑1 (MQ=255)
                          gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTc                                                                           <  2:1522677/100‑1 (MQ=255)
                           agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:435849/100‑1 (MQ=255)
                           agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:457659/100‑1 (MQ=255)
                           agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:992067/100‑1 (MQ=255)
                              aaCATGGATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAATGCTCGTCTCCATCGGTTCCGAt                                                                       >  1:2621085/1‑100 (MQ=255)
                               aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                      <  2:2786381/100‑1 (MQ=255)
                                cATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCTGTTCCGATc                                                                      <  1:277022/99‑1 (MQ=255)
                                   gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  1:2107064/100‑1 (MQ=255)
                                   gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  1:2106573/100‑1 (MQ=255)
                                     aTGGGGGGTATCAGTGTTTGGGAGTTATTGATTATTGCCGTCATCGGTGTACTGCTTTTTGGGACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTg                                                                <  1:848807/100‑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)
                                         gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTc                                                            <  2:2791107/100‑1 (MQ=255)
                                          tggtATCAGTATTTGGCAGTTATTGATTATTTCTGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCg                                                           <  1:2296380/100‑1 (MQ=255)
                                            gtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                          >  2:1315555/1‑99 (MQ=255)
                                             tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        >  2:2740911/1‑100 (MQ=255)
                                             tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        >  2:2140189/1‑100 (MQ=255)
                                             tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        <  2:456920/100‑1 (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                                             >  2:1864428/1‑100 (MQ=255)
                                                        ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  1:1542960/1‑100 (MQ=255)
                                                        ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaa                                             >  2:805919/1‑100 (MQ=255)
                                                             ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                        <  2:2622917/100‑1 (MQ=255)
                                                                ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCaa                                     >  1:790850/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)
                                                                        tGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCAGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGc                                <  1:2366393/97‑1 (MQ=255)
                                                                         gCCGTCATCGTTGTACTGCTTTTTTTTTTTTTTTAGCTCGGCTCCATCGGTTCCGATCTTTGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatg                            >  1:2629263/1‑100 (MQ=255)
                                                                          ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                           >  1:2269831/1‑100 (MQ=255)
                                                                              cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATTTTGGTGCGTAGATCAAAGGCCATAAAAAAGCAATGAACGATGATGaa                       >  1:2273745/1‑100 (MQ=255)
                                                                              cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTCAAAATAGCAATGAGCGAtgatga                        >  1:774155/1‑99 (MQ=255)
                                                                               ctCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAc                      <  2:2823177/99‑1 (MQ=255)
                                                                                 cGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaaaag                    >  1:2823159/1‑97 (MQ=255)
                                                                                   ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCATCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCaa                   >  1:727214/1‑99 (MQ=255)
                                                                                     gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                >  1:2643800/1‑100 (MQ=255)
                                                                                      ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  2:815431/5‑100 (MQ=255)
                                                                                      ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  1:2329875/5‑100 (MQ=255)
                                                                                      ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  1:2413738/5‑100 (MQ=255)
                                                                                       cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTTGTTCAAAGGCTTTCAAACAGCAATGAGCGATGATCAACCAAAGCa               >  1:1201999/4‑99 (MQ=255)
                                                                                       cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg              >  2:1333281/4‑100 (MQ=255)
                                                                                        gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:2724556/3‑100 (MQ=255)
                                                                                        cttctgtttGGCACTAATAGGCGCGGCTACATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             <  1:2456143/94‑1 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:1837111/1‑100 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  2:676950/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: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)
                                                                                                   caccacAAAGGTCGGATCAAGCGGTTCTGCTCTTGGAGCGGGCAACACAGGCTTTAAAAAAGCAAGGAGGCATGATGAGCCAAAGCAGGCTaaaaac     >  1:2816378/1‑95 (MQ=255)
                                                                                                   cACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt  <  2:2457548/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGT  >  NC_000913/4021908‑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: