Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,168 Δ1 bp 6.7% 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. 6.7% 303.9 / 14.7 90coding (2227/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base A (45/39);  new base . (3/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-02

ATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714022‑3714316
                                                                                                                                                   |                                                                                                                                                    
aTTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAg                                                                                                                                                     >  3:403738/1‑148 (MQ=255)
   gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTg                                                                                                                                                  <  4:403738/148‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTggg                                                                                                                                                <  2:11905/149‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgg                                                                                                                                                 <  3:59904/148‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgg                                                                                                                                                 <  3:826565/148‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgg                                                                                                                                                 <  4:1017574/148‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCCAGCTggg                                                                                                                                                <  2:579010/149‑1 (MQ=255)
    cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTAAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgg                                                                                                                                                 <  3:139366/148‑1 (MQ=255)
     gcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTggg                                                                                                                                                <  3:395866/148‑1 (MQ=255)
      cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTgggg                                                                                                                                               <  3:126148/148‑1 (MQ=255)
        tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaat                                                                                                                                            <  1:198523/149‑3 (MQ=255)
        tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCGGCTGGGgaa                                                                                                                                             <  4:928096/148‑2 (MQ=255)
           cAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATc                                                                                                                                          >  4:1057712/1‑148 (MQ=255)
             ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCtt                                                                                                                                        <  4:1034990/148‑1 (MQ=255)
              gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCttt                                                                                                                                       <  3:1057712/148‑1 (MQ=255)
                  aaaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTc                                                                                                                                   <  4:624635/148‑1 (MQ=255)
                    aCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGCCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  2:530584/1‑149 (MQ=255)
                     cccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  4:336847/1‑148 (MQ=255)
                       cAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                             >  2:589224/1‑149 (MQ=255)
                       cAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                             >  1:895726/1‑149 (MQ=255)
                          cTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGATCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGtt                                                                                                                           >  4:730291/1‑148 (MQ=255)
                            aCCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGGCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCTCGAACTGGGGAGATTTTTGGTCAGccgtttca                                                                                                                         <  4:486402/148‑8 (MQ=255)
                             cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCAc                                                                                                                        <  4:575168/148‑1 (MQ=255)
                               cacaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgc                                                                                                                     >  3:1014014/1‑148 (MQ=255)
                                acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                <  2:998564/137‑1 (MQ=255)
                                acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  1:998564/1‑137 (MQ=255)
                                     gATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                                <  3:475789/148‑1 (MQ=255)
                                      aTTATGAGCTGGCCGGTGGTTCATACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCGCCAGCCGCGAGCTGGGGAGAGCTTTGGTCAGCACGTTCAGTGCGCCGtt                                                                                                               >  3:892761/1‑148 (MQ=255)
                                      aTTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTATATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATGCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGTGCCGtt                                                                                                               <  4:896632/148‑1 (MQ=255)
                                        tATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  4:141450/1‑148 (MQ=255)
                                           gaccTGGCCGGTCGTTCAAACGCTGTAAGTGTCATTTCCGGACCGTTGTATTGTTCCAGCCAGGCCAGCGCATTATTACCCGCACCGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGCCTTTAc                                                                                                          >  4:486852/4‑148 (MQ=255)
                                             gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                                   <  2:566286/137‑1 (MQ=255)
                                             gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                                   >  1:566286/1‑137 (MQ=255)
                                             gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaa                                                                                                        <  4:695878/148‑1 (MQ=255)
                                                         ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  1:218612/1‑149 (MQ=255)
                                                             aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  1:781127/1‑149 (MQ=255)
                                                             aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                        >  4:126468/1‑148 (MQ=255)
                                                              aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  4:16692/1‑148 (MQ=255)
                                                               cGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                                      >  3:920685/1‑148 (MQ=255)
                                                               cGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                                      >  4:878390/1‑148 (MQ=255)
                                                                gCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    <  1:893839/149‑1 (MQ=255)
                                                                  tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATcc                                                                                   <  3:141450/148‑1 (MQ=255)
                                                                  tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  <  2:898558/149‑1 (MQ=255)
                                                                          tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCAATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCTGATCCGGc                                                                          >  2:879505/1‑149 (MQ=255)
                                                                               ttcGGACCGTTGTATTTTTCTAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                                     <  2:1001322/147‑1 (MQ=255)
                                                                                  ggACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcccc                                                                  >  2:532697/1‑149 (MQ=255)
                                                                                    aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    >  4:556230/1‑129 (MQ=255)
                                                                                    aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    <  3:556230/129‑1 (MQ=255)
                                                                                     ccGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt                                                                <  4:63147/148‑1 (MQ=255)
                                                                                        ttGTATTTTTCCAGCCATGCCAGCGCCGCATTACCCGCATAGTGATTTACCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGGTGTTACAGATACATTCTGCGTTTAAACCCTGATCCGGCCATGCCTCTTAGTg                                                            >  2:287087/1‑149 (MQ=255)
                                                                                         tGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                            >  3:266816/1‑148 (MQ=255)
                                                                                         tGTATTTATCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa                                                           <  1:935544/149‑1 (MQ=255)
                                                                                           tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  3:1029967/1‑148 (MQ=255)
                                                                                                ttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaaa                                                    <  1:836051/149‑1 (MQ=255)
                                                                                                      ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACg                                               <  3:126468/148‑1 (MQ=255)
                                                                                                               cgcCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTaa                                     <  1:855712/149‑1 (MQ=255)
                                                                                                                gcCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAt                                    <  2:639176/149‑1 (MQ=255)
                                                                                                                 ccGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAt                                    >  4:698924/1‑148 (MQ=255)
                                                                                                                     aTTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGTGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCAGGCCATCCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                                >  4:87491/1‑148 (MQ=255)
                                                                                                                       tACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTc                              >  4:361265/1‑148 (MQ=255)
                                                                                                                          ccGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCGCGCTa                          >  2:618976/1‑149 (MQ=255)
                                                                                                                           cGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTa                          >  3:616655/1‑148 (MQ=255)
                                                                                                                           cGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTa                          <  3:730291/148‑1 (MQ=255)
                                                                                                                           cGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTa                          >  3:426184/1‑148 (MQ=255)
                                                                                                                            gCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTaa                         <  4:516418/148‑1 (MQ=255)
                                                                                                                              acaGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGa                      >  2:904101/1‑149 (MQ=255)
                                                                                                                               caGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          <  3:978260/112‑1 (MQ=255)
                                                                                                                               caGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  4:978260/1‑112 (MQ=255)
                                                                                                                               caGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAc                     >  2:359569/1‑149 (MQ=255)
                                                                                                                                 gCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCg                   >  1:585038/1‑149 (MQ=255)
                                                                                                                                  ccATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc                  >  2:245164/1‑149 (MQ=255)
                                                                                                                                    aTTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTc                >  2:215235/1‑149 (MQ=255)
                                                                                                                                    aTTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTc                >  2:532281/1‑149 (MQ=255)
                                                                                                                                     ttCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTcc               >  1:164487/1‑149 (MQ=255)
                                                                                                                                     ctCCCCAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTc                >  3:207913/2‑148 (MQ=255)
                                                                                                                                      tCCCCAGCCGCGAGCTGGGGAGTTCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                                     >  1:652681/1‑94 (MQ=255)
                                                                                                                                      tCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                                     <  2:652681/94‑1 (MQ=255)
                                                                                                                                      cccccAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCt                               <  4:15676/131‑1 (MQ=255)
                                                                                                                                      cccccAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCt                               >  3:15676/2‑132 (MQ=255)
                                                                                                                                       ccccAGCTCGCGGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTGTAATTCCCTCGCTAATGACCGCTCCg              >  4:437754/1‑148 (MQ=255)
                                                                                                                                       ccccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg              <  4:755304/148‑1 (MQ=255)
                                                                                                                                         ccAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGc            >  4:826792/1‑148 (MQ=255)
                                                                                                                                          cAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCaa          >  2:632615/1‑149 (MQ=255)
                                                                                                                                             ccGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCGCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGaa       <  2:806551/149‑1 (MQ=255)
                                                                                                                                             ccGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGa        <  4:90453/148‑1 (MQ=255)
                                                                                                                                                gcgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAtt     >  2:51543/1‑149 (MQ=255)
                                                                                                                                                gcgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAtt     >  2:117523/1‑149 (MQ=255)
                                                                                                                                                 cgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCaaattct     <  4:266816/148‑6 (MQ=255)
                                                                                                                                                 cgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAtt     <  3:361265/148‑1 (MQ=255)
                                                                                                                                                  gAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAAttt    >  3:496274/1‑148 (MQ=255)
                                                                                                                                                   aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  1:51543/149‑1 (MQ=255)
                                                                                                                                                   aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  1:117523/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714022‑3714316

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