Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,008 A→G 15.7% 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.7% 422.3 / 71.2 199T22A (ACC→GCC) tatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (94/73);  new base G (12/19);  total (107/92)
Fisher's exact test for biased strand distribution p-value = 8.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-02

Only 100 of 206 total aligned reads displayed.

CTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTC  >  NC_000913/4021909‑4022107
                                                                                                   |                                                                                                   
cTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                     <  1:2422235/100‑1 (MQ=255)
cTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                     <  2:4075902/100‑1 (MQ=255)
  ttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAcc                                                                                                   <  1:2021303/100‑1 (MQ=255)
   tGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCa                                                                                                  >  2:1407788/1‑100 (MQ=255)
     tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaa                                                                                                >  2:1010332/1‑100 (MQ=255)
       tAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaaa                                                                                              >  2:3672147/1‑100 (MQ=255)
        aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAg                                                                                             >  1:88910/1‑100 (MQ=255)
         atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGc                                                                                            >  2:2190306/1‑100 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:3331933/100‑2 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:1090722/100‑2 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:4380481/100‑2 (MQ=255)
          tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                           <  2:193743/100‑2 (MQ=255)
          tcatcatcTACAACAAAGGAACATGTATGGGTGATATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCt                                                                                           >  1:2108280/1‑100 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGTCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  1:2482065/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  2:3377536/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                          <  2:727475/100‑3 (MQ=255)
           catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTc                                                                                          >  2:2020147/1‑100 (MQ=255)
             tcatcTTCCGCAGTGGAACATGTATGGTTGGGATCAGGATTTGTCAGTTATTGATTATTGCCGTTATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                        <  2:2562794/100‑5 (MQ=255)
              catcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGc                                                                                       <  1:1010328/100‑1 (MQ=255)
               atcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGc                                                                                       <  2:1989404/99‑1 (MQ=255)
               atcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGGTGTACTGCTTTTTGGCACCAAAAAGCTCGGCt                                                                                      <  2:2617599/100‑1 (MQ=255)
                  tACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtacaacg                                                                                   <  2:3064195/100‑10 (MQ=38)
                  tACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCa                                                                                   <  2:3410742/100‑1 (MQ=255)
                    ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                 >  1:802100/1‑100 (MQ=255)
                     cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCg                                                                                >  2:951408/1‑100 (MQ=255)
                     cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTAGCACCAAGAAGCTCGGCTCCATc                                                                                 >  2:4007128/1‑99 (MQ=255)
                      acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                               >  1:1362355/1‑100 (MQ=255)
                      acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                               >  2:582766/1‑100 (MQ=255)
                         gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGtt                                                                             <  1:1340415/99‑1 (MQ=255)
                         gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTc                                                                            >  2:3575169/1‑100 (MQ=255)
                         gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTc                                                                            >  1:3837257/1‑100 (MQ=255)
                          agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                           <  1:3247630/100‑1 (MQ=255)
                          agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                           <  1:3630836/100‑1 (MQ=255)
                           ggAACATGTATGGGTGGTATCAGTATTTGTCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAATTGCTCGGCTCCATCGGTTCCg                                                                          <  1:1682320/100‑1 (MQ=255)
                           ggAACATGTATGGGTACTATCAGTATTTGGCAGCTATTAATTATTGCCGTCATCGTTGTACTGCTTTTTGTCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                           <  1:2309517/99‑1 (MQ=255)
                            gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCGAAAAAGca                                                                                           <  1:3949035/82‑2 (MQ=255)
                            gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGa                                                                         >  2:3844233/1‑100 (MQ=255)
                             aaCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGAt                                                                        >  1:4205349/1‑100 (MQ=255)
                              aCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                       <  2:3818627/100‑1 (MQ=255)
                               cATGTATGGGTGGTATCAGTCTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                      >  1:2696332/1‑100 (MQ=255)
                                 tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                     >  1:568080/1‑99 (MQ=255)
                                 tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                    <  2:119250/100‑1 (MQ=255)
                                  gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                   >  2:3052089/1‑100 (MQ=255)
                                    aTGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTg                                                                 <  2:506274/100‑1 (MQ=255)
                                      gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCg                                                               >  2:1470900/1‑100 (MQ=255)
                                        gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTc                                                             <  2:2710622/100‑1 (MQ=255)
                                          ggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                           >  2:345268/1‑100 (MQ=255)
                                               cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                      >  2:2227680/1‑100 (MQ=255)
                                                aGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAg                                                     <  2:1108085/100‑1 (MQ=255)
                                                  tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                   <  2:2043019/100‑1 (MQ=255)
                                                  tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGGGCGTCGATCAAAGGc                                                   <  1:736841/100‑1 (MQ=255)
                                                    tttGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCtt                                                 <  2:3762758/100‑1 (MQ=255)
                                                      tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTa                                               >  1:78310/1‑100 (MQ=255)
                                                      tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTa                                               >  2:4172245/1‑100 (MQ=255)
                                                        gCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                             <  2:2371845/100‑1 (MQ=255)
                                                        gCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                             <  2:4227381/100‑1 (MQ=255)
                                                         cAGTTATTGATTGTTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaa                                            >  2:2892462/1‑100 (MQ=255)
                                                         cAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTTTTTTTTTTAAGCTCGGCTCCATCGGTTCCGATCTTGGTGTGTCGATCAAATGCatt                                                >  1:3914795/1‑96 (MQ=255)
                                                          aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaa                                           >  1:2907272/1‑100 (MQ=255)
                                                           gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                          >  1:4192890/1‑100 (MQ=255)
                                                           gTTATTGATTATTGCCGTCATCGTTGTACTGCTTGTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGTTCTTTGTGCTTCGATCAAAGGCTTTaaaaaa                                          <  1:2192678/100‑1 (MQ=255)
                                                            ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                         <  2:3480169/100‑1 (MQ=255)
                                                             tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACTAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGc                                        <  1:172860/100‑1 (MQ=255)
                                                             tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAATTCTCGGCTCTATCGGTTCAGGTCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                         <  1:1097665/99‑1 (MQ=255)
                                                             tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGc                                        <  1:2240050/100‑1 (MQ=255)
                                                              aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCATAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                       >  1:2610725/1‑100 (MQ=255)
                                                              aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCGGATCTTGTTGCGTCGATCAAAGGCTCTCAAAAAGc                                        >  1:2162802/1‑99 (MQ=255)
                                                              aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                       >  1:1160180/1‑100 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGGACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    <  1:1319771/100‑1 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                     >  1:3074762/1‑99 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    >  2:4292645/1‑100 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    >  2:3214690/1‑100 (MQ=255)
                                                                 gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                    <  2:556822/100‑1 (MQ=255)
                                                                   ttattGCCGGCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCTGCTCCATCGTATCCGATCTTGGTGCGTCAATCAAAGGCTTTAAAAAAGCAATg                                    <  1:3282208/98‑1 (MQ=255)
                                                                    tattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGc                                 >  1:3773330/1‑100 (MQ=255)
                                                                     attGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCg                                >  2:2424071/1‑100 (MQ=255)
                                                                      ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                               >  1:4118769/1‑100 (MQ=255)
                                                                       tGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAGAGCTCGGCACCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                              <  1:2821606/100‑1 (MQ=255)
                                                                         ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                            >  2:283526/1‑100 (MQ=255)
                                                                          cGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                            <  1:4292614/99‑1 (MQ=255)
                                                                          cGTCAGCGTTGTACTGCTGTTTGGCATCAAAATGTTCGGCTCCATCTGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatgat                           <  1:2657392/100‑1 (MQ=255)
                                                                            tCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTTGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGTTTTAAAAAAGCAATGAGCGAtgatga                         <  1:930684/100‑1 (MQ=255)
                                                                               ttttggTACTGCTTTTTGGCACCAAAAAGATCGGCTCCATCGGTTCCGATATTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAcc                      <  1:1552181/95‑1 (MQ=255)
                                                                               tCGTTGTACTGCTTTTTTGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAAcc                      >  1:1937604/1‑100 (MQ=255)
                                                                                  ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCaaa                   >  2:1593607/1‑100 (MQ=255)
                                                                                    gccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                 >  2:3611876/6‑100 (MQ=255)
                                                                                    gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                 >  2:2363968/1‑100 (MQ=255)
                                                                                     tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGTACCAAAGCa                >  1:3043925/1‑100 (MQ=255)
                                                                                     ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa                >  1:1138692/5‑100 (MQ=255)
                                                                                      cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg               >  2:1778891/4‑100 (MQ=255)
                                                                                       gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa                >  1:768158/3‑98 (MQ=255)
                                                                                       gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  1:3054631/3‑100 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  1:3143430/1‑100 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  1:1830329/1‑100 (MQ=255)
                                                                                       cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg              >  2:2509307/1‑100 (MQ=255)
                                                                                        tGCTTTTAGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGACATAAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGa             <  1:1538739/100‑1 (MQ=255)
                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt            >  1:3325599/1‑100 (MQ=255)
                                                                                         gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCCATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt            >  2:1712815/1‑100 (MQ=255)
                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc  >  1:3207018/1‑100 (MQ=255)
                                                                                                   aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc  <  2:923258/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: