Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,539,427 T→G S681S (TCT→TCG yhgF → putative transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,4270TG93.1% 481.0 / 15.7 160S681S (TCT→TCGyhgFputative transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base T (7/4);  new base G (61/87);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 161 total aligned reads displayed.

AACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAA  >  NC_000913/3539330‑3539523
                                                                                                 |                                                                                                
aaCCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTc                                                                                                <  1:907880/100‑1 (MQ=255)
 aCCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCa                                                                                               <  1:1887708/100‑1 (MQ=255)
 aCCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCa                                                                                               <  1:1886193/100‑1 (MQ=255)
  ccGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCAGCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGATGGCCTGGTTCACATCTCGTCAt                                                                                              <  2:2552833/100‑1 (MQ=255)
  ccGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCAt                                                                                              <  1:1468585/100‑1 (MQ=255)
   cGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCAtt                                                                                             >  2:3339342/1‑100 (MQ=255)
    gggTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTg                                                                                            <  1:1578673/100‑1 (MQ=255)
     ggTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTgn                                                                                           <  1:3249999/100‑2 (MQ=255)
       tATGATCCTCGAAGGCGCAGTGACCGTCGTCACCACGGTTGGCGCGGTGGTCGATATGGGCGTGCTTCTGGAGGCCCTGGTTCACGTCGGGTCATTGTCg                                                                                         <  2:1265522/100‑1 (MQ=255)
       tATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCg                                                                                         <  1:3706816/100‑1 (MQ=255)
       tATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCg                                                                                         <  2:3663758/100‑1 (MQ=255)
       tATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCg                                                                                         <  2:1251658/100‑1 (MQ=255)
        aTGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGa                                                                                        >  1:3333277/1‑100 (MQ=255)
         tGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGaa                                                                                       >  1:2696240/1‑100 (MQ=255)
          gATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAAc                                                                                      <  2:3777379/100‑1 (MQ=255)
          gATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAAc                                                                                      <  2:2917099/100‑1 (MQ=255)
            tCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACaa                                                                                    >  1:1343021/1‑100 (MQ=255)
            tCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACaa                                                                                    <  1:3181352/100‑1 (MQ=255)
             ccTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAg                                                                                   >  2:1471770/1‑100 (MQ=255)
              cTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGt                                                                                  <  2:3915232/100‑1 (MQ=255)
              cTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGt                                                                                  >  1:3560731/1‑100 (MQ=255)
                 gAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTtg                                                                               <  2:172473/100‑1 (MQ=255)
                 gAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTtg                                                                               >  1:782643/1‑100 (MQ=255)
                   aGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTtgtg                                                                             <  1:2094589/100‑1 (MQ=255)
                    ggCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTgg                                                                            <  1:195430/100‑1 (MQ=255)
                    ggCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTgg                                                                            <  1:544384/100‑1 (MQ=255)
                    ggCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTgg                                                                            <  2:3116781/100‑1 (MQ=255)
                      cgcAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGaa                                                                          <  1:1948829/100‑1 (MQ=255)
                          gTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATc                                                                      <  2:3857806/100‑1 (MQ=255)
                           tGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATcc                                                                     <  1:130407/100‑1 (MQ=255)
                           tGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATcc                                                                     <  2:3448110/100‑1 (MQ=255)
                            gACCACCGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCg                                                                    <  1:415801/100‑1 (MQ=255)
                            gACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCg                                                                    <  2:525816/100‑1 (MQ=255)
                            gACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCg                                                                    <  1:1421541/100‑1 (MQ=255)
                            gACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCg                                                                    >  1:3575105/1‑100 (MQ=255)
                              ccAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCa                                                                  <  1:561461/100‑1 (MQ=255)
                              ccAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCa                                                                  >  1:1518483/1‑100 (MQ=255)
                              ccAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCa                                                                  <  2:2484698/100‑1 (MQ=255)
                                 aCGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATAc                                                               >  2:1530926/1‑100 (MQ=255)
                                  cGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATAcc                                                              <  1:1971504/100‑1 (MQ=255)
                                    tCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgt                                                            >  1:264807/1‑100 (MQ=255)
                                     cACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtg                                                           >  2:2665027/1‑100 (MQ=255)
                                     cACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtg                                                           >  2:2063266/1‑100 (MQ=255)
                                       ccAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtggt                                                         >  2:3013234/1‑100 (MQ=255)
                                        cAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtggtg                                                        >  1:2413133/1‑100 (MQ=255)
                                        cAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtggtg                                                        <  1:2926826/100‑1 (MQ=255)
                                            tttGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAg                                                    <  2:2275478/100‑1 (MQ=255)
                                            tttGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAg                                                    >  2:2635226/1‑100 (MQ=255)
                                            tttGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAg                                                    <  2:3272854/100‑1 (MQ=255)
                                             ttGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGc                                                   <  1:2973940/100‑1 (MQ=255)
                                               ggCGCGTTTGTCGATATTGGCTTGCCTCCGGACGGCCGGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAATATCCGCATACCGTGGTGAAAGCgg                                                 >  2:3974733/1‑100 (MQ=255)
                                                 cgcgTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGc                                               <  1:2545715/100‑1 (MQ=255)
                                                  gcgTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCg                                              <  1:3928585/100‑1 (MQ=255)
                                                  gcgTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCg                                              <  1:720104/100‑1 (MQ=255)
                                                   cgTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGa                                             >  1:2096356/1‑100 (MQ=255)
                                                    gTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGAc                                            <  2:1557596/100‑1 (MQ=255)
                                                       tGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACAtt                                         >  2:3221621/1‑100 (MQ=255)
                                                         tCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATtgt                                       <  1:830745/100‑1 (MQ=255)
                                                          cGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATtgtg                                      >  1:1418820/1‑100 (MQ=255)
                                                            atatTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGaa                                    <  1:1014237/100‑1 (MQ=255)
                                                             tatTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAg                                   <  2:3393178/100‑1 (MQ=255)
                                                              cttGGCGTGCACCAGTACTGCCTGGTTCACATCTCGTCATTGTCCAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAgg                                  <  2:2520223/99‑1 (MQ=255)
                                                              atTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAgg                                  >  2:3142053/1‑100 (MQ=255)
                                                              atTGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAgg                                  <  1:1533172/100‑1 (MQ=255)
                                                               ttGGCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGt                                 >  1:1680874/1‑100 (MQ=255)
                                                                 ggCGTGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGa                               <  1:1631989/100‑1 (MQ=255)
                                                                    gTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAg                            >  2:2069958/1‑100 (MQ=255)
                                                                     tGCATGAGGACGGTCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGt                           <  2:3724073/100‑1 (MQ=255)
                                                                     tGCATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGt                           >  1:535028/1‑100 (MQ=255)
                                                                       cATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGc                         >  1:2276489/1‑100 (MQ=255)
                                                                       cATCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGc                         <  2:1093704/100‑1 (MQ=255)
                                                                       cATCAGGACGGCCTGCTTCACATCTCATCATTGTCGAGCAAGTTAGTGGAAGATCCGCATACCGTCTTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGc                         <  1:2321697/100‑1 (MQ=255)
                                                                         tCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTg                       >  2:2399132/1‑100 (MQ=255)
                                                                         tCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTg                       <  1:2455768/100‑1 (MQ=255)
                                                                         tCAGGACGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTg                       <  1:2739308/100‑1 (MQ=255)
                                                                         tCAGGACGGCCTGGTTCACATCTAGTGATTGTCGAAACCGTCTGTGGCGGATCCGCATACAGTAGTGAAAGCGGGCGCCATTGTGAAGGTGATAGTGCTg                       <  2:3995706/100‑1 (MQ=255)
                                                                               cGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTTAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTg                 <  2:3076419/100‑1 (MQ=255)
                                                                               cGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTg                 <  1:1248740/100‑1 (MQ=255)
                                                                               cGGCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTg                 >  1:2302228/1‑100 (MQ=255)
                                                                                ggCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTgg                <  1:2845627/100‑1 (MQ=255)
                                                                                ggCCTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTgg                <  1:2421218/100‑1 (MQ=255)
                                                                                  ccTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGAt              >  2:3250328/1‑100 (MQ=255)
                                                                                  ccTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGAt              >  2:3251389/1‑100 (MQ=255)
                                                                                   cTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATc             >  1:100966/1‑100 (MQ=255)
                                                                                   cTGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATc             <  2:535629/100‑1 (MQ=255)
                                                                                    tGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCt            >  2:1415261/1‑100 (MQ=255)
                                                                                    tGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCt            <  1:1195190/100‑1 (MQ=255)
                                                                                    tGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCt            <  2:41175/100‑1 (MQ=255)
                                                                                    tGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCt            <  1:3343883/100‑1 (MQ=255)
                                                                                    tGGTTCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCt            <  2:247578/100‑1 (MQ=255)
                                                                                        tCACATCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAg        <  1:3703158/100‑1 (MQ=255)
                                                                                        tCACAGCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAg        <  2:1274684/100‑1 (MQ=255)
                                                                                         cacaTCTCGTCATTGTCGCACAAGGTTGTCGAAGATCCGCATACCGTGGTGAACGCGGGCGACATTGTGAAGGTGAAAGTGCTGGTGGTGGATCTTCAGc       >  2:2237368/1‑100 (MQ=255)
                                                                                           caTCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGt     >  1:3133867/1‑100 (MQ=255)
                                                                                           caTCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGt     >  1:3395215/1‑100 (MQ=255)
                                                                                           caTCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGt     >  1:1087961/1‑100 (MQ=255)
                                                                                           caTCTCGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGt     >  2:1433835/1‑100 (MQ=255)
                                                                                             tctcTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTaa   >  2:2525941/1‑100 (MQ=255)
                                                                                             tctcGTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTaa   >  2:40490/1‑100 (MQ=255)
                                                                                              ctcgTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTaaa  >  2:798189/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAA  >  NC_000913/3539330‑3539523

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