Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,500 T→C 5.0% intergenic (‑113/+54) otsB ← / ← araH trehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,5000TC 5.0% 391.3 / 6.9 119intergenic (‑113/+54)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (61/52);  new base C (3/3);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

Only 100 of 123 total aligned reads displayed.

GGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAA  >  NC_000913/1982352‑1982648
                                                                                                                                                    |                                                                                                                                                     
ggATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGt                                                                                                                                                       >  1:388887/1‑149 (MQ=255)
 gATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                      <  2:814150/149‑1 (MQ=255)
   tAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTAt                                                                                                                                                    <  4:69351/149‑1 (MQ=255)
   tAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTAt                                                                                                                                                    >  3:463274/1‑149 (MQ=255)
     gTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATgg                                                                                                                                                  <  4:133406/149‑1 (MQ=255)
         aGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGCCCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGt                                                                                                                                              <  2:309180/149‑1 (MQ=255)
          ggggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                             <  3:215609/149‑1 (MQ=255)
          ggggTTTCGGTAAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTGTACTGACATGTAAGTATAGACAGGGTGACGGGGGTAGACATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                             <  1:733748/149‑1 (MQ=255)
           gggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTg                                                                                                                                            >  3:781137/1‑149 (MQ=255)
            ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTgg                                                                                                                                           <  4:111265/149‑1 (MQ=255)
             gTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTAGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                          <  2:74948/149‑1 (MQ=255)
              tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTAt                                                                                                                                                    >  2:257760/1‑138 (MQ=255)
              tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                         >  4:441007/1‑149 (MQ=255)
              tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                         >  4:314465/1‑149 (MQ=255)
                   gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                    >  2:620159/1‑149 (MQ=255)
                          gTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaa                                                                                                                            >  4:852639/1‑149 (MQ=255)
                           ttCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaa                                                                                                                           >  2:565494/1‑149 (MQ=255)
                            tCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aacaaa                                                                                                                          <  1:771876/149‑1 (MQ=255)
                            tCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaaa                                                                                                                          <  4:781137/149‑1 (MQ=255)
                                tCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCt                                                                                                                      <  3:441007/149‑1 (MQ=255)
                                tCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCt                                                                                                                      >  2:511212/1‑149 (MQ=255)
                                 cACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTa                                                                                                                     >  1:452955/1‑149 (MQ=255)
                                   cccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTaaa                                                                                                                   >  3:624534/1‑149 (MQ=255)
                                   cccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTaaa                                                                                                                   >  1:287767/1‑149 (MQ=255)
                                           tcctccTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCa                                                                                                           >  1:639299/1‑149 (MQ=255)
                                             ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGa                                                                                                         >  1:275132/1‑149 (MQ=255)
                                             ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCAAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGGTAAATTGTTGCAGa                                                                                                         >  1:608948/1‑149 (MQ=255)
                                              tcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaa                                                                                                        >  4:226551/1‑149 (MQ=255)
                                              tcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaa                                                                                                        >  1:307274/1‑149 (MQ=255)
                                                ctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATATACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAAc                                                                                                                         >  1:356310/1‑130 (MQ=255)
                                                ctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAAc                                                                                                                         <  2:356310/130‑1 (MQ=255)
                                                ctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaaaa                                                                                                      >  1:111754/1‑149 (MQ=255)
                                                 ttCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaaaaa                                                                                                     >  2:157222/1‑149 (MQ=255)
                                                 ttCTCTTTTCTTATTGTGTTCCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCTCCATTTGCAGCGAAGAACAGCGAACCGTAATGGCGGGTGACTATGGTGGG‑TAAAGACGCTAAGTTGGTGTAGaaaaa                                                                                                     >  2:201551/1‑149 (MQ=255)
                                                  tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaaaaaa                                                                                                    <  4:463274/149‑1 (MQ=255)
                                                    tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGc                                                                                                  <  1:430760/149‑1 (MQ=255)
                                                      ttttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAGCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCAt                                                                                                <  3:62446/149‑1 (MQ=255)
                                                      ttttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCAt                                                                                                >  4:503681/1‑149 (MQ=255)
                                                      ttttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCAt                                                                                                >  1:894933/1‑149 (MQ=255)
                                                         tCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAg                                                                                             <  1:157222/149‑1 (MQ=255)
                                                           ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTgg                                                                                                                                           >  4:452924/1‑102 (MQ=255)
                                                            tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaaa                                                                                                                          >  2:526556/1‑117 (MQ=255)
                                                            tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaaa                                                                                                                          <  1:526556/117‑1 (MQ=255)
                                                            tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGaaa                                                                                          <  2:659798/149‑3 (MQ=255)
                                                              ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGt                                                                                        <  2:111754/149‑1 (MQ=255)
                                                              ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGt                                                                                        >  1:781390/1‑149 (MQ=255)
                                                              ttGTGTTGCGCAGAACGAGCTTTACAGACATGTAAGTATAGACAGTGTGTCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACAGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGGATCAGACAGt                                                                                        >  4:252005/1‑149 (MQ=255)
                                                                 tgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCg                                                                                     <  1:855887/149‑1 (MQ=255)
                                                                   ttacgcAGAACGCGCTTTACAGACTTGTAAGTATAGACAGTGTGACGGGGATCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AACAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACCGTGCGtt                                                                                   >  2:351890/4‑149 (MQ=255)
                                                                    tGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCACTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGttt                                                                                  <  3:708982/149‑1 (MQ=255)
                                                                     gcgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc                                                                                 <  2:94951/149‑1 (MQ=255)
                                                                            aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                          <  1:292209/149‑1 (MQ=255)
                                                                            aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                          >  3:315193/1‑149 (MQ=255)
                                                                            aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                          >  1:108487/1‑149 (MQ=255)
                                                                              gcgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                        <  2:728403/149‑1 (MQ=255)
                                                                              gcgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                        <  3:551720/149‑1 (MQ=255)
                                                                               cgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCt                                                                       <  3:150212/149‑1 (MQ=255)
                                                                                gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                      <  3:506245/149‑1 (MQ=255)
                                                                                gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                      <  2:387830/149‑1 (MQ=255)
                                                                                       aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                               >  2:533447/1‑149 (MQ=255)
                                                                                        gACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGt                                                              <  2:856312/149‑1 (MQ=255)
                                                                                           aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                           >  1:543944/1‑149 (MQ=255)
                                                                                           aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                           >  4:819648/1‑149 (MQ=255)
                                                                                           aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                           >  1:863252/1‑149 (MQ=255)
                                                                                               aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGAt                                                       <  3:286434/149‑1 (MQ=255)
                                                                                                  tataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                    >  1:643731/1‑149 (MQ=255)
                                                                                                     aGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGc                                                 >  2:95507/1‑149 (MQ=255)
                                                                                                              gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagcag                                        <  3:503681/149‑1 (MQ=255)
                                                                                                              gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagcag                                        <  1:533447/149‑1 (MQ=255)
                                                                                                                 gggggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAt                                     >  1:162363/1‑149 (MQ=255)
                                                                                                                     gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAg                                 <  2:307274/149‑1 (MQ=255)
                                                                                                                      tCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                                >  2:775965/1‑149 (MQ=255)
                                                                                                                      tCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                                >  1:784077/1‑149 (MQ=255)
                                                                                                                            tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGaa                          >  2:781479/1‑149 (MQ=255)
                                                                                                                            tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGaa                          <  4:461414/149‑1 (MQ=255)
                                                                                                                                cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                      <  2:305917/149‑1 (MQ=255)
                                                                                                                                cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                      >  4:149609/1‑149 (MQ=255)
                                                                                                                                 aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                     >  1:391615/1‑149 (MQ=255)
                                                                                                                                 aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                     >  1:120385/1‑149 (MQ=255)
                                                                                                                                 aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTCTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                     >  3:680334/1‑149 (MQ=255)
                                                                                                                                    gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                  >  1:528305/1‑149 (MQ=255)
                                                                                                                                       taaCAGCCAACGGTTAGGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCCCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgc               >  2:199805/2‑149 (MQ=255)
                                                                                                                                       gaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgc               <  2:894933/149‑1 (MQ=255)
                                                                                                                                       gaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgc               >  4:826493/1‑149 (MQ=255)
                                                                                                                                       gaaCAGCCAACCGTTATGGCGGTTGGCTATGGGGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgc               >  4:230046/1‑149 (MQ=255)
                                                                                                                                         aCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa                                                          >  4:194863/1‑104 (MQ=255)
                                                                                                                                         aCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa                                                          <  3:194863/104‑1 (MQ=255)
                                                                                                                                           aGCCAACCGCCATAACGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGa           >  4:39493/1‑149 (MQ=255)
                                                                                                                                             ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTGGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa         <  3:629653/149‑1 (MQ=255)
                                                                                                                                             ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa         <  4:624534/149‑1 (MQ=255)
                                                                                                                                              cAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg        <  1:565494/149‑1 (MQ=255)
                                                                                                                                               aaCCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAgg       >  3:619704/1‑149 (MQ=255)
                                                                                                                                               aaCCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAgg       >  2:777084/1‑149 (MQ=255)
                                                                                                                                                aCCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGga      <  3:852639/149‑1 (MQ=255)
                                                                                                                                                 ccGTTATGGCGGTTGGCTATGGTGGGAAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGga      <  3:224229/149‑1 (MQ=255)
                                                                                                                                                 ccGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                        >  4:73313/1‑98 (MQ=255)
                                                                                                                                                 ccGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                        <  3:73313/98‑1 (MQ=255)
                                                                                                                                                   gTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaa   <  1:683450/149‑1 (MQ=255)
                                                                                                                                                    ttATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaaa  <  3:447529/149‑1 (MQ=255)
                                                                                                                                                    ttATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaaa  <  3:230046/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAA  >  NC_000913/1982352‑1982648

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