Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,501 T→C 5.1% intergenic (‑114/+53) 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,5010TC 5.1% 382.8 / 6.6 119intergenic (‑114/+53)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (60/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 = 7.69e-01

Only 100 of 123 total aligned reads displayed.

GATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAAT  >  NC_000913/1982353‑1982649
                                                                                                                                                    |                                                                                                                                                     
gATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                       <  2:814150/149‑1 (MQ=255)
  tAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTAt                                                                                                                                                     <  4:69351/149‑1 (MQ=255)
    gTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATgg                                                                                                                                                   <  4:133406/149‑1 (MQ=255)
     ttCAGGGGTTTCGGTTAACGGTTCTGTTACCCGGTTCTCCTCCTTCTCTTATCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGTACAGCCCACCGTTATGGc                                                                                                                                                  <  2:738697/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)
             tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTAt                                                                                                                                                     >  2:257760/1‑138 (MQ=255)
             tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                          >  4:314465/1‑149 (MQ=255)
             tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                          >  2:855887/1‑149 (MQ=255)
             tttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                          >  4:441007/1‑149 (MQ=255)
                  gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                     >  2:620159/1‑149 (MQ=255)
                         gTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaa                                                                                                                             >  4:852639/1‑149 (MQ=255)
                           tCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aacaaa                                                                                                                           <  1:771876/149‑1 (MQ=255)
                           tCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaaa                                                                                                                           <  4:781137/149‑1 (MQ=255)
                               tCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCt                                                                                                                       >  2:511212/1‑149 (MQ=255)
                               tCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCt                                                                                                                       <  3:441007/149‑1 (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)
                                    cGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAAtt                                                                                                                  <  1:482623/149‑1 (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                                                                                                         >  1:307274/1‑149 (MQ=255)
                                             tcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaa                                                                                                         >  4:226551/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)
                                                ttCTCTTTTCTTATTGTGTTCCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCTCCATTTGCAGCGAAGAACAGCGAACCGTAATGGCGGGTGACTATGGTGGG‑TAAAGACGCTAAGTTGGTGTAGaaaaa                                                                                                      >  2:201551/1‑149 (MQ=255)
                                                 tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGaaaaaa                                                                                                     <  4:463274/149‑1 (MQ=255)
                                                   tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGc                                                                                                   <  1:430760/149‑1 (MQ=255)
                                                   tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGc                                                                                                   >  1:817243/1‑149 (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)
                                                        tCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAg                                                                                              >  3:139522/1‑149 (MQ=255)
                                                          ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGAc                                                                                            >  2:371432/1‑149 (MQ=255)
                                                          ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTgg                                                                                                                                            <  3:452924/102‑1 (MQ=255)
                                                          ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTgg                                                                                                                                            >  4:452924/1‑102 (MQ=255)
                                                           tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑aaaaaa                                                                                                                           >  2:526556/1‑117 (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)
                                                                  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:108487/1‑149 (MQ=255)
                                                                           aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                           <  1:292209/149‑1 (MQ=255)
                                                                           aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                           >  3:315193/1‑149 (MQ=255)
                                                                             gcggTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCGATCGTTATGGATGGTGGCTATGGTGGG‑ATAAAACGCTAAATTGTTGCAGAAAAAAGCATAATACAGATCGTTTCGCTTTTTGc                                                                         <  2:543035/149‑1 (MQ=255)
                                                                             gcgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                         <  2:728403/149‑1 (MQ=255)
                                                                             gcgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                         <  3:551720/149‑1 (MQ=255)
                                                                               gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                       <  3:506245/149‑1 (MQ=255)
                                                                                      aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                                >  2:533447/1‑149 (MQ=255)
                                                                                          aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                            >  4:819648/1‑149 (MQ=255)
                                                                                          aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                            >  1:863252/1‑149 (MQ=255)
                                                                                          aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                            >  1:543944/1‑149 (MQ=255)
                                                                                              aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGAt                                                        <  3:286434/149‑1 (MQ=255)
                                                                                                 tataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                     <  3:226551/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                                 >  1:784077/1‑149 (MQ=255)
                                                                                                                     tCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                                 >  2:775965/1‑149 (MQ=255)
                                                                                                                           tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGaa                           <  4:461414/149‑1 (MQ=255)
                                                                                                                           tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGaa                           >  2:781479/1‑149 (MQ=255)
                                                                                                                               cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                       >  4:149609/1‑149 (MQ=255)
                                                                                                                                aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                      >  1:120385/1‑149 (MQ=255)
                                                                                                                                aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                      >  1:391615/1‑149 (MQ=255)
                                                                                                                                aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTCTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                      >  3:680334/1‑149 (MQ=255)
                                                                                                                                 aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGt                     >  1:653374/1‑149 (MQ=255)
                                                                                                                                   gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                   >  1:528305/1‑149 (MQ=255)
                                                                                                                                    caccataGCCAACCGCCATAACGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt                  >  2:116540/7‑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                                                           <  3:194863/104‑1 (MQ=255)
                                                                                                                                        aCAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa                                                           >  4:194863/1‑104 (MQ=255)
                                                                                                                                          aGCCAACCGCCATAACGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGa            >  4:39493/1‑149 (MQ=255)
                                                                                                                                            ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTGGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa          <  3:629653/149‑1 (MQ=255)
                                                                                                                                            ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa          <  1:555552/149‑1 (MQ=255)
                                                                                                                                            ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa          <  2:287767/149‑1 (MQ=255)
                                                                                                                                            ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa          <  4:624534/149‑1 (MQ=255)
                                                                                                                                            ccAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa          <  3:616454/149‑1 (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)
                                                                                                                                                   ttATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaaa   <  3:447529/149‑1 (MQ=255)
                                                                                                                                                    tATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAat  <  3:644345/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAAT  >  NC_000913/1982353‑1982649

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