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.3% 388.2 / ‑1.3 117L107F (TTG→TTTmlaEABC transporter maintaining OM lipid asymmetry, inner membrane permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (49/66);  total (49/68)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 126 total aligned reads displayed.

GCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAG  >  NC_000913/3338831‑3339023
                                                                                                 |                                                                                               
gCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCaaaaaa                                                                                               >  2:1859314/1‑100 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              >  1:2555721/1‑100 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              <  1:158061/100‑1 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              >  2:585908/1‑100 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              >  2:1564993/1‑100 (MQ=255)
 cGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCGGTATCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaac                                                                                              <  1:2013479/100‑1 (MQ=255)
    aTCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaacaac                                                                                           <  1:1364770/100‑1 (MQ=255)
    aTCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaaaacaac                                                                                           >  2:31755/1‑100 (MQ=255)
     tCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAACACAAcg                                                                                          >  2:1528801/1‑100 (MQ=255)
      ccACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGCACCAGCACGCCCGGCAAAAAACAAcgc                                                                                         >  2:2426583/1‑100 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       <  1:2733836/100‑1 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       <  2:459555/100‑1 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       <  1:78711/100‑1 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCgg                                                                                       <  1:2521485/100‑1 (MQ=255)
        aCCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCAAACACAACAAAAAACAACGCgg                                                                                       <  2:1648957/100‑1 (MQ=255)
         ccGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCTGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGc                                                                                      <  2:2134501/100‑1 (MQ=255)
         ccGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTTCGGTTAGCTCCGAACCAGCACGCCCGGCAAAAAACAACGCGGc                                                                                      <  2:2638035/100‑1 (MQ=255)
          cGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGc                                                                                      <  2:2638092/99‑1 (MQ=255)
           gCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCaa                                                                                    <  1:501222/100‑1 (MQ=255)
            ccATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAac                                                                                   <  2:2461777/100‑1 (MQ=255)
             catcatCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAacc                                                                                  <  2:1415800/100‑1 (MQ=255)
                catcTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAaccacc                                                                               >  1:1159762/1‑100 (MQ=255)
                    tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcc                                                                           <  1:528437/100‑1 (MQ=255)
                    tcCATACGGGAGAGTTGCTCTGTAGCGCGCATCAGGCCTATTTCTGCGGTGAGCGCCGAACCAGCACGCTCGTCAAAATAAAACGCGGCAGCCCCCgg                                                                             <  2:1330364/98‑1 (MQ=255)
                      cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcccc                                                                         <  2:787209/100‑1 (MQ=255)
                      cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGcccc                                                                         <  1:538494/100‑1 (MQ=255)
                        tACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAg                                                                       >  1:1099199/1‑100 (MQ=255)
                          cTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGtt                                                                     >  2:2224395/1‑100 (MQ=255)
                          cTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCAGCGGTTCGCGCCGACCCAGCACGCCCGGCAAAAAACCACGCCGCAACCACCGGcccc                                                                         >  2:830567/1‑96 (MQ=255)
                           tGGAGAGTTGCTCTGTAGGGCGCAGCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGt                                                                      <  2:1692923/98‑1 (MQ=255)
                            ggAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAGCCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCa                                                                   <  1:353244/100‑1 (MQ=255)
                            ggAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCa                                                                   <  1:1725354/100‑1 (MQ=255)
                              agagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACg                                                                 <  2:325632/100‑1 (MQ=255)
                              agagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACg                                                                 >  2:1774498/1‑100 (MQ=255)
                               gagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGc                                                                <  1:1193574/100‑1 (MQ=255)
                               gagTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGc                                                                <  2:2251549/100‑1 (MQ=255)
                                agTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCa                                                               <  1:313597/100‑1 (MQ=255)
                                 gTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAg                                                              >  2:2277453/1‑100 (MQ=255)
                                 gTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAg                                                              >  1:1409586/1‑100 (MQ=255)
                                  ttGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGt                                                             <  2:944212/100‑1 (MQ=255)
                                    gCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTa                                                            >  2:1711126/1‑99 (MQ=255)
                                    gCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAg                                                           >  2:2090515/1‑100 (MQ=255)
                                     ctctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGc                                                          <  1:2399618/100‑1 (MQ=255)
                                     ctctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGc                                                          >  1:191307/1‑100 (MQ=255)
                                      tctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCg                                                         <  2:1370452/100‑1 (MQ=255)
                                      tctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCg                                                         <  1:40297/100‑1 (MQ=255)
                                      tctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCAAATCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCg                                                         <  2:1927961/100‑1 (MQ=255)
                                       ctGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGa                                                        >  2:1605967/1‑100 (MQ=255)
                                       ctGTAGCGCGCATCAGGCCGATTCCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGc                                                          <  2:1323692/98‑1 (MQ=255)
                                         gTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATa                                                      <  2:505451/100‑1 (MQ=255)
                                           aGCGCGCATCCGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAAc                                                    >  2:1634441/1‑100 (MQ=255)
                                            gcgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACg                                                   >  2:906207/1‑100 (MQ=255)
                                            gcgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACg                                                   <  1:2503730/100‑1 (MQ=255)
                                            gcgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACg                                                   <  1:1746733/100‑1 (MQ=255)
                                                 cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                              <  1:935147/100‑1 (MQ=255)
                                                 cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCAGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                              <  2:2721686/100‑1 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCAt                                         <  1:2392808/100‑1 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCAt                                         >  1:1757756/1‑100 (MQ=255)
                                                      ggCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCAt                                         >  2:1380314/1‑100 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        >  2:2061643/1‑100 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        >  1:529527/1‑100 (MQ=255)
                                                       gCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        <  1:1474733/100‑1 (MQ=255)
                                                        ccGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                        >  1:1136996/1‑99 (MQ=255)
                                                        ccGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATAc                                       >  1:1202451/1‑100 (MQ=255)
                                                           aTTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCa                                    <  1:635559/100‑1 (MQ=255)
                                                            tttCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAg                                   <  2:483484/100‑1 (MQ=255)
                                                               cTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACt                                <  1:2200372/100‑1 (MQ=255)
                                                                    gTTAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGttt                           <  1:1026174/100‑1 (MQ=255)
                                                                      tAGCGCCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTcc                         <  1:2709216/100‑1 (MQ=255)
                                                                        gcgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGc                       >  2:1778395/1‑100 (MQ=255)
                                                                        gcgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGc                       >  1:2420381/1‑100 (MQ=255)
                                                                        gcgcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGACACCAGCATACCCAGACTGGTTTccac                       >  2:1636864/1‑98 (MQ=255)
                                                                          gcCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCAc                     <  1:612841/100‑1 (MQ=255)
                                                                          gcCGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCAc                     <  1:1362520/100‑1 (MQ=255)
                                                                           ccGAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACt                    >  1:2244776/1‑100 (MQ=255)
                                                                             gAACCAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACtat                  <  2:29726/100‑1 (MQ=255)
                                                                                ccAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAg               <  2:726298/100‑1 (MQ=255)
                                                                                ccAGCACACCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAg               <  2:2649642/100‑1 (MQ=255)
                                                                                 cAGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGt              <  1:2423094/100‑1 (MQ=255)
                                                                                  aGCACGCCCGGCAAAAAAGAACGCGGCAACCCCCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAg               >  2:653953/1‑98 (MQ=255)
                                                                                  aGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTg             <  1:670617/100‑1 (MQ=255)
                                                                                  aGCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTg             >  1:1185108/1‑100 (MQ=255)
                                                                                   gCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTgg            <  2:917972/100‑1 (MQ=255)
                                                                                   gCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTgg            <  1:1605974/100‑1 (MQ=255)
                                                                                   gCACGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTg             >  2:1141579/1‑99 (MQ=255)
                                                                                     aCGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTc          <  2:2486597/100‑1 (MQ=255)
                                                                                     aCGCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTc          <  2:1572066/100‑1 (MQ=255)
                                                                                       gCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAg        <  1:914097/100‑1 (MQ=255)
                                                                                       gCCCGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAg        <  1:1013636/100‑1 (MQ=255)
                                                                                        cccGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGa       >  2:1517443/1‑100 (MQ=255)
                                                                                         ccGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGaa      >  2:1063629/1‑100 (MQ=255)
                                                                                         ccGCCAAAAAACAACGAGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCa         <  2:1201984/97‑1 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     >  1:1492220/1‑100 (MQ=255)
                                                                                          cGGCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc     <  2:1433909/100‑1 (MQ=255)
                                                                                           ggCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAcc    <  1:2090528/100‑1 (MQ=255)
                                                                                           ggCAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAcc    >  1:218317/1‑100 (MQ=255)
                                                                                             cAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAg  >  1:1910547/1‑100 (MQ=255)
                                                                                             cAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAg  <  1:177168/100‑1 (MQ=255)
                                                                                             cAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAg  >  1:1271797/1‑100 (MQ=255)
                                                                                             cAAAAAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGAGAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGCCAGAACCAg  >  2:325741/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: