Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,500 T→C 10.6% 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,5000TC10.6% 239.7 / 16.4 94intergenic (‑113/+54)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/49);  new base C (4/6);  total (39/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.95e-01

CGGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAG  >  NC_000913/1982351‑1982645
                                                                                                                                                     |                                                                                                                                                  
cgGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGt                                                                                                                                                    >  2:1545370/1‑150 (MQ=255)
 ggATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                   <  2:1669951/150‑1 (MQ=255)
 ggATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                   <  2:54251/150‑1 (MQ=255)
  gATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTa                                                                                                                                                  >  1:761589/1‑150 (MQ=255)
      gTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGc                                                                                                                                              <  2:382743/150‑1 (MQ=255)
      gTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGATTTACTGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGACCAGCCATCCGGTATGGc                                                                                                                                              >  2:918530/1‑150 (MQ=255)
       ttCTGGGGTTTCGGTTATCGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGGTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCg                                                                                                                                             <  2:155062/150‑1 (MQ=255)
        tCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTGTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCgg                                                                                                                                            <  2:322179/150‑1 (MQ=255)
             ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAATAGCCAACCGCCATAACGGTTGGc                                                                                                                                       <  1:1339936/150‑1 (MQ=255)
             ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGc                                                                                                                                       >  2:1668852/1‑150 (MQ=255)
             ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGc                                                                                                                                       <  1:1714860/150‑1 (MQ=255)
             ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                       <  1:29741/150‑1 (MQ=255)
             ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                       <  2:437818/150‑1 (MQ=255)
              gtatCGGTTAACGGTGCTGTCACCCGGTTCTCCTCCATCTCTTTTCTTATTGTGATGCGCAGAACGCGCATTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCt                                                                                                                                      <  2:751583/147‑1 (MQ=255)
                ttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGCtgt                                                                                                                                    <  2:482124/150‑3 (MQ=255)
                 tCGGTTCACGGCTCTGTCACCCGGTGCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATTGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGCtgtt                                                                                                                                   <  2:1556408/150‑4 (MQ=255)
                 tCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtg                                                                                                                                   <  1:588897/150‑1 (MQ=255)
                  cGGTTTACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAATAGCCAACCGCCATAACGGTTGGCtgttc                                                                                                                                  <  2:1235507/150‑5 (MQ=255)
                  cGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtgg                                                                                                                                  >  1:411011/1‑150 (MQ=255)
                   tgTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                 >  2:1324285/2‑150 (MQ=255)
                   ggTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGTGGTTGGCTAtggt                                                                                                                                 >  2:159418/1‑150 (MQ=255)
                    gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggtg                                                                                                                                >  2:189217/1‑150 (MQ=255)
                     ttAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggtgg                                                                                                                               <  2:1445001/150‑1 (MQ=255)
                     ttAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggtgg                                                                                                                               <  1:726224/150‑1 (MQ=255)
                        aCGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGCCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaa                                                                                                                            <  1:1005945/150‑1 (MQ=255)
                        aCGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaa                                                                                                                            <  2:868843/150‑1 (MQ=255)
                        aCGGTTCAGGCACCCGGTGCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCGTTACAGACATGTAAGTATCGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGACCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaa                                                                                                                            <  1:169122/150‑1 (MQ=255)
                           gTTGTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTCTGTTGCGCAGAACGCGCTTTACAGACATGTTAGTATAGACAGTGAGCCGGGGGTCGCCATTTGCAACGTAGTACAGCCAACCGTTATGTCGGTTGGCTATTGTGGGaaaat                                                                                                                         <  2:689507/150‑2 (MQ=255)
                                      cGGGTCCCCTCCAACTCATTTCTTATTGTGGTGCCCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAAttg                                                                                                              <  1:143092/150‑1 (MQ=255)
                                             cctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑a                                                                                                      >  1:817862/1‑150 (MQ=255)
                                              ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑aa                                                                                                     >  2:1173616/1‑150 (MQ=255)
                                                   tctcTTTTCTTATTGTGTTGCGCAGAACTCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGAAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAg                                                                                                <  1:522985/150‑1 (MQ=255)
                                                   tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAg                                                                                                <  2:1431544/150‑1 (MQ=255)
                                                   tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAg                                                                                                <  2:1212484/150‑1 (MQ=255)
                                                     tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCa                                                                                              <  2:566862/150‑1 (MQ=255)
                                                      ctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCAt                                                                                             >  2:618284/1‑150 (MQ=255)
                                                      ctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCAt                                                                                             <  2:761589/150‑1 (MQ=255)
                                                        tttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCa                                                                                           >  2:666121/1‑150 (MQ=255)
                                                         ttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAg                                                                                          >  1:488622/1‑150 (MQ=255)
                                                          tCTTATTGTGTTGCGCAGAACGCGCTTTACTGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGa                                                                                         <  1:945969/150‑1 (MQ=255)
                                                            ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACa                                                                                       >  1:951172/1‑150 (MQ=255)
                                                               ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTg                                                                                    <  1:451603/150‑1 (MQ=255)
                                                               ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTg                                                                                    <  2:484564/150‑1 (MQ=255)
                                                               ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTg                                                                                    <  2:411011/150‑1 (MQ=255)
                                                                     tGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTc                                                                              <  2:975020/150‑1 (MQ=255)
                                                                      gcgcTGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCg                                                                             <  1:28807/150‑1 (MQ=255)
                                                                      gcgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCg                                                                             <  2:899816/150‑1 (MQ=255)
                                                                      gcgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCg                                                                             >  1:1527354/1‑150 (MQ=255)
                                                                         cAGAACGCGCTTTACAGACATGTACGTATCGACAGTGTGCCGGGGGTCACCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAACCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCtt                                                                          <  2:951121/150‑1 (MQ=255)
                                                                             aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCATCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                      <  2:149503/150‑1 (MQ=255)
                                                                             aCGCGATTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCAATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                      <  2:1255991/150‑1 (MQ=255)
                                                                                cgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                   <  2:601077/150‑1 (MQ=255)
                                                                                 gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGTGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTg                                                                  <  1:159418/150‑1 (MQ=255)
                                                                                 gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTg                                                                  >  2:1712815/1‑150 (MQ=255)
                                                                                      aCAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGCGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTACGCTTTTTGCTTGTAACg                                                             >  2:366578/1‑150 (MQ=255)
                                                                                       cAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                            <  1:1094597/150‑1 (MQ=255)
                                                                                       cAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCAGTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCTGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                            >  1:216196/1‑150 (MQ=255)
                                                                                        aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGt                                                           <  2:1625884/150‑1 (MQ=255)
                                                                                          aCATTTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCg                                                         <  2:1519294/150‑1 (MQ=255)
                                                                                             tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAg                                                      >  1:1552457/1‑150 (MQ=255)
                                                                                             tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAg                                                      <  1:1596128/150‑1 (MQ=255)
                                                                                              ttAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                     <  1:1690339/149‑1 (MQ=255)
                                                                                              gttaGTATCGACAGTGTGCCGGGGGTCGGCATTTGCCACGTAGAATAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCAGCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                     <  1:1681323/147‑1 (MQ=255)
                                                                                              gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                     >  1:1269524/1‑150 (MQ=255)
                                                                                                 aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCa                                                  >  1:491463/1‑150 (MQ=255)
                                                                                                 aGTATAGACAGTGGGCCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACAGTTATGTCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCa                                                  <  1:1308444/150‑1 (MQ=255)
                                                                                                  gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCTACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                 >  2:110979/1‑150 (MQ=255)
                                                                                                    taaGACAGTGTGACGTGGGTCGCCATTTGCAACGAAGAACAGCCACCCGTAATGGCGGTTGGCTATGGTGTGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCg                                               <  2:255931/148‑1 (MQ=255)
                                                                                                    ataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCg                                               <  1:1019435/150‑1 (MQ=255)
                                                                                                    ataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCg                                               <  1:1025696/150‑1 (MQ=255)
                                                                                                     taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGc                                              >  1:1324575/1‑150 (MQ=255)
                                                                                                       gACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTg                                            >  1:1630086/1‑150 (MQ=255)
                                                                                                       gACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTg                                            >  1:1286798/1‑150 (MQ=255)
                                                                                                           gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcag                                        >  2:388708/1‑150 (MQ=255)
                                                                                                               gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAgg                                    <  1:1545370/150‑1 (MQ=255)
                                                                                                                     ggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATAACCGCTGCCAGCAGGATTAAg                              >  2:971540/1‑150 (MQ=255)
                                                                                                                      gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGTTTCGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                             >  1:1559394/1‑150 (MQ=255)
                                                                                                                      gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                             <  2:1652061/150‑1 (MQ=255)
                                                                                                                      gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGAGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                             <  1:258007/150‑1 (MQ=255)
                                                                                                                         gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc                          >  1:25271/1‑150 (MQ=255)
                                                                                                                            aTTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑‑AAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAAc                      <  1:1668852/150‑1 (MQ=255)
                                                                                                                              ttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAAcc                     <  1:1637797/150‑1 (MQ=255)
                                                                                                                               tGCAACGAAGAACAGCCAACAGTTATGGCGGTTGGCTATGTTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACTGTCGAAGATCACCTCTGCCAGCAGGATAAAGCCGCGAAccc                    >  2:1125197/1‑149 (MQ=255)
                                                                                                                                gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGaacaac                   >  1:515436/1‑150 (MQ=255)
                                                                                                                                gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGAGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                   >  2:272266/1‑150 (MQ=255)
                                                                                                                                 cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                  <  2:1527354/150‑1 (MQ=255)
                                                                                                                                 cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                  >  1:1247872/1‑150 (MQ=255)
                                                                                                                                  aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                  <  2:1662604/150‑1 (MQ=255)
                                                                                                                                   aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCAACAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCTCGAACCACGTa                >  1:1255806/1‑150 (MQ=255)
                                                                                                                                     gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCCGGATGAAGCCGCGAACCACGTACt              >  2:1535751/1‑150 (MQ=255)
                                                                                                                                         aaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgcgc          <  1:1497270/150‑1 (MQ=255)
                                                                                                                                          ataGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGAAGGATTAAGCCGCGAACCACGTACTGcgcg         >  1:899163/3‑150 (MQ=255)
                                                                                                                                          ataGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCAAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGcgcg         >  2:94506/3‑150 (MQ=255)
                                                                                                                                           cAGCCAACCGTTATGGCTGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGa        <  2:1667993/150‑1 (MQ=255)
                                                                                                                                            aGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa       >  2:1377421/1‑150 (MQ=255)
                                                                                                                                             gCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGGTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAGTAAGCAGCGAACCACGTACTGAGCGaaa      >  2:134254/1‑150 (MQ=255)
                                                                                                                                                 aCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag  >  1:1400693/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAG  >  NC_000913/1982351‑1982645

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