Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,150 T→A 16.4% intergenic (+38/‑170) ihfB → / → ycaI integration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1500TA16.4% 308.6 / 44.6 122intergenic (+38/‑170)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/61);  new base A (9/11);  total (50/72)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 123 total aligned reads displayed.

ATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGC  >  NC_000913/964052‑964245
                                                                                                  |                                                                                               
ntACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGt                                                                                                 >  2:1194401/2‑99 (MQ=255)
aTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTAATCAAACTTGAACAAGAGAAAAGCACCTGTc                                                                                                >  2:2393389/1‑100 (MQ=255)
  aCGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                                              <  1:392267/100‑1 (MQ=255)
   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                             >  2:2405845/1‑100 (MQ=255)
   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                             <  1:2408803/100‑1 (MQ=255)
   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                             >  1:1339361/1‑100 (MQ=255)
   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                             >  1:1350495/1‑100 (MQ=255)
     ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTg                                                                                           <  1:2522592/100‑1 (MQ=255)
      tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCACACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                          >  2:1101467/1‑100 (MQ=255)
      tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                          <  2:1108429/100‑1 (MQ=255)
      tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                          >  2:141337/1‑100 (MQ=255)
        cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtt                                                                                        >  1:1140371/1‑100 (MQ=255)
         tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttt                                                                                       >  1:51897/1‑100 (MQ=255)
         tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttt                                                                                       <  2:2505135/100‑1 (MQ=255)
          cACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtttt                                                                                      <  2:1651559/100‑1 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                     <  2:2483608/100‑1 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                     <  1:1744013/100‑1 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                     <  1:494380/100‑2 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                     <  1:1960166/100‑2 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                     <  1:2235356/100‑2 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                     <  2:827370/100‑2 (MQ=255)
           aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                     <  2:1882925/100‑2 (MQ=255)
            cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                    >  1:158376/1‑100 (MQ=255)
             tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTc                                                                                   <  1:2509400/100‑1 (MQ=255)
              ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCa                                                                                  <  2:863533/100‑1 (MQ=255)
               tAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCa                                                                                  <  1:204862/99‑1 (MQ=255)
               tAAACCAGGTAAAGAACTGCGCGATCGCGCCAATATTGACGGTTAAGTTGTTTACTCAACCTTGAACGAGAGAAAAGCACCCGAAAGGTGCTTTtctctc                                                                                 <  1:1081929/100‑3 (MQ=255)
                aaaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                <  2:1041036/100‑1 (MQ=255)
                 aaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                               <  1:2390271/100‑1 (MQ=255)
                   ccTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtct                                                                             >  2:161353/1‑100 (MQ=255)
                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                            <  2:132052/100‑1 (MQ=255)
                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                            <  2:51898/100‑1 (MQ=255)
                     tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                           >  1:1450636/1‑100 (MQ=255)
                     tGGTAAAGAACTGCGCGATCGCGACAATATTTACGGTTAAGTTTTTTACTCAATCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                           <  1:1451087/100‑1 (MQ=255)
                      ggTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTa                                                                          >  1:8398/1‑100 (MQ=255)
                       gTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTaa                                                                         >  1:227799/1‑100 (MQ=255)
                         aaaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                       <  2:2042681/100‑1 (MQ=255)
                         aaaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                       <  2:1533309/100‑1 (MQ=255)
                          aaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                      >  2:737776/1‑100 (MQ=255)
                           aGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                                     >  1:66980/1‑100 (MQ=255)
                           aGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                                     >  1:1351359/1‑100 (MQ=255)
                            gAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgg                                                                    <  2:1471991/100‑1 (MQ=255)
                             aaCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGTTTTTTTCATTTCTCTAATctgg                                                                    <  1:1004922/99‑1 (MQ=255)
                               cTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                 <  2:985197/100‑1 (MQ=255)
                                tgtgcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaact                                                                <  2:1569035/97‑1 (MQ=255)
                                    cgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgaaa                                                            >  2:840564/1‑100 (MQ=255)
                                         gcgcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTg                                                        <  2:2055730/99‑1 (MQ=255)
                                                 atTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCaga                                               >  2:867930/1‑100 (MQ=255)
                                                    tACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtt                                            <  2:515347/100‑1 (MQ=255)
                                                    tACGGCTAAGATTTTTACTCAAACTTCAACGAGAGACAAGCACCTGTCGAGTGCGTTTTTGATTTCTCTAATCTTGAACTGGAAGCAGCTTCGcggcatt                                            >  1:1542204/1‑94 (MQ=255)
                                                      cGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                          <  1:2163965/100‑1 (MQ=255)
                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                         >  2:2178567/1‑100 (MQ=255)
                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                         <  2:2571199/100‑1 (MQ=255)
                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                         >  1:306665/1‑100 (MQ=255)
                                                        gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGa                                        <  2:1685720/100‑1 (MQ=255)
                                                        gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGa                                        <  2:1971167/100‑1 (MQ=255)
                                                         ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGaa                                       >  2:1413067/1‑100 (MQ=255)
                                                         ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGaa                                       <  2:684974/100‑1 (MQ=255)
                                                           aaGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACa                                     >  2:1204847/1‑100 (MQ=255)
                                                            aGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAg                                    <  1:361384/100‑1 (MQ=255)
                                                             gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGt                                   >  2:1333605/1‑100 (MQ=255)
                                                             gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGt                                   >  2:2298828/1‑100 (MQ=255)
                                                             gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGt                                   >  2:2266005/1‑100 (MQ=255)
                                                                 tttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTc                               <  2:1677877/100‑1 (MQ=255)
                                                                  ttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCa                              <  2:1205695/100‑1 (MQ=255)
                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGTACTGGAAGTTGCCTCGCAGAGTGTTTAACAGTTTTCAcc                            >  2:765176/1‑100 (MQ=255)
                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGTCTCGCAGAGTTTTGAACAGTTTTCAcc                            <  1:346267/100‑1 (MQ=255)
                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                            >  1:1909442/1‑100 (MQ=255)
                                                                     cTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAccc                           >  1:833556/1‑100 (MQ=255)
                                                                     cTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAccc                           <  1:827740/100‑1 (MQ=255)
                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCtt                         >  1:2395302/1‑100 (MQ=255)
                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCtt                         >  1:29283/1‑100 (MQ=255)
                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCtt                         >  2:2229783/1‑100 (MQ=255)
                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCtt                         <  2:227801/100‑1 (MQ=255)
                                                                          aCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                      <  2:1693610/100‑1 (MQ=255)
                                                                           cTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGt                     <  2:185499/100‑1 (MQ=255)
                                                                           cTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGt                     <  1:2393374/100‑1 (MQ=255)
                                                                              gAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaa                  >  2:1729789/1‑100 (MQ=255)
                                                                                tcttgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAt                <  1:2124662/93‑1 (MQ=38)
                                                                                aCGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAt                <  1:1080047/100‑1 (MQ=255)
                                                                                 cGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtt               >  2:2378536/1‑100 (MQ=255)
                                                                                  atgaaaaaaGCACCCGACAGGTGCTTTTCTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGCACAGTTTTCACCCTTTCGTTAAAttc              >  1:422932/6‑100 (MQ=37)
                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttct             >  1:4708/5‑100 (MQ=38)
                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtaga             >  2:2124676/5‑97 (MQ=38)
                                                                                   agagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttct             >  1:1404901/1‑100 (MQ=255)
                                                                                    gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt            >  2:382696/1‑100 (MQ=255)
                                                                                    gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAattatt            >  2:604207/4‑100 (MQ=38)
                                                                                    gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt            >  1:1310180/4‑100 (MQ=39)
                                                                                    gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt            >  1:2470170/4‑100 (MQ=39)
                                                                                    gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt            >  1:854867/4‑100 (MQ=39)
                                                                                     agaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc           <  2:594727/100‑1 (MQ=255)
                                                                                     agaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc           <  2:1951258/100‑1 (MQ=255)
                                                                                      gaAAAGCACCTGTCGGGTGCTTTTTTTATTTCTTTAATCTGGAACTGGAAGCTAACCCGCAGTGTCTTGAACAGAATTCACCCTTTCGTTAAAttctt            >  1:1220150/1‑98 (MQ=255)
                                                                                      gaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct          >  2:2040487/1‑100 (MQ=255)
                                                                                       aaaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGTTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg         <  2:204689/100‑1 (MQ=255)
                                                                                          aGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAat      <  2:1439457/100‑1 (MQ=255)
                                                                                           gCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAata     <  2:194976/100‑1 (MQ=255)
                                                                                             aCCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTGTGAACAGTTTTCACCCTTTCGTGAAATTCTTCTGAATATg   <  1:819317/100‑1 (MQ=255)
                                                                                              ccTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATg   <  1:665158/99‑1 (MQ=255)
                                                                                              ccTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc  <  2:2581463/100‑1 (MQ=255)
                                                                                                  |                                                                                               
ATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGC  >  NC_000913/964052‑964245

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: