Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,838,723 T→C 16.7% L159S (TTG→TCG)  yicL → EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,7230TC16.7% 185.8 / 29.4 83L159S (TTG→TCG) yicLEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (6/8);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

CATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAG  >  NC_000913/3838625‑3838822
                                                                                                  |                                                                                                   
cATTCTGGTTTTCTTCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCgttg                                                                                                    <  2:443143/100‑1 (MQ=255)
cATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCgttg                                                                                                    >  2:2123339/1‑100 (MQ=255)
  ttCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCgttgtt                                                                                                  <  1:152780/100‑1 (MQ=255)
   tCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTc                                                                                                 <  2:628108/100‑1 (MQ=255)
    cTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCt                                                                                                >  1:414046/1‑100 (MQ=255)
    cTGGTTTTCTGCGCGATTTTGACATCGCTGGTCGGGACTTTTTCATTGGTGACACACGGTAATCTGACGTCATTATCGATCTCTCCTGCCGCGTTGATCt                                                                                                >  2:58583/1‑100 (MQ=255)
    atGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCt                                                                                                >  2:1132744/2‑100 (MQ=255)
     tGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTg                                                                                               >  1:2131634/1‑100 (MQ=255)
       gTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTggg                                                                                             <  2:1627008/100‑1 (MQ=255)
        ttttCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGGCATTATCGATCTCTCCTGCCGCGTTGTTCTgggg                                                                                            >  1:575450/1‑100 (MQ=255)
         tttCTGCGCTATTTTGACATCGCTGGGCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTgggg                                                                                            <  1:1137311/99‑1 (MQ=255)
            cTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCAtt                                                                                        >  1:248642/1‑100 (MQ=255)
             tGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTg                                                                                       >  1:981293/1‑100 (MQ=255)
              gcgcTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGc                                                                                      >  2:767301/1‑100 (MQ=255)
              gcgcTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGCGGCAATg                                                                                       >  2:2115553/1‑99 (MQ=255)
               cgcTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGcc                                                                                     <  1:1518764/100‑1 (MQ=255)
                gcTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCt                                                                                    <  2:957993/100‑1 (MQ=255)
                gcTATTTTGACATCGCTGGTCGGGACTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTc                                                                                   <  1:2322297/100‑1 (MQ=255)
                 cTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTc                                                                                   >  2:2284750/1‑100 (MQ=255)
                  tATTTTGACATCGCTTGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCg                                                                                  <  2:1087437/100‑1 (MQ=255)
                  tATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCg                                                                                  <  1:2367273/100‑1 (MQ=255)
                  tATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCg                                                                                  >  2:1367205/1‑100 (MQ=255)
                   aTTTTGACATCGCTGTTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCgg                                                                                 <  1:2347044/100‑1 (MQ=255)
                   aTTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCgg                                                                                 >  1:532736/1‑100 (MQ=255)
                     tttGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCa                                                                               <  1:1882442/100‑1 (MQ=255)
                        gACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCAttt                                                                            <  1:786583/100‑1 (MQ=255)
                         aCATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtg                                                                           <  1:1696927/100‑1 (MQ=255)
                          cATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgc                                                                          <  1:2074956/100‑1 (MQ=255)
                          cATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgc                                                                          >  2:268215/1‑100 (MQ=255)
                            tCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgct                                                                         <  1:1804099/99‑1 (MQ=255)
                             cGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgctgc                                                                       <  2:2001910/100‑1 (MQ=255)
                             cGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgctgc                                                                       <  2:2005549/100‑1 (MQ=255)
                             cGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTtgctgc                                                                       <  2:974079/100‑1 (MQ=255)
                                tGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTtgctgctgc                                                                    >  1:2305074/1‑100 (MQ=255)
                                      ggACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCGCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCta                                                              <  2:2235782/100‑1 (MQ=255)
                                       gACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTCCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCtat                                                             <  2:1275367/100‑1 (MQ=255)
                                       gACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGGTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCtat                                                             <  2:1470022/100‑1 (MQ=255)
                                         cTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATac                                                           >  2:1047884/1‑100 (MQ=255)
                                         cTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATac                                                           >  1:1854418/1‑100 (MQ=255)
                                         cTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATac                                                           >  1:2181155/1‑100 (MQ=255)
                                          ttttttATTGTTGACACACGGTAATCCGACGTCATGATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATacc                                                          <  1:2130603/100‑1 (MQ=255)
                                          ttttttATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATacc                                                          <  1:1382967/100‑1 (MQ=255)
                                           tttttATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATacca                                                         >  2:616619/1‑100 (MQ=255)
                                             tttATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATaccacc                                                       <  2:1289628/100‑1 (MQ=255)
                                              ttattaGTCACACACGGTAATCTGCCGTACTTATCGATCTCTCCTGCCTCGTCGTTCTGGGGCATTGCCTCGGCTTTTGCTGCTGCATACTGTactccc                                                       <  2:730118/99‑5 (MQ=255)
                                              ttATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCt                                                      <  2:95577/100‑1 (MQ=255)
                                                aTTGGTGACACACAGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTAt                                                    >  2:1472251/1‑100 (MQ=255)
                                                 ttGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATc                                                   <  1:268854/100‑1 (MQ=255)
                                                  tGGTGACACACGGTAATCCGCCGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATcc                                                  >  1:326562/1‑100 (MQ=255)
                                                   ggTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATccc                                                 <  1:942985/100‑1 (MQ=255)
                                                   ggTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATccc                                                 <  1:2284771/100‑1 (MQ=255)
                                                     tGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTc                                               <  1:393109/100‑1 (MQ=255)
                                                     tGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTc                                               >  2:1180474/1‑100 (MQ=255)
                                                     tGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTc                                               <  1:563838/100‑1 (MQ=255)
                                                      gccactcGGTAATCCTACGTCATTATCGTTCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCa                                              <  2:632962/98‑1 (MQ=255)
                                                      gACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCa                                              >  2:188922/1‑100 (MQ=255)
                                                        cacacGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATCCCACCTATCCCTCAAc                                            >  1:2233002/1‑100 (MQ=255)
                                                        cacacGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAAc                                            >  2:347205/1‑100 (MQ=255)
                                                        cacacGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAAc                                            <  1:1976395/100‑1 (MQ=255)
                                                        cacacGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAAc                                            >  2:1613030/1‑100 (MQ=255)
                                                        cacacGGTAATCCGACGTCATTATCGATCTATCATGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCAATCCCTCAAc                                            >  2:1572085/1‑100 (MQ=255)
                                                                aaTCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTg                                    >  2:1556392/1‑100 (MQ=255)
                                                                aaTCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTg                                    >  1:962527/1‑100 (MQ=255)
                                                                 aTCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGc                                   <  1:2367201/100‑1 (MQ=255)
                                                                  tCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGcc                                  >  2:2138246/1‑100 (MQ=255)
                                                                   ccGACGTCATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCCTTTGCTGCTGCATTCTATACCACCTATCCCTCAACACTAATTGccc                                 >  2:687099/1‑100 (MQ=255)
                                                                     gACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGAGCATTGCCTCGGCATTTGCTGATGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGc                               >  1:1388554/1‑100 (MQ=255)
                                                                       cGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTa                             <  1:2137478/100‑1 (MQ=255)
                                                                         tCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATg                           <  1:782539/100‑1 (MQ=255)
                                                                          cATTATCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATgg                          >  1:2318428/1‑100 (MQ=255)
                                                                           aTTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGc                         <  2:2049114/100‑1 (MQ=255)
                                                                            ttATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCa                        <  1:2008845/100‑1 (MQ=255)
                                                                             tATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCAc                       <  1:80687/100‑1 (MQ=255)
                                                                              aTCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACg                      >  1:993932/1‑100 (MQ=255)
                                                                              aTCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACg                      <  2:1551282/100‑1 (MQ=255)
                                                                              aTCGATCTCTCCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACg                      <  2:1786677/100‑1 (MQ=255)
                                                                               tCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGt                     >  1:1229908/1‑100 (MQ=255)
                                                                                cGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGtt                    >  2:92946/1‑100 (MQ=255)
                                                                                 gATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTa                   >  1:2279355/1‑100 (MQ=255)
                                                                                 gATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTa                   >  1:1526918/1‑100 (MQ=255)
                                                                                   tctctcCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTAcc                 <  1:891202/100‑1 (MQ=255)
                                                                                     tctcCTGCCGCGTCGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAg               >  1:1069786/1‑100 (MQ=255)
                                                                                      ctcCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAgt              >  1:409788/1‑100 (MQ=255)
                                                                                      ctcCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAg               <  2:701394/99‑1 (MQ=255)
                                                                                              gcgTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCt      >  2:303336/1‑100 (MQ=255)
                                                                                               cgTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTg     <  1:2138263/100‑1 (MQ=255)
                                                                                                gttgttCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTgg    >  1:1368950/1‑100 (MQ=255)
                                                                                                  tgttCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAg  <  1:767303/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CATTCTGGTTTTCTGCGCTATTTTGACATCGCTGGTCGGGACTTTTTTATTGGTGACACACGGTAATCCGACGTCATTATCGATCTCTCCTGCCGCGTTGTTCTGGGGCATTGCCTCGGCATTTGCTGCTGCATTCTATACCACCTATCCCTCAACGCTAATTGCCCGCTATGGCACGTTACCAGTCGTCGGCTGGAG  >  NC_000913/3838625‑3838822

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: