Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,501 T→C 6.6% 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 6.6% 164.8 / 2.7 65intergenic (‑114/+53)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (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 = 8.82e-01

GATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGA  >  NC_000913/1982353‑1982646
                                                                                                                                                    |                                                                                                                                                 
gATAGTTAAGGGTTGTCGGTTAACGGTGCTGGCACCCGTTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                   <  1:684113/149‑1 (MQ=255)
     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)
                                                                                                                                                    |                                                                                                                                                 
GATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGA  >  NC_000913/1982353‑1982646

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: