Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,505 G→A 6.1% intergenic (‑118/+49) 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,5050GA 6.1% 161.0 / 2.3 65intergenic (‑118/+49)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/2);  total (35/30)
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.48e-01

TTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTA  >  NC_000913/1982358‑1982653
                                                                                                                                                   |                                                                                                                                                    
ttCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGCACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGc                                                                                                                                                     >  2:484172/1‑149 (MQ=255)
    ggggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                                 >  2:395007/1‑149 (MQ=255)
    ggggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                                 <  1:289744/149‑1 (MQ=255)
     gtgTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGAGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTAAGGGCGGTTg                                                                                                                                                >  2:972787/3‑149 (MQ=255)
     gggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATAGTGTTGCGAAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGCCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTg                                                                                                                                                >  1:8674/1‑149 (MQ=255)
       ggttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACGGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                              <  1:415446/147‑1 (MQ=255)
       gTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                              <  2:287243/149‑1 (MQ=255)
       gTTTCGGTTAACGGTTCTGAAACCCGGTTCTCCTACTTATCTTTTCTTATTGTGTTGCGCACAAAGCGCTGTACAGACATTTAAGTATAGACAGTGTGACGGGGGTCGCCATGTGCAACGCAGAACAGCCAAACGATCTGGCGGTGGGc                                                                                                                                              >  2:926743/1‑149 (MQ=255)
           ctgTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATTGACAGTGTGACGGGTGTCGCCATTTGCAACGAAGATCAGCCAACCGTTATGGCGGTTGGCTAtg                                                                                                                                          <  1:747633/147‑1 (MQ=255)
           cGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTCGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtg                                                                                                                                          <  1:289694/149‑1 (MQ=255)
            ggTTCACGGTGCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtgg                                                                                                                                         <  1:43501/149‑1 (MQ=255)
             gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTTTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCACAGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                        >  1:844045/1‑149 (MQ=255)
             gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                        >  2:246771/1‑149 (MQ=255)
             gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTATCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                        >  2:91526/1‑149 (MQ=255)
             gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTGTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCGTTTGCAAGGAAGAACAGCCACCCGTTATGGCGGTAGGCTAgggt                                                                                                                                        >  2:580288/1‑149 (MQ=255)
                    gTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaa                                                                                                                                 >  2:433411/1‑149 (MQ=255)
                     ttCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGTAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaaa                                                                                                                                >  2:683545/1‑149 (MQ=255)
                             cccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTaaa                                                                                                                        >  1:94061/1‑149 (MQ=255)
                                   tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGCCATGTAAGTATAGACAGTGGGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAAttgttg                                                                                                                  <  2:577097/149‑1 (MQ=255)
                                   tctcCTCCTTCTCTTTTCTTATTGTCTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAAttgttg                                                                                                                  <  2:220483/149‑1 (MQ=255)
                                     tcctccTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGAAACGAAGAACAGCCATCCGTTATGGCGCTTGGCTATGGTGGGCAAAAACGCTAAATTGTTCCa                                                                                                                >  2:807959/1‑149 (MQ=255)
                                           ttCTGTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCAAACCGTTATGGCGTTTGGCTATGGTGGGAAAACACGCTAAATTGTAGCAGaaaaa                                                                                                          >  1:574512/1‑149 (MQ=255)
                                                    cTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGa                                                                                                 <  2:770576/149‑1 (MQ=255)
                                                       aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCTACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAg                                                                                              >  1:284483/1‑149 (MQ=255)
                                                       aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAg                                                                                              >  1:288433/1‑149 (MQ=255)
                                                             ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt                                                                                        <  1:520891/149‑1 (MQ=255)
                                                             ttGCGCAGAACGCGCGTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt                                                                                        >  1:448442/1‑149 (MQ=255)
                                                               gcgcAGAAGGCGCTTTACAGAAATGTAAGTATCGACAGTGTGACGGGGTTCGCGATTTGCAACGAAGAACAGCCAACCGTTTTGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc                                                                                      <  1:588360/149‑1 (MQ=255)
                                                                cgcAGCTCGCGCGATACAGCAAGGCAAGTATAGACAGTTGGAGCGGGTGCGCCTTTTGCAACGAAGCACAGCCAACCGTTATGGCGGTTTGCTATGGGGGGAAAACTCGCTAAATTGTTGCAGAAAAAAGCAGCCGACAGTGCGtttgg                                                                                     <  2:673394/149‑3 (MQ=255)
                                                                cgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCg                                                                                     <  2:201512/149‑1 (MQ=255)
                                                                   acaaCGCGCTTTACAGACATGTAAGTCTAGACAGGGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAATAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtt                                                                                  <  2:597816/147‑1 (MQ=255)
                                                                    gAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGTAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttt                                                                                 <  1:683545/149‑1 (MQ=255)
                                                                     aaCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtttt                                                                                >  1:678023/1‑149 (MQ=255)
                                                                     aaCGCGCTTTACAGACATGAAAGTATAGACAGTCAGACCGGCCTCGCCATATGCAACGAAGAACAGACAACAGTTATGGCGGTTGGCTATGGTGTGAATAACCGCTAAATTGTTGCATAAAAAAGCATCAGACAGTACCATTCGCtttt                                                                                >  2:465706/1‑149 (MQ=255)
                                                                      aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                               <  1:555546/149‑1 (MQ=255)
                                                                          gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                           <  1:895083/149‑1 (MQ=255)
                                                                            tttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGt                                                                         <  1:246771/149‑1 (MQ=255)
                                                                                 aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                                    >  1:243737/1‑149 (MQ=255)
                                                                                    cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGGTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCTGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCg                                                                 <  1:77531/149‑1 (MQ=255)
                                                                                     aTGTAAGAATAGACAGTGTGACGGGGGTCCACATATGCAACGAAGAACAGCCAACCGATATGGAGGTTGGCTATGGAGGGAAAAAACGCTAAATTGATGCAGAAAAAAGCATCAGACAGAGCGTTTCCCTATTAGCCTGTAACGGACGa                                                                >  1:234481/1‑149 (MQ=255)
                                                                                      tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa                                                               <  2:390171/149‑1 (MQ=255)
                                                                                        tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                             >  1:169433/1‑149 (MQ=255)
                                                                                         aattataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGAt                                                            <  2:844045/146‑1 (MQ=255)
                                                                                         aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACCAAGAACAGCCAACCGTAATGGCGGGTGGCTATGGTGGGAAAAAACGCTAAATTGTTTCAGAAAAAAGCATCAGACAGTGCGTTGCGCTTTTTGCTTGTAACGGACGAAGAt                                                            <  2:97321/149‑1 (MQ=255)
                                                                                                          cGGGCGTCGCCATTTGCAACGAAGAACAGCCCAACGTTATGGCGGGTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAACGCCTCAGACAGCGCGGTTCGCTTTTTGCTTGTAACGGGCGAAGAGCACCGCTGCCAGCAGGa                                           <  2:151314/149‑1 (MQ=255)
                                                                                                             gggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                        >  2:432597/1‑149 (MQ=255)
                                                                                                               gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAg                                      <  1:484172/149‑1 (MQ=255)
                                                                                                                  gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcg                                   <  1:844545/149‑1 (MQ=255)
                                                                                                                   ccATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc                                  >  1:882350/1‑149 (MQ=255)
                                                                                                                         gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAATCCTCGAACCa                            >  1:1024998/1‑149 (MQ=255)
                                                                                                                          cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                           >  2:423347/1‑149 (MQ=255)
                                                                                                                          cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                           <  1:91526/149‑1 (MQ=255)
                                                                                                                          cAACGAAGAACAGCCAACCGTAATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGAAGGATTAAGCCCCGaaaaac                           >  2:425507/1‑145 (MQ=255)
                                                                                                                           aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCCCg                          >  2:472469/1‑149 (MQ=255)
                                                                                                                             cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGCCAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTa                        >  2:483004/1‑149 (MQ=255)
                                                                                                                             cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTa                        >  2:1024174/1‑149 (MQ=255)
                                                                                                                               aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAATACGCGAAATTGTTGCAGAAAAAAACATCAGACAGTGCGTTTCGGTGTTTGCTTGTAACGATCGAACTTCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt                      >  1:984769/1‑149 (MQ=255)
                                                                                                                               aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt                      <  2:956960/149‑1 (MQ=255)
                                                                                                                                accataGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTg                     >  2:126825/6‑149 (MQ=255)
                                                                                                                                 ccataGCCAACCGCCATAACGGTGGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGGCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCCCGTACTgc                    >  1:429259/5‑149 (MQ=255)
                                                                                                                                  ctgagtcAAGCGTTATTGCGGTTGGATTTGGTGGGAAAAAACGCTATATTGTTGCAGAATAAAGCATCAGTCAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAGGTACTgcg                   <  2:1024998/143‑1 (MQ=255)
                                                                                                                                     aGCCACCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGCACCCGTTACTGCGCGa                >  1:159364/1‑149 (MQ=255)
                                                                                                                                       ccAACAGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa              <  1:545175/149‑1 (MQ=255)
                                                                                                                                        cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg             <  2:94061/149‑1 (MQ=255)
                                                                                                                                            cGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTCTAACGGTACAACAACACCTCTCCCAGCAGGATTAAGCCGAGAACCACGTACTTCGCGAAAGgaga         >  1:154692/1‑149 (MQ=255)
                                                                                                                                                   gCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGAAGGATTAAGCCGCTAACCACGTACTGCGCGAAAGGAGAAATATTa  >  2:515054/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTA  >  NC_000913/1982358‑1982653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: