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

GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714021‑3714316
                                                                                                                                                    |                                                                                                                                                    
gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAg                                                                                                                                                     >  2:64788/1‑149 (MQ=255)
 aTTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGc                                                                                                                                                    >  3:308976/1‑149 (MQ=255)
   tGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTg                                                                                                                                                  <  2:469554/149‑1 (MQ=255)
    gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgg                                                                                                                                                 >  1:694044/1‑149 (MQ=255)
     cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTggg                                                                                                                                                <  3:289810/149‑1 (MQ=255)
      gcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCTGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgggg                                                                                                                                               <  3:642258/149‑1 (MQ=255)
       cATGGCAGGATGAGACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGga                                                                                                                                              <  3:63004/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                              >  3:754808/1‑149 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                              <  1:37301/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGCGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                              <  2:389197/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACAAGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                              <  4:542804/149‑1 (MQ=255)
        aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaa                                                                                                                                             <  1:601948/149‑2 (MQ=255)
         tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaat                                                                                                                                            <  1:43218/149‑3 (MQ=255)
          ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAACCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGAt                                                                                                                                           <  4:464241/149‑1 (MQ=255)
            cAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCt                                                                                                                                         >  3:328425/1‑149 (MQ=255)
             aGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCtt                                                                                                                                        >  3:469072/1‑149 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                                       <  1:670530/149‑1 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                                       <  3:88248/149‑1 (MQ=255)
              ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTGGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                                       <  2:750530/149‑1 (MQ=255)
                 tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                                    <  3:10788/149‑1 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   <  4:389160/149‑1 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   >  2:332601/1‑149 (MQ=255)
                  gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCAGGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   >  3:525842/1‑149 (MQ=255)
                    aaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTTCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                                 >  4:711408/1‑149 (MQ=255)
                    aaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                                 >  2:661390/1‑149 (MQ=255)
                     aCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  4:629321/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  3:619419/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  4:390684/1‑149 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  4:71541/1‑149 (MQ=255)
                            tACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGATATCTTTGGTCAGCACGTTCa                                                                                                                         >  4:438169/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCGGCACGTTCACTgcg                                                                                                                    >  3:244759/1‑149 (MQ=255)
                                      gATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                               <  4:308976/149‑1 (MQ=255)
                                        ttATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  2:401642/1‑149 (MQ=255)
                                              gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                       >  4:92979/1‑149 (MQ=255)
                                              gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                       <  4:570730/149‑1 (MQ=255)
                                                      gtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCa                                                                                               <  1:332601/149‑1 (MQ=255)
                                                          ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  2:288431/1‑149 (MQ=255)
                                                          ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  3:140347/1‑149 (MQ=255)
                                                             aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                        >  3:333843/1‑149 (MQ=255)
                                                             aaaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                        >  3:201337/1‑149 (MQ=255)
                                                              aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  2:698271/1‑149 (MQ=255)
                                                               aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                                      <  1:64788/149‑1 (MQ=255)
                                                                     tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                                >  1:326820/1‑149 (MQ=255)
                                                                     gaaGTGTCAGTTCCGGAGCGTTGTATTTTTCCAGCCAGGCGAGCGCCGTATTCCCCGCACATCCATTCCTCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAGACCGTCAGCGGTTAAATCCAGa                                                                                <  4:54378/148‑1 (MQ=255)
                                                                          gtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCgg                                                                           <  3:153823/149‑1 (MQ=255)
                                                                           tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCACCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                                                                          >  1:710417/1‑149 (MQ=255)
                                                                              gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                                       <  2:566824/149‑1 (MQ=255)
                                                                                 ccGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                        <  1:270854/145‑1 (MQ=255)
                                                                                 ccGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                        >  2:270854/1‑145 (MQ=255)
                                                                                  cGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGccc                                                                   <  3:426103/149‑1 (MQ=255)
                                                                                    gACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             <  4:529128/105‑1 (MQ=255)
                                                                                    gACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  3:529128/1‑105 (MQ=255)
                                                                                       cGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCCCAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  <  2:675613/129‑1 (MQ=255)
                                                                                       cGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  >  1:675613/1‑129 (MQ=255)
                                                                                        gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                             <  2:425120/149‑1 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  2:262090/1‑149 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  1:547262/1‑149 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  3:246961/1‑149 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  3:408313/1‑149 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  2:744701/1‑149 (MQ=255)
                                                                                            tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt                                                         <  3:71541/149‑1 (MQ=255)
                                                                                             aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                        >  2:627240/1‑149 (MQ=255)
                                                                                              tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGc                                                       >  1:749677/1‑149 (MQ=255)
                                                                                                 ttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaaa                                                    <  4:754808/149‑1 (MQ=255)
                                                                                                    cAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCa                                                 <  2:115138/149‑1 (MQ=255)
                                                                                                          tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAg                                           <  4:408313/149‑1 (MQ=255)
                                                                                                          tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAg                                           <  2:352341/149‑1 (MQ=255)
                                                                                                           gCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAgg                                          >  3:2117/1‑149 (MQ=255)
                                                                                                             cAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                                >  2:230186/1‑45 (MQ=25)
                                                                                                             cAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                                <  1:230186/45‑1 (MQ=25)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                                >  2:681032/1‑149 (MQ=255)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTACACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                                >  2:632807/1‑149 (MQ=255)
                                                                                                                      aTTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCt                               <  2:547262/149‑1 (MQ=255)
                                                                                                                       ttACCCGCACAGCCATTCCCCAGCCGCAAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATACCCTc                              >  2:755479/1‑149 (MQ=255)
                                                                                                                             gCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                        >  1:204615/1‑149 (MQ=255)
                                                                                                                                caGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAc                     <  3:642858/149‑1 (MQ=255)
                                                                                                                                  gCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCg                   >  2:206390/1‑149 (MQ=255)
                                                                                                                                   ccATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc                  >  2:282140/1‑149 (MQ=255)
                                                                                                                                    cATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCt                 >  3:679212/1‑149 (MQ=255)
                                                                                                                                       tCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg              >  2:653501/1‑149 (MQ=255)
                                                                                                                                          ccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa           >  2:39870/1‑149 (MQ=255)
                                                                                                                                           cAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCaa          <  1:698271/149‑1 (MQ=255)
                                                                                                                                                    aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  2:694044/149‑1 (MQ=255)
                                                                                                                                                    aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATACCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  3:645331/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  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: