Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,247,812 C→T 35.9% Q14* (CAG→TAG)  yqjA → general envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,247,8120CT35.9% 115.5 / 107.3 114Q14* (CAG→TAG) yqjAgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (22/19);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

Only 100 of 118 total aligned reads displayed.

TTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGC  >  NC_000913/3247714‑3247909
                                                                                                  |                                                                                                  
ttGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCa                                                                                                   <  2:2682243/100‑1 (MQ=255)
 tGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAg                                                                                                  <  1:402054/100‑1 (MQ=255)
  gCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAgg                                                                                                 <  2:1591546/100‑1 (MQ=255)
      tAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGAttt                                                                                             >  1:1532843/1‑100 (MQ=255)
      tAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGAttt                                                                                             <  2:362041/100‑1 (MQ=255)
       aGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGAtttt                                                                                            <  2:2007394/100‑1 (MQ=255)
        gTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTg                                                                                           >  2:1574053/1‑100 (MQ=255)
        gTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTg                                                                                           <  2:2111583/100‑1 (MQ=255)
         tGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGa                                                                                          <  2:263969/100‑1 (MQ=255)
          gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGaa                                                                                         >  1:275688/1‑100 (MQ=255)
          gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGaa                                                                                         >  2:1793292/1‑100 (MQ=255)
           cATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGaaa                                                                                        <  2:2830650/100‑1 (MQ=255)
            aTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAAc                                                                                       >  2:2612355/1‑100 (MQ=255)
              cTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAAccc                                                                                     <  2:2164043/100‑1 (MQ=255)
              cTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAAccc                                                                                     <  1:1515905/100‑1 (MQ=255)
               ttttAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCt                                                                                    >  1:1113383/1‑100 (MQ=255)
                tttAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTg                                                                                   <  1:1654259/100‑1 (MQ=255)
                  tAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGc                                                                                 <  2:1667550/100‑1 (MQ=255)
                   aaaaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGcc                                                                                <  1:182241/100‑1 (MQ=255)
                   aaaaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGcc                                                                                <  2:1854814/100‑1 (MQ=255)
                    aaaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCa                                                                               >  1:431832/1‑100 (MQ=255)
                     aaCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCaa                                                                              >  1:792556/1‑100 (MQ=255)
                       ccAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATc                                                                            >  2:731900/1‑100 (MQ=255)
                       ccAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATc                                                                            <  1:423064/100‑1 (MQ=255)
                           cGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                        >  2:2222043/1‑100 (MQ=255)
                           cGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                        <  2:1294412/100‑1 (MQ=255)
                           cGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                        <  2:1066701/100‑1 (MQ=255)
                           cGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATc                                                                        <  2:1692227/100‑1 (MQ=255)
                            gTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCg                                                                       <  2:759689/100‑1 (MQ=255)
                              ggCGTTAGCCGGTTCACCAGGAATAATGAATGGAACTTTAGACCCAATTGCTGCAAGCCCTGTGGCCGCAGGACGTTGAAACCCTGGCCAATCCATCgat                                                                     >  1:1525113/1‑100 (MQ=255)
                                 gTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCATCgatggat                                                                <  1:260128/100‑4 (MQ=255)
                                  ttAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGAtt                                                                <  1:2139093/100‑1 (MQ=255)
                                   tAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTg                                                               >  2:182447/1‑100 (MQ=255)
                                   tAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTg                                                               <  1:371880/100‑1 (MQ=255)
                                    aaCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTgg                                                              >  1:1035363/1‑100 (MQ=255)
                                      ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                            >  2:197079/1‑100 (MQ=255)
                                      ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                            >  2:600441/1‑100 (MQ=255)
                                      ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCa                                                            <  2:893681/100‑1 (MQ=255)
                                        gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTTGGCGCAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                          <  1:695317/100‑4 (MQ=255)
                                        gATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccag                                                          <  2:38055/100‑4 (MQ=255)
                                        gATGCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCAt                                                           <  1:2648678/99‑1 (MQ=255)
                                         aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgt                                                         <  2:205741/100‑1 (MQ=255)
                                         aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCATCGATGGATTGGccagg                                                         <  2:1378470/100‑5 (MQ=255)
                                           tCACCAGGAATAATGAATGCAACTTTTGACCCAATTGCTGCAAGCCTTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCAGCGATCATTGGCAggttg                                                       <  1:2728875/100‑1 (MQ=255)
                                            cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCAtgttgt                                                      >  1:2858241/1‑100 (MQ=255)
                                               cAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtt                                                   <  2:710959/100‑1 (MQ=255)
                                               cAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtt                                                   >  2:657015/1‑100 (MQ=255)
                                                 ggAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTAtttt                                                 >  1:1597168/1‑100 (MQ=255)
                                                  gAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTg                                                <  2:1423190/100‑1 (MQ=255)
                                                   aataatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGt                                               <  2:2897526/100‑1 (MQ=255)
                                                     taatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCt                                             <  1:2495564/100‑1 (MQ=255)
                                                     taatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCt                                             <  2:2252897/100‑1 (MQ=255)
                                                      aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCtt                                            <  1:1651003/100‑1 (MQ=255)
                                                      aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCtt                                            >  2:623180/1‑100 (MQ=255)
                                                         gaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtt                                         <  2:1446677/100‑1 (MQ=255)
                                                         gaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtt                                         <  2:1891414/100‑1 (MQ=255)
                                                         gaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtt                                         <  2:1744192/100‑1 (MQ=255)
                                                          aatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgttt                                        >  1:1654548/1‑100 (MQ=255)
                                                            tgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCttgtttgt                                      >  2:1431486/1‑100 (MQ=255)
                                                              gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTaa                                    >  1:1973424/1‑100 (MQ=255)
                                                              gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGAATTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTaa                                    >  1:1591259/1‑100 (MQ=255)
                                                               aaCTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                   <  1:1869419/100‑1 (MQ=255)
                                                               aaCTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                   >  1:858382/1‑100 (MQ=255)
                                                                      gACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTc                            >  1:707770/1‑100 (MQ=255)
                                                                        cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                          >  2:434123/1‑100 (MQ=255)
                                                                        cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                          >  2:1315790/1‑100 (MQ=255)
                                                                        cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                          >  1:2108590/1‑100 (MQ=255)
                                                                        cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                          >  2:658436/1‑100 (MQ=255)
                                                                         ccAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCtt                         <  2:543885/100‑1 (MQ=255)
                                                                          cAATTGCTGCAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGCTGATTGGCAGGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                        >  2:75312/1‑100 (MQ=255)
                                                                          cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                        <  1:1017036/100‑1 (MQ=255)
                                                                            aTTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaa                      >  1:478218/1‑100 (MQ=255)
                                                                            aTTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAACCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaa                      <  1:1072688/100‑1 (MQ=255)
                                                                               gctgcAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                   >  2:1759633/1‑100 (MQ=255)
                                                                               gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                   <  2:2032695/100‑1 (MQ=255)
                                                                               gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                   >  2:1823484/1‑100 (MQ=255)
                                                                               gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                   >  1:113097/1‑100 (MQ=255)
                                                                                 tgcAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                 >  2:172074/1‑100 (MQ=255)
                                                                                  gcAAGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                >  1:2456147/1‑100 (MQ=255)
                                                                                  gcAAGCCCTGTGGGCGTAGGATTCTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTCTTTTGAATTGTTTGTCATTTTGTTGCTTGAAAACGgc                >  1:1983484/1‑100 (MQ=255)
                                                                                  gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgc                >  2:2251269/1‑100 (MQ=255)
                                                                                   cAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgct               >  1:839797/1‑100 (MQ=255)
                                                                                   cAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgct               >  1:2539619/1‑100 (MQ=255)
                                                                                     aGCCCTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttg             >  2:171708/1‑100 (MQ=255)
                                                                                     aGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttg             >  1:1071562/1‑100 (MQ=255)
                                                                                      gCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgc            >  2:212646/1‑100 (MQ=255)
                                                                                      gCCCTGTGGGCGCAGGATTTTGAAACCCTGACCAATTCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgc            <  1:2696050/100‑1 (MQ=255)
                                                                                      gCCCTGTGGGCGCAGGATTTTGAAACCCTCGCCACTCCATCGATGATTCGCATGTTGTATTTCGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgc            >  1:978985/1‑100 (MQ=255)
                                                                                         cTGTGGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTc         >  2:1197631/1‑100 (MQ=255)
                                                                                         cTGTGGGCGCAGGATCTGAAAACACTGGTTACTCTATCGATGATTGGCATGTTGTATTTTGTTTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgctt          <  1:2231427/99‑1 (MQ=255)
                                                                                          tgtgGGCGCAGGATTTTGAAACCCTGACCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTcc        <  1:2536533/100‑1 (MQ=255)
                                                                                           gtgGGCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcg       >  2:1467215/1‑100 (MQ=255)
                                                                                           gtgGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcg       >  1:1822935/1‑100 (MQ=255)
                                                                                            tgGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg      <  2:389848/100‑1 (MQ=255)
                                                                                              ggCGTAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcg    >  1:2707833/1‑100 (MQ=255)
                                                                                              ggCGCAGGATTTTGAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg   <  2:2246760/100‑1 (MQ=255)
                                                                                               gcgtAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg   <  2:32414/100‑1 (MQ=255)
                                                                                               gcgtAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcgg   >  2:2561261/1‑100 (MQ=255)
                                                                                                cgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc  <  2:286883/100‑1 (MQ=255)
                                                                                                cgcAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc  >  1:1645668/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGATGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGC  >  NC_000913/3247714‑3247909

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: