Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,168 Δ1 bp 7.3% coding (2227/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,1680A. 7.3% 271.0 / 13.7 82coding (2227/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base A (44/32);  new base . (2/4);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714021‑3714316
                                                                                                                                                    |                                                                                                                                                    
gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAg                                                                                                                                                     >  2:64788/1‑149 (MQ=255)
 aTTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGc                                                                                                                                                    >  3:308976/1‑149 (MQ=255)
   tGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTg                                                                                                                                                  <  2:469554/149‑1 (MQ=255)
    gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                                 >  1:694044/1‑149 (MQ=255)
     cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTggg                                                                                                                                                <  3:289810/149‑1 (MQ=255)
      gcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCTGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggg                                                                                                                                               <  3:642258/149‑1 (MQ=255)
       cATGGCAGGATGAGACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGga                                                                                                                                              <  3:63004/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGga                                                                                                                                              >  3:754808/1‑149 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGga                                                                                                                                              <  1:37301/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGCGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGga                                                                                                                                              <  2:389197/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACAAGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGga                                                                                                                                              <  4:542804/149‑1 (MQ=255)
        aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaa                                                                                                                                             <  1:601948/149‑2 (MQ=255)
         tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaat                                                                                                                                            <  1:43218/149‑3 (MQ=255)
          ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAACCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGAt                                                                                                                                           <  4:464241/149‑1 (MQ=255)
            cAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCt                                                                                                                                         >  3:328425/1‑149 (MQ=255)
             aGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCtt                                                                                                                                        >  3:469072/1‑149 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                                       <  1:670530/149‑1 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                                       <  3:88248/149‑1 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTGGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                                       <  2:750530/149‑1 (MQ=255)
                 tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGt                                                                                                                                    <  3:10788/149‑1 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   <  4:389160/149‑1 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   >  2:332601/1‑149 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCAGGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   >  3:525842/1‑149 (MQ=255)
                    aaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTTCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                                 >  4:711408/1‑149 (MQ=255)
                    aaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                                 >  2:661390/1‑149 (MQ=255)
                     aCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  4:629321/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  3:619419/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  4:390684/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  4:71541/1‑149 (MQ=255)
                            tACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGATATCTTTGGTCAGCACGTTCa                                                                                                                         >  4:438169/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCGGCACGTTCACTgcg                                                                                                                    >  3:244759/1‑149 (MQ=255)
                                      gATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                               <  4:308976/149‑1 (MQ=255)
                                        ttATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  2:401642/1‑149 (MQ=255)
                                              gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                       >  4:92979/1‑149 (MQ=255)
                                              gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                       <  4:570730/149‑1 (MQ=255)
                                                      gtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCa                                                                                               <  1:332601/149‑1 (MQ=255)
                                                          ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  2:288431/1‑149 (MQ=255)
                                                          ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  3:140347/1‑149 (MQ=255)
                                                             aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                        >  3:333843/1‑149 (MQ=255)
                                                             aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                        >  3:201337/1‑149 (MQ=255)
                                                              aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  2:698271/1‑149 (MQ=255)
                                                               aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                                      <  1:64788/149‑1 (MQ=255)
                                                                     tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                                >  1:326820/1‑149 (MQ=255)
                                                                     gaaGTGTCAGTTCCGGAGCGTTGTATTTTTCCAGCCAGGCGAGCGCCGTATTCCCCGCACATCCATTCCTCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAGACCGTCAGCGGTTAAATCCAGa                                                                                <  4:54378/148‑1 (MQ=255)
                                                                          gtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCgg                                                                           <  3:153823/149‑1 (MQ=255)
                                                                           tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCACCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                                                                          >  1:710417/1‑149 (MQ=255)
                                                                              gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                                       <  2:566824/149‑1 (MQ=255)
                                                                                 ccGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                        <  1:270854/145‑1 (MQ=255)
                                                                                 ccGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                        >  2:270854/1‑145 (MQ=255)
                                                                                  cGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGccc                                                                   <  3:426103/149‑1 (MQ=255)
                                                                                    gACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             <  4:529128/105‑1 (MQ=255)
                                                                                    gACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  3:529128/1‑105 (MQ=255)
                                                                                       cGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCCCAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  <  2:675613/129‑1 (MQ=255)
                                                                                       cGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  >  1:675613/1‑129 (MQ=255)
                                                                                        gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                             <  2:425120/149‑1 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  2:262090/1‑149 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  1:547262/1‑149 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  3:246961/1‑149 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  3:408313/1‑149 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  2:744701/1‑149 (MQ=255)
                                                                                            tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt                                                         <  3:71541/149‑1 (MQ=255)
                                                                                             aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                        >  2:627240/1‑149 (MQ=255)
                                                                                              tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGc                                                       >  1:749677/1‑149 (MQ=255)
                                                                                                 ttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaaa                                                    <  4:754808/149‑1 (MQ=255)
                                                                                                    cAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCa                                                 <  2:115138/149‑1 (MQ=255)
                                                                                                          tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAg                                           <  4:408313/149‑1 (MQ=255)
                                                                                                          tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAg                                           <  2:352341/149‑1 (MQ=255)
                                                                                                           gCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAgg                                          >  3:2117/1‑149 (MQ=255)
                                                                                                             cAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTggg                                                                                                                                                >  2:230186/1‑45 (MQ=25)
                                                                                                             cAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTggg                                                                                                                                                <  1:230186/45‑1 (MQ=25)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                                >  2:681032/1‑149 (MQ=255)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTACACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                                >  2:632807/1‑149 (MQ=255)
                                                                                                                      aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCt                               <  2:547262/149‑1 (MQ=255)
                                                                                                                       ttACCCGCACAGCCATTCCCCAGC‑CGCAAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATACCCTc                              >  2:755479/1‑149 (MQ=255)
                                                                                                                             gCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                        >  1:204615/1‑149 (MQ=255)
                                                                                                                                caGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAc                     <  3:642858/149‑1 (MQ=255)
                                                                                                                                  gCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCg                   >  2:206390/1‑149 (MQ=255)
                                                                                                                                   ccATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc                  >  2:282140/1‑149 (MQ=255)
                                                                                                                                    cATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCt                 >  3:679212/1‑149 (MQ=255)
                                                                                                                                       tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg              >  2:653501/1‑149 (MQ=255)
                                                                                                                                          ccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa           >  2:39870/1‑149 (MQ=255)
                                                                                                                                           cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCaa          <  1:698271/149‑1 (MQ=255)
                                                                                                                                                    aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  2:694044/149‑1 (MQ=255)
                                                                                                                                                    aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATACCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  3:645331/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714021‑3714316

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑