Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,007 C→T 18.2% 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,0070CT18.2% 442.1 / 75.7 202G21G (GGC→GGTtatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base C (92/71);  new base T (15/22);  total (107/95)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

Only 100 of 207 total aligned reads displayed.

GCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGT  >  NC_000913/4021908‑4022106
                                                                                                   |                                                                                                   
gCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGc                                                                                                     <  2:1260884/100‑1 (MQ=255)
 cTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                    >  2:4035428/1‑100 (MQ=255)
 cTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCa                                                                                                    <  1:2422235/100‑1 (MQ=255)
   ttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCAcc                                                                                                  <  1:2021303/100‑1 (MQ=255)
    tGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCa                                                                                                 >  2:1407788/1‑100 (MQ=255)
     gTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaa                                                                                                >  2:3861590/1‑100 (MQ=255)
      tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaa                                                                                               >  2:1010332/1‑100 (MQ=255)
       ttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaa                                                                                              <  2:3430819/100‑1 (MQ=255)
        tAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaaa                                                                                             >  2:3672147/1‑100 (MQ=255)
         aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAg                                                                                            >  1:88910/1‑100 (MQ=255)
          atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCGGTCATCGTGGTACTGCTTTTTGGCACCAAAAAGc                                                                                           >  2:732085/1‑100 (MQ=255)
          atcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                           <  2:3469643/100‑1 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:3598127/100‑2 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:2387137/100‑2 (MQ=255)
           tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                          <  2:1090722/100‑2 (MQ=255)
           tcatcatcTACAACAAAGGAACATGTATGGGTGATATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCt                                                                                          >  1:2108280/1‑100 (MQ=255)
            catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGTCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcag                                                                                         <  1:2482065/100‑3 (MQ=255)
              tcatcTTCCGCAGTGGAACATGTATGGTTGGGATCAGGATTTGTCAGTTATTGATTATTGCCGTTATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagta                                                                                       <  2:2562794/100‑5 (MQ=255)
               catcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGc                                                                                      <  1:1010328/100‑1 (MQ=255)
                atcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGGTGTACTGCTTTTTGGCACCAAAAAGCTCGGCt                                                                                     <  2:2617599/100‑1 (MQ=255)
                  cTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGcagtacaac                                                                                   <  2:1848182/100‑9 (MQ=38)
                     ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  2:4200774/1‑100 (MQ=255)
                     ccACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATc                                                                                >  1:802100/1‑100 (MQ=255)
                      cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCg                                                                               >  2:3816681/1‑100 (MQ=255)
                      cacaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCg                                                                               >  2:951408/1‑100 (MQ=255)
                       acaGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCgg                                                                              >  1:1362355/1‑100 (MQ=255)
                          gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGtt                                                                            <  1:1340415/99‑1 (MQ=255)
                          gagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTc                                                                           >  1:3837257/1‑100 (MQ=255)
                           agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:3630836/100‑1 (MQ=255)
                           agGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:3247630/100‑1 (MQ=255)
                            ggAACATGTATGGGTGGTATCAGTATTTGTCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAATTGCTCGGCTCCATCGGTTCCg                                                                         <  1:1682320/100‑1 (MQ=255)
                            ggAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCg                                                                         >  2:3239037/1‑100 (MQ=255)
                            ggAACATGTATGGGTACTATCAGTATTTGGCAGCTATTAATTATTGCCGTCATCGTTGTACTGCTTTTTGTCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                          <  1:2309517/99‑1 (MQ=255)
                             gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCGAAAAAGca                                                                                          <  1:3949035/82‑2 (MQ=255)
                              aaCATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGAt                                                                       >  1:4205349/1‑100 (MQ=255)
                                cATGTATGGGTGGTATCAGTCTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCt                                                                     >  1:2696332/1‑100 (MQ=255)
                                 aTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                    <  2:553514/100‑1 (MQ=255)
                                  tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                   <  1:2190830/100‑1 (MQ=255)
                                  tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                   <  2:119250/100‑1 (MQ=255)
                                  tGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                   <  2:2025765/100‑1 (MQ=255)
                                   gTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                  <  2:3058117/100‑1 (MQ=255)
                                    tATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                 >  2:2237993/1‑100 (MQ=255)
                                       gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCg                                                              >  2:1470900/1‑100 (MQ=255)
                                        ggtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGt                                                             >  2:2367763/1‑100 (MQ=255)
                                             tATCAGTATTTGGCAGTTATTGATGATTGCCGTCAGCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                        <  1:3181072/100‑1 (MQ=255)
                                                cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                     >  2:2133444/1‑100 (MQ=255)
                                                cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                     >  2:4341398/1‑100 (MQ=255)
                                                cAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                     >  2:441952/1‑100 (MQ=255)
                                                 aGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAg                                                    <  2:1108085/100‑1 (MQ=255)
                                                     tttGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCtt                                                <  2:4251900/100‑1 (MQ=255)
                                                       tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTa                                              >  1:78310/1‑100 (MQ=255)
                                                       tGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTa                                              <  2:3042393/100‑1 (MQ=255)
                                                         gCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                            <  2:638367/100‑1 (MQ=255)
                                                          cAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTTTTTTTTTTAAGCTCGGCTCCATCGGTTCCGATCTTGGTGTGTCGATCAAATGCatt                                               >  1:3914795/1‑96 (MQ=255)
                                                           aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaa                                          <  2:2204512/100‑1 (MQ=255)
                                                           aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaa                                          >  1:2907272/1‑100 (MQ=255)
                                                            gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                         >  2:201894/1‑100 (MQ=255)
                                                            gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                         >  1:4192890/1‑100 (MQ=255)
                                                             ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                        <  2:2721141/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)
                                                               aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTACTAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                      >  1:742641/1‑100 (MQ=255)
                                                               aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCATAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                      >  1:2610725/1‑100 (MQ=255)
                                                               aTTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                      >  1:1160180/1‑100 (MQ=255)
                                                               aTTGATTATTACCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                      >  1:610144/1‑100 (MQ=255)
                                                                  gATTATTGCCGTCATCGTTGTACTGCTTTTTGGGACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATg                                   <  1:1319771/100‑1 (MQ=255)
                                                                  gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAAt                                    >  1:3074762/1‑99 (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)
                                                                        tGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                             >  2:341048/1‑100 (MQ=255)
                                                                          ccGTCATCGTTGTTCTGCTTTTTGGCATCAAAAAGCTCGGCTTCATCGGTTCCGATCTTGGTGCGTCTATCAAAGGCTTTAAAAAAGCAATGAGCg                               <  1:768745/96‑1 (MQ=255)
                                                                          ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgatga                           >  2:1985935/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)
                                                                                      tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGTACCAAAGCa               >  1:3043925/1‑100 (MQ=255)
                                                                                      tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               <  1:3472959/100‑1 (MQ=255)
                                                                                      ccgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCa               >  1:1138692/5‑100 (MQ=255)
                                                                                       cgagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg              >  2:3832518/4‑100 (MQ=255)
                                                                                       aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAg              <  2:2835595/100‑1 (MQ=255)
                                                                                        gagCTTTTTGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:3054631/3‑100 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             <  2:3629703/100‑1 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  2:2509307/1‑100 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             >  1:1830329/1‑100 (MQ=255)
                                                                                        cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTACCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             <  2:3096287/100‑1 (MQ=255)
                                                                                         tGCTTTTAGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGACATAAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGa            <  1:1538739/100‑1 (MQ=255)
                                                                                          gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAgg             <  2:3542554/98‑1 (MQ=255)
                                                                                          gCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt           >  1:3325599/1‑100 (MQ=255)
                                                                                            tttttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaa         >  2:558050/1‑100 (MQ=255)
                                                                                               ttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAc      <  2:3358021/100‑1 (MQ=255)
                                                                                               ttGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAAc      <  1:2227671/100‑1 (MQ=255)
                                                                                                 ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCa    >  2:1713065/1‑100 (MQ=255)
                                                                                                   cACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGt  >  2:4134852/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: