Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,338,928 C→A 100% L107F (TTG→TTT mlaE ← ABC transporter maintaining OM lipid asymmetry, inner membrane permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,338,9280CA98.7% 521.8 / ‑1.5 153L107F (TTG→TTTmlaEABC transporter maintaining OM lipid asymmetry, inner membrane permease protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (70/79);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 159 total aligned reads displayed.

GCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAG  >  NC_000913/3338831‑3339023
                                                                                                 |                                                                                               
gCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCaaaaaa                                                                                               >  1:1459531/1‑100 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              >  2:210001/1‑100 (MQ=255)
  ggATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaaca                                                                                             >  2:1123766/1‑100 (MQ=255)
  ggATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaaca                                                                                             <  1:2742062/100‑1 (MQ=255)
  ggATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGAGCCGAACCAGCACGCCCGGCAaaaaaca                                                                                             >  2:847971/1‑100 (MQ=255)
    aTCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaacaac                                                                                           <  2:3381274/100‑1 (MQ=255)
      ncACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAAcgc                                                                                         >  1:3774896/2‑100 (MQ=255)
       cACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAAcgcg                                                                                        >  1:3802821/1‑100 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       <  2:3155318/100‑1 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       >  2:432504/1‑100 (MQ=255)
          cGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCa                                                                                     <  1:3594994/100‑1 (MQ=255)
             catcatCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAacc                                                                                  >  2:1785509/1‑100 (MQ=255)
                  tctcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCgg                                                                             >  2:2348222/1‑100 (MQ=255)
                   ctcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGc                                                                            >  2:1509748/1‑100 (MQ=255)
                   ctcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGc                                                                            >  2:3039552/1‑100 (MQ=255)
                    tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCATCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcc                                                                           <  2:2733084/100‑1 (MQ=255)
                    tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcc                                                                           >  2:3654378/1‑100 (MQ=255)
                    tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcc                                                                           >  2:2383993/1‑100 (MQ=255)
                    tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcc                                                                           >  1:1837933/1‑100 (MQ=255)
                     ccATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGccc                                                                          >  1:1935424/1‑100 (MQ=255)
                      cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAaccacc                                                                               <  1:1842146/94‑1 (MQ=255)
                      cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcccc                                                                         <  1:2221619/100‑1 (MQ=255)
                      cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcccc                                                                         <  1:3878809/100‑1 (MQ=255)
                        tACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAg                                                                       <  1:2291826/100‑1 (MQ=255)
                         aCTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAg                                                                       <  2:3446872/99‑1 (MQ=255)
                         aCTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGt                                                                      <  2:2092855/100‑1 (MQ=255)
                          cTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGt                                                                      <  2:3282918/99‑1 (MQ=255)
                            ggAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCa                                                                   <  2:625908/100‑1 (MQ=255)
                            ggAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCa                                                                   >  2:1195940/1‑100 (MQ=255)
                            ggAGAGTTGCTCTGTACCGCGCATCAGGACGATTTCTGCGGGTAGCGCCGAACCAGAAAAAAAAAAAAAAAACAACGCGGCAACCACCGGcc                                                                           <  1:1046318/92‑1 (MQ=255)
                              agagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACg                                                                 <  1:1695973/100‑1 (MQ=255)
                               gagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCAAGAAAAGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGc                                                                <  1:1401271/100‑1 (MQ=255)
                                agTTGCTCTGTAGCGCGCATCAGGCGGATTTCTGAGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGc                                                                <  1:1361813/99‑1 (MQ=255)
                                  ttGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGt                                                             >  2:3099055/1‑100 (MQ=255)
                                     ctctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCCAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGc                                                          <  1:777629/100‑1 (MQ=255)
                                       ctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTa                                                            <  1:2759046/96‑1 (MQ=255)
                                       ctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGa                                                        <  2:354269/100‑1 (MQ=255)
                                         gTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATa                                                      <  1:3227584/100‑1 (MQ=255)
                                            gcgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACg                                                   <  2:570921/100‑1 (MQ=255)
                                             cgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGc                                                  >  1:1647267/1‑100 (MQ=255)
                                             cgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGc                                                  >  1:1625978/1‑100 (MQ=255)
                                             cgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGc                                                  >  2:112878/1‑100 (MQ=255)
                                                gcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccac                                               <  2:3774927/100‑1 (MQ=255)
                                                gcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccac                                               <  1:112876/100‑1 (MQ=255)
                                                 cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                              <  1:1009857/100‑1 (MQ=255)
                                                 cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                              >  2:3584174/1‑100 (MQ=255)
                                                 cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                              <  2:1053733/100‑1 (MQ=255)
                                                   tCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAg                                            <  2:837276/100‑1 (MQ=255)
                                                     aGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACTGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCa                                          <  1:3356581/100‑1 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCAt                                         >  2:1597426/1‑100 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCAt                                         >  2:609916/1‑100 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACACCACCAACa                                          >  2:1349547/1‑99 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        <  2:1647276/100‑1 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        >  1:3549028/1‑100 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCAGACCCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCa                                          <  1:968555/98‑1 (MQ=255)
                                                        ccGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAc                                       >  1:3362496/1‑100 (MQ=255)
                                                        ccGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAc                                       <  1:1073815/100‑1 (MQ=255)
                                                        ccGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAc                                       <  1:105720/100‑1 (MQ=255)
                                                          gATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAccc                                     >  2:2945052/1‑100 (MQ=255)
                                                          gATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAccc                                     >  2:358437/1‑100 (MQ=255)
                                                           aTTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCa                                    <  2:2570609/100‑1 (MQ=255)
                                                            tttCTGCGGTTAGCGCCGAACCAGCACGCCCGTCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAg                                   >  2:3676543/1‑100 (MQ=255)
                                                              tCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGan                                 <  1:3810606/100‑2 (MQ=255)
                                                              tCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGAc                                 >  1:1205498/1‑100 (MQ=255)
                                                                 gCGGTTAGCGCCGAACCAGGACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTgg                              <  1:2366800/100‑1 (MQ=255)
                                                                 gCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTgg                              <  2:3046293/100‑1 (MQ=255)
                                                                 gCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTgg                              <  2:2883425/100‑1 (MQ=255)
                                                                  cGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGt                             <  2:3552543/100‑1 (MQ=255)
                                                                   ggTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGtt                            >  1:2912926/1‑100 (MQ=255)
                                                                    gTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGttt                           <  1:1547287/100‑1 (MQ=255)
                                                                      tAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTcc                         <  1:3035329/100‑1 (MQ=255)
                                                                        gcgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGc                       <  2:1837942/100‑1 (MQ=255)
                                                                        gcgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGc                       >  1:1563360/1‑100 (MQ=255)
                                                                         cgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCa                      <  2:3686193/100‑1 (MQ=255)
                                                                         cgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCa                      >  1:3596673/1‑100 (MQ=255)
                                                                         cgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGAAGCGATAACGCCACCAGCATACCCAAACTGGTTTCCGCa                      >  1:1004831/1‑100 (MQ=255)
                                                                          gcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCAc                     <  1:988491/100‑1 (MQ=255)
                                                                           ccGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACt                    >  1:333152/1‑100 (MQ=255)
                                                                           ccGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACt                    <  2:1068165/100‑1 (MQ=255)
                                                                           ccGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACt                    >  1:2858934/1‑100 (MQ=255)
                                                                           ccGAACCAGCACGCCCGGCAAAAAACAACCCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAGCGCCACCAGCATACCCAGACTGGTTTCCGCACt                    <  1:2328679/100‑1 (MQ=255)
                                                                             gAACCATCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCACTAGCGCTAACGCCACCAGCATACCCAGACTGCCTTCCGc                       >  1:2143201/1‑95 (MQ=255)
                                                                               aCCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATaa                <  2:1016594/100‑1 (MQ=255)
                                                                                  aGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTg             >  2:1362724/1‑100 (MQ=255)
                                                                                   gCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTgg            >  1:1752134/1‑100 (MQ=255)
                                                                                   gCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTgg            <  2:1691826/100‑1 (MQ=255)
                                                                                    cACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCCGTTCCCGCAATAGCGATAACGACACCAGCATACCCAGACTGGTTTCCCAACTATAAGTGGt           >  1:1688662/1‑100 (MQ=255)
                                                                                      cGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCCGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCa         <  2:2725562/100‑1 (MQ=255)
                                                                                       gCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAg        <  2:2452010/100‑1 (MQ=255)
                                                                                        cccGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGa       >  1:708569/1‑100 (MQ=255)
                                                                                        cccGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGa       >  2:2248468/1‑100 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     <  2:1648701/100‑1 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     <  2:2504321/100‑1 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     <  1:1471779/100‑1 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     <  2:939521/100‑1 (MQ=255)
                                                                                           ggCAAAAAACAACGCGGCAACCACCGGCCTCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTCGTTTCCGCACTATAAGTGGTCAGAAc     <  1:2249408/99‑1 (MQ=255)
                                                                                           ggCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAcc    >  1:3747845/1‑100 (MQ=255)
                                                                                            gCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCa   <  2:2858957/100‑1 (MQ=255)
                                                                                            gCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCa   >  1:1884358/1‑100 (MQ=255)
                                                                                             cAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAg  <  1:1141567/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAG  >  NC_000913/3338831‑3339023

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: