Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 17.5% intergenic (+57/+27) chiQ → / ← fur chitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,1730AG17.5% 353.7 / 69.5 166intergenic (+57/+27)chiQ/furchitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (71/66);  new base G (14/15);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 166 total aligned reads displayed.

GTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGAC  >  NC_000913/710076‑710270
                                                                                                 |                                                                                                 
gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAgg                                                                                                 <  2:1937234/100‑1 (MQ=255)
 tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGt                                                                                                >  1:2721795/1‑100 (MQ=255)
 tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGt                                                                                                >  2:16508/1‑100 (MQ=255)
 tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCCCAgg                                                                                                 >  2:3274920/1‑99 (MQ=255)
  gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGtt                                                                                               >  1:200221/1‑100 (MQ=255)
  gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGtt                                                                                               >  2:387023/1‑100 (MQ=255)
  gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGGAACTCACGCTATTCATATTTCCCGTGCATAAAAAAGCCAACCTGCAGGtt                                                                                               >  1:3466384/1‑100 (MQ=255)
        gTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttt                                                                                         <  1:1541210/100‑1 (MQ=255)
        gTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttt                                                                                         <  1:373851/100‑1 (MQ=255)
        gTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCttt                                                                                         <  2:300543/100‑1 (MQ=255)
          cTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttt                                                                                       <  2:3809889/100‑2 (MQ=255)
          cTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                       <  2:296879/100‑1 (MQ=255)
          cTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                       >  1:2049352/1‑100 (MQ=255)
           tGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtct                                                                                      <  2:2385315/100‑1 (MQ=255)
            ggCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCtttttta                                                                                     <  2:2051376/100‑4 (MQ=255)
            ggCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCtttttta                                                                                     <  1:1968038/100‑4 (MQ=39)
             gCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCg                                                                                    >  2:2706621/1‑100 (MQ=255)
              cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttatg                                                                                   <  1:1052321/100‑6 (MQ=38)
              cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttatg                                                                                   <  1:1981328/100‑6 (MQ=38)
              cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttatg                                                                                   <  1:1633252/100‑6 (MQ=38)
              cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCg                                                                                    >  1:356111/1‑99 (MQ=255)
              cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                   <  2:2499609/100‑1 (MQ=255)
               aGGAAATACGCAGTAATAACAAATAAGTTAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGtt                                                                                  <  1:1624711/100‑1 (MQ=255)
                ggAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTc                                                                                 <  1:2684879/100‑1 (MQ=255)
                 gAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCa                                                                                <  1:631896/100‑1 (MQ=255)
                 gAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCa                                                                                <  2:3412088/100‑1 (MQ=255)
                  aaaTACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAg                                                                               <  2:1621406/100‑1 (MQ=255)
                   aaTACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAg                                                                               <  1:33172/99‑1 (MQ=255)
                    aTACGCAGTAATACCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggc                                                                             <  1:3337246/100‑1 (MQ=255)
                    aTACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggc                                                                             <  2:3455080/100‑1 (MQ=255)
                       cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctgg                                                                          <  2:1579012/100‑1 (MQ=255)
                       cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctgg                                                                          >  1:3415367/1‑100 (MQ=255)
                        gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctggc                                                                         <  2:3165256/100‑1 (MQ=255)
                         cAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctggct                                                                        >  2:3747889/1‑100 (MQ=255)
                               taaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                  <  1:387018/100‑1 (MQ=255)
                               taaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                  >  2:288369/1‑100 (MQ=255)
                               taaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                  >  2:1735663/1‑100 (MQ=255)
                                   aaaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCtt                                                              >  1:2992894/1‑100 (MQ=255)
                                         gTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgcg                                                        <  2:882263/100‑1 (MQ=255)
                                          tGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgcgc                                                       >  2:406407/1‑100 (MQ=255)
                                           gagaGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCa                                                      >  2:1886461/1‑100 (MQ=255)
                                               gCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGtt                                                  >  2:3503973/1‑100 (MQ=255)
                                                cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTc                                                 >  1:1044478/1‑100 (MQ=255)
                                                 tGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCa                                                <  2:471887/100‑1 (MQ=255)
                                                 tGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCa                                                >  2:889921/1‑100 (MQ=255)
                                                 tGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCa                                                <  2:768551/100‑1 (MQ=255)
                                                  gTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                               >  2:1275170/1‑100 (MQ=255)
                                                  gTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                               <  1:1804123/100‑1 (MQ=255)
                                                  gTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTAGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                               >  1:83190/1‑100 (MQ=255)
                                                   tAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATc                                              >  2:2121141/1‑100 (MQ=255)
                                                     nctctCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCtt                                            >  2:3148150/2‑100 (MQ=255)
                                                     aCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCtt                                            <  1:1733519/100‑1 (MQ=255)
                                                          cGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                       <  2:3559042/100‑1 (MQ=255)
                                                          cGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                       <  2:3411605/100‑1 (MQ=255)
                                                          cGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                       <  1:1501006/100‑1 (MQ=255)
                                                          cGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                       <  1:193988/100‑1 (MQ=255)
                                                          cGCTTTTATTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATTTCATCTTCGCGGc                                      <  1:557270/100‑1 (MQ=255)
                                                             ttttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCaa                                    >  2:2182337/1‑100 (MQ=255)
                                                              tttCTTATTTCCCTTGCATCACAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGc                                      <  1:2272576/97‑1 (MQ=255)
                                                              tttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTTCGCATGTTCATCTTCGCGGCaa                                    >  1:1256511/1‑99 (MQ=255)
                                                                tCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCg                                 >  1:1761371/1‑100 (MQ=255)
                                                                 cTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGc                                <  1:3652947/100‑1 (MQ=255)
                                                                   tATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCt                              >  2:132995/1‑100 (MQ=255)
                                                                     tttCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTc                            <  1:1066508/100‑1 (MQ=255)
                                                                      ttCGCTTGTATAAAAAAGACAACCCGCAGTTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCTTTCGTGTGCATGTTAATCTTCGCGGCAATCGCCTTCg                           <  1:1410590/100‑1 (MQ=255)
                                                                       tCCCTTGCATAAAAAAGCCAACCCGCAGGTTGTCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCgg                          <  1:6513/100‑1 (MQ=255)
                                                                       tCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCgg                          <  1:2328159/100‑1 (MQ=255)
                                                                        cccTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGc                         <  2:1891002/100‑1 (MQ=255)
                                                                         cctgaacgagaaaaGCCAACCTGCGGGTTGGCTTTTTTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGca                        >  2:2016594/11‑100 (MQ=25)
                                                                          cTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcac                       >  2:3503307/1‑100 (MQ=255)
                                                                           tgaacgagaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcaca                      >  1:3614906/9‑100 (MQ=37)
                                                                              acgagaaaaGCCAACCTGCGGTTTGGCTTTTCTCGGTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTg                   >  1:1509074/6‑100 (MQ=38)
                                                                               cgagaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGa                  >  2:3128184/5‑100 (MQ=38)
                                                                               aTAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGa                  >  2:2974468/1‑100 (MQ=255)
                                                                                tAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAc                 >  1:3373974/1‑100 (MQ=255)
                                                                                gagaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAc                 >  2:1456440/4‑100 (MQ=39)
                                                                                 agaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAcc                >  2:3538131/3‑100 (MQ=255)
                                                                                 aaaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAcc                >  2:2592534/1‑100 (MQ=255)
                                                                                  gaaaaGCCAAGCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTGGCCTTCGTGCGCATGTTCAGCTTCGCGGCAATCGCCTTCGGCACAGTGACCg               >  1:2617513/2‑100 (MQ=255)
                                                                                  aaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCg               >  2:3185696/1‑100 (MQ=255)
                                                                                  aaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCg               >  1:416530/1‑100 (MQ=255)
                                                                                   aaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGt              >  2:3201897/1‑100 (MQ=255)
                                                                                   aaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGt              >  2:576101/1‑100 (MQ=255)
                                                                                   aaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGt              >  1:37768/1‑100 (MQ=255)
                                                                                    aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTGAGGCGGGTTTATTTGCGTTGGTGCGCATGTTCATATTCGGGGCAATCGCCTTCGGCACATTGACCGt              >  2:521610/1‑99 (MQ=25)
                                                                                    aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTa             >  2:924680/1‑100 (MQ=255)
                                                                                    aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTa             >  2:176651/1‑100 (MQ=255)
                                                                                    aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGACTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTCCGGCACAGTGCCCg               >  1:600138/1‑98 (MQ=38)
                                                                                       gCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAg          >  2:1827720/1‑100 (MQ=255)
                                                                                        ccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGGCTTCGTGCGCGTGTTCATCTTCGCGGCAATCGCGTTGGGCACAGGGGGCGTTAAGa         >  1:3147807/1‑100 (MQ=255)
                                                                                         cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCCATCGCCTTCGGCACAGTGACCGTaa            >  1:3245073/1‑96 (MQ=255)
                                                                                         cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGAt        <  2:3581439/100‑1 (MQ=255)
                                                                                           aCCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATag      <  2:83192/100‑1 (MQ=255)
                                                                                            cccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga     <  2:2035118/100‑1 (MQ=255)
                                                                                            cccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga     <  2:1539928/100‑1 (MQ=255)
                                                                                            cccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga     >  1:3018309/1‑100 (MQ=255)
                                                                                            cccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga     >  1:193759/1‑100 (MQ=255)
                                                                                             ccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATagag    <  2:671678/100‑1 (MQ=255)
                                                                                              cGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATagag    <  1:2413733/99‑1 (MQ=255)
                                                                                               gCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGAc  <  1:1886452/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
GTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGAC  >  NC_000913/710076‑710270

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: