Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,573,528 A→G 100% pseudogene (1115/1503 nt) yjiT → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,5280AG99.0% 309.1 / ‑5.3 90pseudogene (1115/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (33/56);  total (33/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCA  >  NC_000913/4573432‑4573623
                                                                                                |                                                                                               
aaaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAg                                                                                              <  1:2084211/100‑1 (MQ=255)
 aaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGt                                                                                             <  2:1710264/100‑1 (MQ=255)
 aaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAAGTTGATGGCAg                                                                                              >  1:1705131/1‑99 (MQ=255)
  aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                            >  1:2125186/1‑100 (MQ=255)
  aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                            >  2:927295/1‑100 (MQ=255)
    cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaa                                                                                          >  2:2287942/1‑100 (MQ=255)
      aaGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAAt                                                                                        <  1:2545158/100‑1 (MQ=255)
        gTGAGATTTGGCAGGCAAAATCCATCTGTAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                      <  1:795821/100‑1 (MQ=255)
        gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGGTTGTTGGGAGCATTAAACTTGATGGCAGTGAAAt                                                                                        >  2:411290/1‑98 (MQ=255)
        gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                      >  2:858865/1‑100 (MQ=255)
         tGAGATTTGGCAGGCAAAATCCATCGGAAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGc                                                                                     <  1:128090/100‑1 (MQ=255)
          gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                    <  1:120973/100‑1 (MQ=255)
          gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                    >  1:1315613/1‑100 (MQ=255)
          gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                    <  1:411288/100‑1 (MQ=255)
            gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCa                                                                                   >  2:407391/1‑99 (MQ=255)
            gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                  <  1:1460066/100‑1 (MQ=255)
              tttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTg                                                                                <  1:2277951/100‑1 (MQ=255)
               ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                               <  1:858862/100‑1 (MQ=255)
               ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTg                                                                                <  1:2253542/99‑1 (MQ=255)
                 ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                             >  1:2034664/1‑100 (MQ=255)
                   caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                           >  2:2452018/1‑100 (MQ=255)
                   caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                           <  1:794736/100‑1 (MQ=255)
                    ctgcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                          <  2:1379370/98‑1 (MQ=255)
                    aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                          <  1:4871/100‑1 (MQ=255)
                     ggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTCAACTTGGTGGCAGTGAAATTGCGATTGGt                                                                              >  2:1969388/1‑95 (MQ=255)
                     ggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGt                                                                         <  2:1908083/100‑1 (MQ=255)
                      gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAggag                                                                        >  1:1379362/1‑98 (MQ=255)
                      gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                        >  2:2366001/1‑100 (MQ=255)
                      gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                        >  1:498161/1‑100 (MQ=255)
                       caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGc                                                                       <  1:2533044/100‑1 (MQ=255)
                          aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                    <  2:499953/100‑1 (MQ=255)
                           aTCCATCGGGCCGTTTGACCTTAGTGGCTCGTGCTGGAAGGATGATTGCTGGGCGAATTACACTTCATGGCAGTGAAATAGCCATTggaccag                                                                          <  1:519427/93‑6 (MQ=255)
                            tCCATCGGCAGGTTTTTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTa                                                                  <  2:2553282/100‑1 (MQ=255)
                              cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAc                                                                <  2:1646500/100‑1 (MQ=255)
                               aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               <  1:2590299/100‑1 (MQ=255)
                               aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               >  1:1263887/1‑100 (MQ=255)
                               aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               <  2:1529965/100‑1 (MQ=255)
                                tCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCt                                                              <  2:87102/100‑1 (MQ=255)
                                 cGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCtt                                                             <  2:2152982/100‑1 (MQ=255)
                                    cAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCAt                                                          <  1:1338294/100‑1 (MQ=255)
                                     aGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATc                                                         <  2:990057/100‑1 (MQ=255)
                                      ggTTTGTCGGTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAACCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCg                                                        <  1:1935218/100‑1 (MQ=255)
                                         ttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTg                                                     <  2:2016077/100‑1 (MQ=255)
                                          tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGa                                                    <  2:1951009/100‑1 (MQ=255)
                                          tGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGa                                                    >  1:2422683/1‑100 (MQ=255)
                                             cgttcggtgCTCGTGTTGAAGGGATGACGGTTGGGAGCATTAAACTTGATGGCGGTGAAATTGCCATTCGTGAGGTGACGTTAGCCTTCATCGTTGa                                                    >  1:1591826/9‑97 (MQ=255)
                                              gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                                <  2:2284470/100‑1 (MQ=255)
                                              gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                                <  1:2362614/100‑1 (MQ=255)
                                                tcgTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGa                                              <  1:1566928/98‑1 (MQ=255)
                                                 aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                             <  1:875516/100‑1 (MQ=255)
                                                      cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        <  1:838777/100‑1 (MQ=255)
                                                      cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        <  1:196725/100‑1 (MQ=255)
                                                      cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        >  1:2546547/1‑100 (MQ=255)
                                                        cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        <  1:12631/98‑1 (MQ=255)
                                                          tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                    <  1:2583295/100‑1 (MQ=255)
                                                          tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                    <  1:973205/100‑1 (MQ=255)
                                                            cTGGCGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                  >  2:754621/1‑100 (MQ=255)
                                                            cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAACTTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                  <  1:896190/100‑1 (MQ=255)
                                                            cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                  <  2:1293008/100‑1 (MQ=255)
                                                              ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                                >  2:376075/1‑100 (MQ=255)
                                                              ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                                >  2:1936401/1‑100 (MQ=255)
                                                               gAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                                <  2:1730928/99‑1 (MQ=255)
                                                                 gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                             <  1:1935680/100‑1 (MQ=255)
                                                                  ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTGACCTTCATCGTTGCTGCAGATCAATGGCTGTTACAggg                            >  1:2483950/1‑100 (MQ=255)
                                                                  ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAggg                            >  1:1297385/1‑100 (MQ=255)
                                                                  ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAggg                            >  1:2463728/1‑100 (MQ=255)
                                                                    atgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAc                          <  2:2421075/100‑1 (MQ=255)
                                                                    atgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAc                          >  1:2583025/1‑100 (MQ=255)
                                                                     tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACa                         <  2:2463745/100‑1 (MQ=255)
                                                                       atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAgg                       <  1:1793839/100‑1 (MQ=255)
                                                                        ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGc                      <  1:225211/100‑1 (MQ=255)
                                                                        ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGAGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACCGGc                      >  1:2582948/1‑100 (MQ=255)
                                                                           ttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttc                   <  2:547576/100‑1 (MQ=255)
                                                                            tgGGAGCATTAAACTTGATGGCAGTGAAATTGCTATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttct                  <  2:1598032/100‑1 (MQ=255)
                                                                             gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                 <  2:1266732/100‑1 (MQ=255)
                                                                              ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg                <  2:1343528/100‑1 (MQ=255)
                                                                              ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg                >  2:1282309/1‑100 (MQ=255)
                                                                              ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg                >  2:900827/1‑100 (MQ=255)
                                                                              ggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTg                >  2:2300314/1‑100 (MQ=255)
                                                                               gAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc               <  1:2179170/100‑1 (MQ=255)
                                                                               gAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc               <  2:1765703/100‑1 (MQ=255)
                                                                                 gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg             >  2:62711/1‑100 (MQ=255)
                                                                                  cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGATGCCGTTAACCTTCATCGTTGATGCAGATTAATGGCTATTACAACGACAGGCTTCTTGCAg             >  1:1341488/1‑99 (MQ=255)
                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGTGACAGGCTTCTTGCAgtg           <  2:246421/100‑1 (MQ=255)
                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           <  1:1302167/100‑1 (MQ=255)
                                                                                    ttAAACTTGATGGCAGTGAAATTGCCATTGGGTAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtgt          <  1:2066806/100‑1 (MQ=255)
                                                                                     tAAACTTGATGGCAGTGAACTTGCCATTGGTGAGGTGCCGTTAACCGGCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt            <  1:1280110/97‑1 (MQ=255)
                                                                                     tAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTTCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg         <  2:1815732/100‑1 (MQ=255)
                                                                                       aaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCg       >  1:2556112/1‑100 (MQ=255)
                                                                                        aCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa      >  2:702023/1‑100 (MQ=255)
                                                                                            gatgaCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagca  >  2:2098743/1‑100 (MQ=255)
                                                                                                |                                                                                               
AAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCA  >  NC_000913/4573432‑4573623

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: