Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,163:1 +T 16.1% 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.T16.1% 573.1 / 128.1 186coding (2232/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base . (83/73);  new base T (17/13);  total (100/86)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 193 total aligned reads displayed.

TGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATG  >  NC_000913/3714064‑3714262
                                                                                                    |                                                                                                   
tGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc                                                                                                      >  1:1350040/1‑100 (MQ=255)
  aGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cg                                                                                                   <  2:2927728/100‑1 (MQ=255)
   gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑cgc                                                                                                  >  1:2764229/1‑100 (MQ=255)
     tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGctcgcg                                                                                                 <  1:3053454/100‑6 (MQ=255)
     tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGa                                                                                                >  2:3856601/1‑100 (MQ=255)
       gCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGc                                                                                              <  2:2124499/100‑1 (MQ=255)
        ccGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCt                                                                                             >  2:1296679/1‑100 (MQ=255)
         cGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCt                                                                                             <  2:1044315/99‑1 (MQ=255)
         cGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTg                                                                                            >  1:1373908/1‑100 (MQ=255)
          ggtggtTCATACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                           <  2:2445087/100‑1 (MQ=255)
           gtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTggg                                                                                          <  1:2705261/100‑1 (MQ=255)
           gtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTggg                                                                                          <  1:1851432/100‑1 (MQ=255)
            tggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggg                                                                                         >  1:1087176/1‑100 (MQ=255)
            tggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggg                                                                                         <  1:3206355/100‑1 (MQ=255)
             gttgcacactcTGTAAGTGACAGTTCCGGACCGTTGTACTTTTCCACCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCAAGCTGGGga                                                                                        <  2:2791723/90‑1 (MQ=255)
             ggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTTCAGCCATGCCAGCGCCGTATTACCCGCAAAGCCATTCCCCAGC‑CGCGAGCTGgaga                                                                                        >  2:1179939/1‑100 (MQ=255)
              gtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggggg                                                                                       >  1:616917/1‑98 (MQ=255)
              gtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgag                                                                                       <  2:3720065/100‑1 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaat                                                                                      <  1:549910/100‑3 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgaga                                                                                      >  1:1915358/1‑100 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgaga                                                                                      >  1:487846/1‑100 (MQ=255)
               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCAATCCCCGGC‑CGCGTGCTGGGgag                                                                                       <  2:3167097/99‑1 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGTACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  1:912948/100‑4 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                     <  1:3655504/100‑4 (MQ=255)
                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGAt                                                                                     >  1:1498287/1‑100 (MQ=255)
                 cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATc                                                                                    <  1:3139059/100‑1 (MQ=255)
                    aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatggctg                                                                                 <  2:3628533/100‑8 (MQ=38)
                    aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                 <  1:1743403/100‑1 (MQ=255)
                      gCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTgg                                                                               <  2:1306866/100‑1 (MQ=255)
                      gCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTgg                                                                               <  1:794833/100‑1 (MQ=255)
                       cTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCCGC‑CGCGAGCTGGGGAGATCTTTGGt                                                                              >  2:192594/1‑100 (MQ=255)
                        tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                             <  2:442668/100‑1 (MQ=255)
                        tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                             <  2:2257496/100‑1 (MQ=255)
                         gTAAGTGTCAGTTCCGGACCGTTGTATTCTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTTCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                             <  1:1085715/99‑1 (MQ=255)
                              tgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                       <  1:962061/100‑1 (MQ=255)
                               gtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGt                                                                      >  2:2836687/1‑100 (MQ=255)
                                tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGt                                                                      >  2:2495459/1‑99 (MQ=255)
                                 cAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                    >  1:2473543/1‑100 (MQ=255)
                                  aGTTCCGGACCGTTGTATTTTTTTTTTTTTTTCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                    >  2:3914209/1‑99 (MQ=255)
                                  aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCa                                                                   >  1:1180447/1‑100 (MQ=255)
                                   gTTCCGGACCGTTGTATTTTTCCAGCCATGGCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCGGCACGTTCAc                                                                  >  2:2587198/1‑100 (MQ=255)
                                   gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCAc                                                                  >  2:2911208/1‑100 (MQ=255)
                                     tCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                <  1:1281487/100‑1 (MQ=255)
                                       cGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                              <  1:3301681/100‑1 (MQ=255)
                                        ggACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                             >  2:3501795/1‑100 (MQ=255)
                                           ccGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                          >  1:1224678/1‑100 (MQ=255)
                                             gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                        >  2:220666/1‑100 (MQ=255)
                                             gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                        >  1:3268129/1‑100 (MQ=255)
                                                 tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                    <  2:2055731/100‑1 (MQ=255)
                                                 tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                    >  2:3803689/1‑100 (MQ=255)
                                                 tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                    <  1:466725/100‑1 (MQ=255)
                                                  aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                   >  2:1532906/1‑100 (MQ=255)
                                                  aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                   >  1:1332912/1‑100 (MQ=255)
                                                     tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAt                                                >  1:2906929/1‑100 (MQ=255)
                                                     tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGACCTTTGTTACGAACGTTCACTCCGCCGTTTTTCCa                                                   >  2:1443427/1‑97 (MQ=255)
                                                         cAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCg                                            >  2:2923935/1‑100 (MQ=255)
                                                         cAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCg                                            >  1:3647538/1‑100 (MQ=255)
                                                          aGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGt                                           <  1:2276021/100‑1 (MQ=255)
                                                           gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGTGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTc                                          >  1:3190508/1‑100 (MQ=255)
                                                           gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTc                                          >  2:3362413/1‑100 (MQ=255)
                                                            ccATGCCAGCGCCGTATTACCCGCGCAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCa                                         <  2:839043/100‑1 (MQ=255)
                                                            ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCa                                         <  1:2042829/100‑1 (MQ=255)
                                                                gCCAACGACTTATTACACGAACAGCCATGCCCCCGC‑CGCAAGCTGGGGAGCTCAATGGTCAGCACGTTCACAGACCCGTTTTTACAAATACCGTCAGCgg                                     <  2:2557912/100‑1 (MQ=255)
                                                                 ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGt                                    >  2:70397/1‑100 (MQ=255)
                                                                   aGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                  <  1:1103324/100‑1 (MQ=255)
                                                                    gcgcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                 >  2:98283/1‑100 (MQ=255)
                                                                     cgcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                >  2:808378/1‑100 (MQ=255)
                                                                      gcCGTATTACCCGCACAGCCCTTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAAt                               >  2:2409371/1‑100 (MQ=255)
                                                                         gTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                            <  2:1284644/100‑1 (MQ=255)
                                                                         ccgatctccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                            <  1:2205718/93‑1 (MQ=255)
                                                                           aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                          >  2:1671522/1‑100 (MQ=255)
                                                                              acacGCACAGCCATTCCCCAGC‑AGCGAGGTGAGGAGATCTTTGGTCAGAACGATCACTGCGCCGTTTTTACAAAGACCGTAATGGGTTAAATCCAGCTc                        >  2:3427659/4‑99 (MQ=255)
                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                      >  2:952968/1‑100 (MQ=255)
                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                      >  1:3799413/1‑100 (MQ=255)
                                                                               cccGCACAGCCAGTCCCCAGC‑CGCGAGGTGGGGAGCTCTTTGGTCAGCACGTTCACTGCGCCGTGTTTACTAATACCGTCAGTGGTTAAATCCAGATCCg                      <  2:263093/100‑1 (MQ=255)
                                                                                  gCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGcc                   <  1:1772328/100‑1 (MQ=255)
                                                                                   cacaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                  <  1:1662345/100‑1 (MQ=255)
                                                                                    caaagatctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                 >  1:2177069/9‑100 (MQ=38)
                                                                                     caGCCATTCCCCAGC‑CGCGAGCGGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                  <  2:2873907/98‑1 (MQ=255)
                                                                                           ctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt          >  2:875399/2‑100 (MQ=255)
                                                                                           ctCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCt           >  1:3622136/2‑99 (MQ=255)
                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  1:819311/1‑100 (MQ=255)
                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  2:2777134/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  2:1241172/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:2623388/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  2:2973928/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  1:3025336/1‑100 (MQ=255)
                                                                                             ccccAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTCCGCCGTTGTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  2:3331716/1‑99 (MQ=255)
                                                                                             ccccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc         >  2:3373620/1‑99 (MQ=255)
                                                                                             ccccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg        >  2:60367/1‑100 (MQ=255)
                                                                                              cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt       >  1:1733182/1‑100 (MQ=255)
                                                                                              cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt       >  1:158395/1‑100 (MQ=255)
                                                                                              cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt       >  1:3224681/1‑100 (MQ=255)
                                                                                                cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa     <  1:1671531/100‑1 (MQ=255)
                                                                                                cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa     >  2:493544/1‑100 (MQ=255)
                                                                                                cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa     <  1:372243/100‑1 (MQ=255)
                                                                                                 aGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa    <  2:1541733/100‑1 (MQ=255)
                                                                                                  gC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt   <  2:306599/100‑1 (MQ=255)
                                                                                                  gC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt   <  1:1456102/100‑1 (MQ=255)
                                                                                                   c‑cGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg  <  1:3042307/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑