Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,163:1 +T 17.3% coding (2232/2334 nt) bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,1631.T17.3% 684.2 / 174.8 231coding (2232/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base . (97/94);  new base T (20/20);  total (117/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 247 total aligned reads displayed.

TGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATG  >  NC_000913/3714064‑3714262
                                                                                                    |                                                                                                   
 gAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc‑c                                                                                                    <  2:4067288/100‑1 (MQ=255)
  aGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cg                                                                                                   >  2:127076/1‑100 (MQ=255)
  aGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cg                                                                                                   >  1:615437/1‑100 (MQ=255)
  aGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTAACCGCACGGCCATTCCCCAGC‑cg                                                                                                   <  2:2254409/100‑1 (MQ=255)
   gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cgc                                                                                                  >  2:4111631/1‑100 (MQ=255)
    cTGGCCGGTGGTTCAAACGCCGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cgcg                                                                                                 >  2:2566633/1‑100 (MQ=255)
      ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAg                                                                                               >  1:3271893/1‑100 (MQ=255)
       gCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGc                                                                                              <  2:1455175/100‑1 (MQ=255)
          ggtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                           <  1:2916166/100‑1 (MQ=255)
            tggtTCAAACGCTGTAATTGTCATTTCCTTACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggg                                                                                         <  1:1895986/100‑1 (MQ=255)
            tggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGCGCTgggg                                                                                         >  2:2631937/1‑100 (MQ=255)
            tggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggg                                                                                         <  2:3150950/100‑1 (MQ=255)
             ggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGga                                                                                        <  2:2203371/100‑1 (MQ=255)
              gtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgag                                                                                       <  1:1146008/100‑1 (MQ=255)
              gtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgag                                                                                       >  1:3444832/1‑100 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCCTTCCCCAGC‑CGCGAGCTGGGgaga                                                                                      >  2:1551861/1‑100 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgaga                                                                                      <  2:2007541/100‑1 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgaga                                                                                      <  2:30871/100‑1 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGcaga                                                                                      >  2:1165504/1‑100 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  1:1704268/100‑4 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  2:2858929/100‑4 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  1:144662/100‑4 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  2:495862/100‑4 (MQ=255)
                  aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCCGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGGTCt                                                                                   >  1:2495778/1‑100 (MQ=255)
                  aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatggc                                                                                   <  1:2389476/100‑6 (MQ=38)
                  aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCt                                                                                   >  1:3477555/1‑100 (MQ=255)
                  aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCt                                                                                   <  1:1055684/100‑1 (MQ=255)
                  aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCAGGCCAGCGCTGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  2:780241/98‑4 (MQ=255)
                   aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGCGATCtt                                                                                  >  2:1932380/1‑100 (MQ=255)
                         gTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCa                                                                            <  1:1010988/100‑1 (MQ=255)
                            aGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCa                                                                         <  1:4085122/100‑1 (MQ=255)
                             gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                        <  2:1426710/100‑1 (MQ=255)
                             gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                        <  2:2135281/100‑1 (MQ=255)
                                 cAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGCGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                    >  2:3926923/1‑100 (MQ=255)
                                 cAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                    >  2:3369869/1‑100 (MQ=255)
                                 cAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                    >  1:433377/1‑100 (MQ=255)
                                  aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCCCAGCCCTTCCCCAGC‑CGCGAGCTGGGGAGCTCTTTGGTCAGCACGt                                                                      >  2:3907022/1‑97 (MQ=255)
                                  aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCa                                                                   >  2:3175065/1‑100 (MQ=255)
                                       cGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                              >  1:1929776/1‑100 (MQ=255)
                                         gACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                            <  1:3101467/100‑1 (MQ=255)
                                          aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCg                                                           <  1:2434084/100‑1 (MQ=255)
                                             gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGTCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                        >  2:2892931/1‑100 (MQ=255)
                                             gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                        >  1:113280/1‑100 (MQ=255)
                                              ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                       <  1:3518688/100‑1 (MQ=255)
                                               tGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttttt                                                      >  1:2666451/1‑100 (MQ=255)
                                                 tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                    <  2:347452/100‑1 (MQ=255)
                                                  aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                    <  1:3899224/99‑1 (MQ=255)
                                                  aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                   <  2:545341/100‑1 (MQ=255)
                                                  aTTTTTCCAGCCATGCAAGCGCCGTATTACCCGCAAAGCCATTCCACAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTAACTGCGCCGTTTTTAc                                                    >  1:2089468/1‑99 (MQ=255)
                                                    ttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCATTTTTACaaa                                                 >  1:3797486/1‑100 (MQ=255)
                                                     tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAt                                                >  1:2405716/1‑100 (MQ=255)
                                                       tcaaGCCATGCGAGCGCCGTAGTACCCGCACATCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                   <  2:3986329/92‑1 (MQ=255)
                                                        ccAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTCTGGTCTGCACGTTCACTGCGCCGTTTTTACAAATAcc                                             >  2:2821418/1‑100 (MQ=255)
                                                         cAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCg                                            >  1:2796399/1‑100 (MQ=255)
                                                           gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTCTTACAAATACCGTc                                          >  1:3069404/1‑100 (MQ=255)
                                                             cATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAg                                        <  2:3418353/100‑1 (MQ=255)
                                                              aTGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGc                                       >  1:395037/1‑100 (MQ=255)
                                                              aTGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGc                                       >  1:1531636/1‑100 (MQ=255)
                                                               tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCg                                      <  1:319563/100‑1 (MQ=255)
                                                                 ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGt                                    >  1:582408/1‑100 (MQ=255)
                                                                  cAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                   <  2:1659969/100‑1 (MQ=255)
                                                                  cAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                   >  1:2518614/1‑100 (MQ=255)
                                                                  cAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                   >  1:1669554/1‑100 (MQ=255)
                                                                   aGCGCCGTATTACCCGCACAGCCATTCCCCAGT‑CGCGAGGTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                  <  2:2939374/100‑1 (MQ=255)
                                                                       ccGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                              <  1:362938/100‑1 (MQ=255)
                                                                       ccGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                              >  2:2978931/1‑100 (MQ=255)
                                                                        cGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATcc                             <  1:3758330/100‑1 (MQ=255)
                                                                           aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                          <  1:758691/100‑1 (MQ=255)
                                                                           aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                          <  2:1803170/100‑1 (MQ=255)
                                                                            ttACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                         <  1:4158060/100‑1 (MQ=255)
                                                                             tACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTTGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATc                        <  2:1037078/100‑1 (MQ=255)
                                                                             tACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATc                        >  2:1349617/1‑100 (MQ=255)
                                                                              aCCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATcc                       >  1:2197745/1‑100 (MQ=255)
                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATGTTTGGTCAGGACGGTCACTGCGCGGTGTTTACAAATACAGTAAGCGGGTAAAGCCa                            >  2:3790371/1‑94 (MQ=255)
                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                      <  1:922283/100‑1 (MQ=255)
                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGCTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                      <  1:2953258/100‑1 (MQ=255)
                                                                                ccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTCAATCCAGATccag                     <  1:2840052/100‑3 (MQ=255)
                                                                                 cGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                    >  2:2296153/1‑100 (MQ=255)
                                                                                 cGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                    >  1:2701664/1‑100 (MQ=255)
                                                                                   cacaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCaagatcggaag                  >  2:3594655/1‑90 (MQ=255)
                                                                                   cacaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                  <  2:3570661/100‑1 (MQ=255)
                                                                                       gCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcc              >  1:1963222/1‑100 (MQ=255)
                                                                                       gCCATGCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcc              <  2:2399526/100‑1 (MQ=255)
                                                                                        gatctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGccc             >  2:93770/5‑100 (MQ=255)
                                                                                           ctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt          >  2:2482443/2‑100 (MQ=255)
                                                                                           ctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt          >  1:986629/2‑100 (MQ=255)
                                                                                           ctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt          >  1:1257856/2‑100 (MQ=255)
                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  2:3260518/1‑100 (MQ=255)
                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  1:3335079/1‑100 (MQ=255)
                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAACACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCGATGCCCCTTc         >  2:2762725/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:3921861/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:3786807/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:3677063/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:2622251/1‑100 (MQ=255)
                                                                                             ccccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        <  1:3578018/100‑1 (MQ=255)
                                                                                              cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt       <  1:2296158/100‑1 (MQ=255)
                                                                                               ccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt       >  1:862377/1‑99 (MQ=255)
                                                                                               ccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg      >  2:1120252/1‑100 (MQ=255)
                                                                                                cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa     >  2:4013976/1‑100 (MQ=255)
                                                                                                   c‑cGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg  >  2:306389/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATG  >  NC_000913/3714064‑3714262

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: ‑