Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,163:1 +T 5.5% 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 5.5% 604.6 / 31.0 145coding (2232/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base . (70/67);  new base T (2/6);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 148 total aligned reads displayed.

TTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCA  >  NC_000913/3714016‑3714307
                                                                                                                                                    |                                                                                                                                                
ttGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc‑c                                                                                                                                                 >  3:299208/1‑149 (MQ=255)
         gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                        <  3:716992/149‑1 (MQ=255)
          cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCCGCTCGCG‑GCTggg                                                                                                                                       <  3:1291943/149‑1 (MQ=255)
          cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTggg                                                                                                                                       <  2:652600/149‑1 (MQ=255)
            cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGga                                                                                                                                     <  2:1020367/149‑1 (MQ=255)
             aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaa                                                                                                                                    <  3:1065699/149‑2 (MQ=255)
             aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgag                                                                                                                                    >  1:1489625/1‑149 (MQ=255)
               ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaatg                                                                                                                                  <  1:206855/149‑4 (MQ=255)
               ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGAt                                                                                                                                  <  3:963645/149‑1 (MQ=255)
                gCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATc                                                                                                                                 >  3:706239/1‑149 (MQ=255)
                 cAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCt                                                                                                                                >  3:1206945/1‑149 (MQ=255)
                  aGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCtt                                                                                                                               >  3:1479875/1‑149 (MQ=255)
                    gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                             <  4:1443767/149‑1 (MQ=255)
                    gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                             <  3:888105/149‑1 (MQ=255)
                    gATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                             <  3:747167/149‑1 (MQ=255)
                    gATAAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                             <  2:234368/149‑1 (MQ=255)
                        aaaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTACGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCa                                                                                                                         <  4:1233295/149‑1 (MQ=255)
                          aCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                       >  3:128005/1‑149 (MQ=255)
                            ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                     >  2:1271171/1‑149 (MQ=255)
                            ccAACTACCCACATGGATAATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                     >  4:1311842/1‑149 (MQ=255)
                              aaCTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACCGCCATTCCCCAGC‑CGCGAGCTGGGGAGACCTTTGGTCAGCACga                                                                                                                   >  4:69208/1‑148 (MQ=255)
                                    ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                                                             >  1:1245058/1‑149 (MQ=255)
                                      acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGctcgc                                                                                                                                               <  3:803843/114‑5 (MQ=255)
                                      acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                           >  1:563234/1‑149 (MQ=255)
                                      acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                           >  3:385631/1‑149 (MQ=255)
                                        aTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                                                                         >  3:1442810/1‑149 (MQ=255)
                                          ggATTATGAGCTGGCCGGTGGTTCAAACGCTGTACGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGTACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                       >  2:14490/1‑149 (MQ=255)
                                          ggATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                       >  1:1108896/1‑149 (MQ=255)
                                           gATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                      <  2:1349187/149‑1 (MQ=255)
                                            aTTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                                                                     >  1:996007/1‑149 (MQ=255)
                                             ttATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                    >  2:644325/1‑149 (MQ=255)
                                              tATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttttt                                                                                                   <  3:681695/149‑1 (MQ=255)
                                               aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTa                                                                                                  <  2:535936/149‑1 (MQ=255)
                                               aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTa                                                                                                  <  1:381642/149‑1 (MQ=255)
                                               aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTa                                                                                                  <  4:1479875/149‑1 (MQ=255)
                                                   gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                              <  4:256160/149‑1 (MQ=255)
                                                      ggCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAc                                                                                           <  3:721576/149‑1 (MQ=255)
                                                               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                  <  4:475907/149‑1 (MQ=255)
                                                               ttCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                  <  1:148913/149‑1 (MQ=255)
                                                                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGtt                                                                                                                  <  3:19937/116‑1 (MQ=255)
                                                                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGtt                                                                                                                  >  4:19937/1‑116 (MQ=255)
                                                                tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGt                                                                                 >  3:656823/1‑149 (MQ=255)
                                                                 cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                                                                <  3:278648/149‑1 (MQ=255)
                                                                 cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGtt                                                                                <  3:1442496/149‑1 (MQ=255)
                                                                   aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                              >  1:1362151/1‑101 (MQ=255)
                                                                   aaCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                              <  2:1362151/101‑1 (MQ=255)
                                                                    aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaaa                                                                             <  1:1391136/149‑1 (MQ=255)
                                                                     cGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAAt                                                                            <  3:996546/149‑1 (MQ=255)
                                                                      gCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                           <  2:1489625/149‑1 (MQ=255)
                                                                        tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGCGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                         <  4:1013942/149‑1 (MQ=255)
                                                                         gTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                        >  2:229658/1‑149 (MQ=255)
                                                                          tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                                                                     >  1:132004/1‑119 (MQ=255)
                                                                           aaGTGTCAGTTCCGGACCGTTGTATTTTTGCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTACATCCAGAt                                                                      >  3:882995/1‑149 (MQ=255)
                                                                             gtgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATcc                                                                    >  4:299943/1‑149 (MQ=255)
                                                                              tgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                                                                   >  4:771085/1‑149 (MQ=255)
                                                                                tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                              >  2:47120/1‑88 (MQ=255)
                                                                                tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                              <  1:47120/88‑1 (MQ=255)
                                                                                tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                                                                 <  4:706239/149‑1 (MQ=255)
                                                                                   gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                              >  1:1431936/1‑149 (MQ=255)
                                                                                   gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                              >  1:894349/1‑149 (MQ=255)
                                                                                    ttCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                             >  1:1043311/1‑149 (MQ=255)
                                                                                     tCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGc                                                            <  4:745093/149‑1 (MQ=255)
                                                                                             gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                    <  1:198383/149‑1 (MQ=255)
                                                                                               tGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGTCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGa                                                  <  4:888290/149‑1 (MQ=255)
                                                                                               tGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttttt                                                                                                   <  3:1449296/100‑1 (MQ=255)
                                                                                               tGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttttt                                                                                                   >  4:1449296/1‑100 (MQ=255)
                                                                                                gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                 <  4:258493/149‑1 (MQ=255)
                                                                                                 tATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt                                                <  3:618051/149‑1 (MQ=255)
                                                                                                  aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATg                                               <  1:1271171/149‑1 (MQ=255)
                                                                                                   tttttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGc                                              <  1:804541/149‑1 (MQ=255)
                                                                                                     tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaa                                            >  2:973693/1‑149 (MQ=255)
                                                                                                       tCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAAt                                          >  1:24475/1‑149 (MQ=255)
                                                                                                          aGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCAc                                       >  3:1224682/1‑149 (MQ=255)
                                                                                                           gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACg                                      >  4:13270/1‑149 (MQ=255)
                                                                                                              aTGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                   <  3:532607/149‑1 (MQ=255)
                                                                                                               tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAg                                  <  1:283416/149‑1 (MQ=255)
                                                                                                                 ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                >  2:610214/1‑149 (MQ=255)
                                                                                                                 ccAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGt                                >  4:659872/1‑149 (MQ=255)
                                                                                                                  cAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGtt                               <  2:111357/149‑1 (MQ=255)
                                                                                                                          tATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                       <  4:506391/149‑1 (MQ=255)
                                                                                                                           aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGTCAGGTTTAAtt                          <  3:958172/145‑1 (MQ=255)
                                                                                                                           aTTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCt                      <  1:504751/149‑1 (MQ=255)
                                                                                                                            ttACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTc                     <  3:299943/149‑1 (MQ=255)
                                                                                                                            ttACCCGCACAGCCATTACCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGATTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCGTTCGTGAATGCAAATCACGCCAGTTTTAATTCCCTc                     <  4:1253343/149‑1 (MQ=255)
                                                                                                                              aCCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGc                   <  3:659872/149‑1 (MQ=255)
                                                                                                                               cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCt                  <  2:483399/149‑1 (MQ=255)
                                                                                                                                ccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAGCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTa                 >  1:1430790/1‑149 (MQ=255)
                                                                                                                                  gCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                        >  2:115658/1‑92 (MQ=255)
                                                                                                                                  gCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                        <  1:115658/92‑1 (MQ=255)
                                                                                                                                  gCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCATCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt               <  3:918796/149‑1 (MQ=255)
                                                                                                                                    acaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGt                                                                                                       <  4:1150440/59‑1 (MQ=255)
                                                                                                                                     caGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAc            <  4:644493/149‑1 (MQ=255)
                                                                                                                                        ccATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGACAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc         >  2:551761/1‑149 (MQ=255)
                                                                                                                                          aTTCCCCAGC‑CGCGAGCTTGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                       <  3:889445/133‑1 (MQ=255)
                                                                                                                                          aTTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                       >  4:889445/1‑133 (MQ=255)
                                                                                                                                           ttCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTcc      >  2:1175696/1‑149 (MQ=255)
                                                                                                                                            tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg     >  4:356082/1‑149 (MQ=255)
                                                                                                                                               ccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt               <  1:1314326/136‑1 (MQ=255)
                                                                                                                                               ccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt               >  2:1314326/1‑136 (MQ=255)
                                                                                                                                               ccAGC‑CGCGAGCTGGGGAGAACTTTGGTCAGCACGTTCACTGCGCCGTTTTTACATATACCGTCAGCGGTTTAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa  >  2:609499/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCA  >  NC_000913/3714016‑3714307

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: ‑