Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,163:1 +T 7.4% 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.T 7.4% 714.9 / 55.0 176coding (2232/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base . (78/85);  new base T (6/7);  total (84/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 177 total aligned reads displayed.

TCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAA  >  NC_000913/3714019‑3714308
                                                                                                                                                  |                                                                                                                                                 
 cGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGa                                                                                                                                              <  4:1560926/149‑1 (MQ=255)
    ttGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCt                                                                                                                                           <  3:467499/149‑1 (MQ=255)
       cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                        <  3:1393742/149‑1 (MQ=255)
       cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTggg                                                                                                                                        <  1:475182/149‑1 (MQ=255)
        gcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgggg                                                                                                                                       <  1:123454/149‑1 (MQ=255)
         cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGga                                                                                                                                      <  1:1005209/149‑1 (MQ=255)
          aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGGAAGTTTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGgag                                                                                                                                     <  4:1500047/149‑1 (MQ=255)
           tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaat                                                                                                                                    <  1:1495527/149‑3 (MQ=255)
                ggATGAAACCCAACTACCCACATGGATTATGTGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                               <  1:1637069/149‑1 (MQ=255)
                ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                               <  3:667442/149‑1 (MQ=255)
                 gATGAAACCCAACTACCCACATGGCTTATGAGCTGGCCGCTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGAGGAGATCTTTg                                                                                                                              <  3:1110666/149‑1 (MQ=255)
                 gATGAAACCCAACTACCCACATGGATTATGTGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTg                                                                                                                              <  1:257598/149‑1 (MQ=255)
                 gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTg                                                                                                                              <  1:182378/149‑1 (MQ=255)
                   tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                            <  1:565266/149‑1 (MQ=255)
                   tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                            >  4:1650814/1‑149 (MQ=255)
                   tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                            >  1:1637934/1‑149 (MQ=255)
                    gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                           >  3:225119/1‑149 (MQ=255)
                     aaaCCCAACTGCCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGCCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCTCAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCa                                                                                                                          >  4:1096258/1‑149 (MQ=255)
                       aCCCAACGACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                        >  2:932686/1‑149 (MQ=255)
                          cAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                     >  2:1356100/1‑149 (MQ=255)
                                cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                        <  1:1119450/140‑1 (MQ=255)
                                cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACt                                                                                                               <  1:67289/149‑1 (MQ=255)
                                cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGCAGATCTTTGGTCAGc                                                                                                                        >  2:1119450/1‑140 (MQ=255)
                                 ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCtt                                                                                                                                >  4:1203868/1‑131 (MQ=255)
                                 ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                                                              >  2:468293/1‑149 (MQ=255)
                                  cacaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgc                                                                                                             <  2:866009/149‑1 (MQ=255)
                                   acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                            >  1:835827/1‑149 (MQ=255)
                                    caTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                           >  4:147899/1‑149 (MQ=255)
                                     aTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                                                                          <  1:913902/149‑1 (MQ=255)
                                     aTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                                                                          <  1:1208291/149‑1 (MQ=255)
                                     aTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                                                                          >  3:943954/1‑149 (MQ=255)
                                                  tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATa                                                                                             >  4:811637/1‑149 (MQ=255)
                                                  tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCAGCCCAGGCCAGCGCGGTATTACCGCACCAGCAATCCCCCAGCCGCGAGAGGGGGATATCTTTGGTCACCACGTTCACGGCCCCGTTTTCACaattc                                                                                             >  2:1545536/1‑146 (MQ=255)
                                                   ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAc                                                                                            >  2:261520/1‑149 (MQ=255)
                                                   ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAc                                                                                            >  3:1672325/1‑149 (MQ=255)
                                                    gCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAcc                                                                                           >  4:947752/1‑149 (MQ=255)
                                                       ggtggtTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTc                                                                                        >  4:1011352/1‑149 (MQ=255)
                                                             tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                        <  3:969933/127‑1 (MQ=255)
                                                                    cTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATcc                                                                           <  3:811637/149‑1 (MQ=255)
                                                                     tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                          >  3:486806/1‑149 (MQ=255)
                                                                     tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCATATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                          <  1:1332965/149‑1 (MQ=255)
                                                                      gTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                         <  3:613580/149‑1 (MQ=255)
                                                                       tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                            >  2:1384997/1‑97 (MQ=255)
                                                                       tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGt                                                                                                                            <  1:1384997/97‑1 (MQ=255)
                                                                       tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                        <  2:1330383/149‑1 (MQ=255)
                                                                        aaGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                       >  1:1039261/1‑149 (MQ=255)
                                                                        aaGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                       >  1:1359765/1‑149 (MQ=255)
                                                                          gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                       <  2:643373/147‑1 (MQ=255)
                                                                          gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                       >  1:643373/1‑147 (MQ=255)
                                                                          gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATcc                                                                     >  3:62082/1‑149 (MQ=255)
                                                                           tgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                                                                    >  4:1505235/1‑149 (MQ=255)
                                                                            gtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCgg                                                                   >  1:335312/1‑149 (MQ=255)
                                                                               aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                              <  4:712242/135‑1 (MQ=255)
                                                                               aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                              >  3:712242/1‑135 (MQ=255)
                                                                               aGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                <  1:468293/149‑1 (MQ=255)
                                                                                gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                               >  3:530994/1‑149 (MQ=255)
                                                                                 ttCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                              >  4:447512/1‑149 (MQ=255)
                                                                                 ttCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                              >  2:1207973/1‑149 (MQ=255)
                                                                                 ttCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                              >  3:163200/1‑149 (MQ=255)
                                                                                  tCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                             <  3:121234/149‑1 (MQ=255)
                                                                                  tCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                             <  1:1307763/149‑1 (MQ=255)
                                                                                   ccGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcc                                                            >  3:105975/1‑149 (MQ=255)
                                                                                        ccGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc                                                       <  2:335312/149‑1 (MQ=255)
                                                                                          gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                     <  4:1308938/149‑1 (MQ=255)
                                                                                           ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                    >  3:588999/1‑149 (MQ=255)
                                                                                           ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                    >  1:55922/1‑149 (MQ=255)
                                                                                             gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                  <  4:976516/149‑1 (MQ=255)
                                                                                              tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                  <  1:1078088/148‑1 (MQ=255)
                                                                                              tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt                                                 <  2:55922/149‑1 (MQ=255)
                                                                                               aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                                >  2:1172934/1‑149 (MQ=255)
                                                                                                tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                     <  2:679208/143‑1 (MQ=255)
                                                                                                   ttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaaa                                            <  3:1011352/149‑1 (MQ=255)
                                                                                                       aGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCAc                                        <  3:1526010/149‑1 (MQ=255)
                                                                                                       aGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTGTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCAc                                        <  4:850057/149‑1 (MQ=255)
                                                                                                         ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGc                                      <  3:924784/149‑1 (MQ=255)
                                                                                                            tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                  >  4:1324642/1‑134 (MQ=255)
                                                                                                            tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                  <  3:1324642/134‑1 (MQ=255)
                                                                                                              ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                 <  3:877628/149‑1 (MQ=255)
                                                                                                                aGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACTGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGttt                               >  3:1625048/1‑149 (MQ=255)
                                                                                                                aGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGttt                               >  4:1398984/1‑149 (MQ=255)
                                                                                                                   gcCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAt                            <  1:796942/149‑1 (MQ=255)
                                                                                                                    ccGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                 <  3:446492/143‑1 (MQ=255)
                                                                                                                    ccGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                 >  4:446492/1‑143 (MQ=255)
                                                                                                                       tATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                        >  4:945956/1‑149 (MQ=255)
                                                                                                                       tATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                        <  2:290409/149‑1 (MQ=255)
                                                                                                                              cGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTaa                 <  2:1435333/149‑1 (MQ=255)
                                                                                                                               gCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                >  2:779770/1‑149 (MQ=255)
                                                                                                                               gCACAGCCATTCCCCAGCAGTGAGCTGGGGAGATCTTTGGTCTGTACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                <  2:1328710/149‑1 (MQ=255)
                                                                                                                                      ccATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCa                                                                <  2:1046241/95‑1 (MQ=255)
                                                                                                                                       cATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCt         >  4:488888/1‑149 (MQ=255)
                                                                                                                                       atctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                   >  2:296077/4‑75 (MQ=37)
                                                                                                                                       atctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                   <  1:296077/72‑1 (MQ=37)
                                                                                                                                          tCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg      >  2:911308/1‑149 (MQ=255)
                                                                                                                                           ccccAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCgg     >  3:1513784/1‑149 (MQ=255)
                                                                                                                                            cccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc    <  4:225119/149‑1 (MQ=255)
                                                                                                                                            cccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc    <  3:1142977/149‑1 (MQ=255)
                                                                                                                                            cccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc    <  1:1356100/149‑1 (MQ=255)
                                                                                                                                             ccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa   >  4:1021957/1‑149 (MQ=255)
                                                                                                                                             ccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa   >  1:1619721/1‑149 (MQ=255)
                                                                                                                                              gaGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCCCTGCGCCGTTCTTACAAGTACCGTGAGAGGTTAAATGCAGATCCGGCCGTGCCGCTTCGGGAATGCCAGTCACGCCAGGTTTAATTCCCTCGCTCATGACCGCTCCGGCaa  >  2:979064/2‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAA  >  NC_000913/3714019‑3714308

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