Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,501 T→C 6.0% 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.0% 196.7 / 3.3 67intergenic (‑114/+53)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

AGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAG  >  NC_000913/1982356‑1982645
                                                                                                                                                 |                                                                                                                                                
aGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATgg                                                                                                                                              >  1:1029439/1‑150 (MQ=255)
aGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATgg                                                                                                                                              >  1:1023491/1‑150 (MQ=255)
        ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGc                                                                                                                                      <  2:16833/150‑1 (MQ=255)
        ggTTTCGGTTAACGGTTCTGTCACCCGGTTCTACTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGAGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                      <  2:1351517/150‑1 (MQ=255)
                ttAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggcgg                                                                                                                              <  2:352061/150‑1 (MQ=255)
                 tAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTggg                                                                                                                             <  1:438531/150‑1 (MQ=255)
                          tGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGc                                                                                                                    >  1:388452/1‑150 (MQ=255)
                           gTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCt                                                                                                                   >  2:790775/1‑150 (MQ=255)
                             cACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTaa                                                                                                                 >  2:450086/1‑150 (MQ=255)
                              aCCCGGTTCTCCTCCTTCTCTTTTCTTAATGTGTTGCGCAGAAAGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAAACGTTATGGCGGTTGGCTAGTGTGGGAAAAACAGCTaaa                                                                                                                >  1:84095/1‑150 (MQ=255)
                                     tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCTCGTTACAGACATGTAAGTATAGACAGTGGGACGGGTGTAGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGc                                                                                                         <  2:81086/150‑1 (MQ=255)
                                     tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGc                                                                                                         <  1:15865/150‑1 (MQ=255)
                                      ctcctcCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTACGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCa                                                                                                        <  2:83149/150‑1 (MQ=255)
                                       tcctccTTCTCTTTTCTTATTGTGTTCCTCAGAACGAGCTTTACAGACATGTAAGTATAGACAGTGTGACGCGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAg                                                                                                       <  2:934713/150‑1 (MQ=255)
                                        cctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCACCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGa                                                                                                      >  2:693224/1‑150 (MQ=255)
                                        cctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTCACAGACATGTAAGTATAGACAGTGTGACGGGGGAAGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGa                                                                                                      >  1:853092/1‑150 (MQ=255)
                                              tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAg                                                                                                <  2:765819/150‑1 (MQ=255)
                                              tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAg                                                                                                >  2:877647/1‑150 (MQ=255)
                                               ctctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGc                                                                                               <  1:295436/150‑1 (MQ=255)
                                                tctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCa                                                                                              >  1:781863/1‑150 (MQ=255)
                                                          ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTg                                                                                    <  2:1029439/150‑1 (MQ=255)
                                                          ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTg                                                                                    <  2:1023491/150‑1 (MQ=255)
                                                            gtgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCg                                                                                  >  2:551088/1‑150 (MQ=255)
                                                              tttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCAAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt                                                                                <  2:739566/149‑1 (MQ=255)
                                                               ttGCGCAGAACGCGCTTTACAGACATGTTAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGttt                                                                               <  1:999952/150‑1 (MQ=255)
                                                                tGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc                                                                              <  1:581435/150‑1 (MQ=255)
                                                                tGCGCAGAACGCGCTATACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc                                                                              <  2:952891/150‑1 (MQ=255)
                                                                   gcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCt                                                                           <  2:549725/150‑1 (MQ=255)
                                                                    cAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtt                                                                          >  2:1127802/1‑150 (MQ=255)
                                                                        aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                      <  1:551088/150‑1 (MQ=255)
                                                                         cgcgCTTTACAGACATGTAAGTATAGAGAGTGTGACGGGGGTCGCCATTTGCAACGAAGACCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGc                                                                     <  1:1214676/150‑1 (MQ=255)
                                                                           cgcTTTACAGACATGTTAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                   <  1:20828/150‑1 (MQ=255)
                                                                             cTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGt                                                                 >  1:856876/1‑150 (MQ=255)
                                                                                tataGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc                                                              >  2:923428/4‑150 (MQ=255)
                                                                                  cAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCCTTATGCCGGTTGGCAATGGTGGGAAAAAACGCTAAATAGTTGCAGAAAAAAGCATCAGCCAGTGCGTTACGCTATTTGTTTGTAACgc                                                            >  2:1036878/1‑149 (MQ=255)
                                                                                      cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                        >  2:837977/1‑150 (MQ=255)
                                                                                      cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa                                                        >  1:1138796/1‑150 (MQ=255)
                                                                                         gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAGCGTTCGAAGa                                                     >  1:84135/1‑150 (MQ=255)
                                                                                           aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATc                                                   <  1:332646/150‑1 (MQ=255)
                                                                                            aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCa                                                  >  1:14305/1‑150 (MQ=255)
                                                                                             gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                 >  1:518142/1‑150 (MQ=255)
                                                                                             gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAACCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                 >  1:604623/1‑150 (MQ=255)
                                                                                              tataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAcc                                                >  1:1136485/1‑150 (MQ=255)
                                                                                                taGACAGTGTGACGTGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAATAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATAACCGc                                              >  2:1052745/1‑150 (MQ=255)
                                                                                                  gACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGTCGAAGATCACCGCTg                                            >  2:44181/1‑150 (MQ=255)
                                                                                                          gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAgg                                    <  2:14305/150‑1 (MQ=255)
                                                                                                          gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAgg                                    >  1:599698/1‑150 (MQ=255)
                                                                                                             gggggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAtt                                 <  1:347445/150‑1 (MQ=255)
                                                                                                              ggggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAATCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                <  1:790775/150‑1 (MQ=255)
                                                                                                              ggggGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                <  1:1043300/150‑1 (MQ=255)
                                                                                                                ggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACAGCTGCCAGCTGGATTAAg                              >  2:112001/1‑150 (MQ=255)
                                                                                                                 gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc                             <  1:628320/150‑1 (MQ=255)
                                                                                                                 gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGGCGAAGATCACCGCTGCCAGCAGGATTAAGc                             <  1:379110/150‑1 (MQ=255)
                                                                                                                  tCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGATAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGcc                            <  1:149683/150‑1 (MQ=255)
                                                                                                                   cGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcg                           <  1:877647/150‑1 (MQ=255)
                                                                                                                           gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                   >  1:1350930/1‑150 (MQ=255)
                                                                                                                               cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTTCCAGCAGGATTAAGCCGCGAACCACGTAc               >  1:80231/1‑150 (MQ=255)
                                                                                                                               cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc               >  1:727709/1‑150 (MQ=255)
                                                                                                                                gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt              >  1:1079382/1‑150 (MQ=255)
                                                                                                                                  accataGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgc            >  1:361506/6‑150 (MQ=255)
                                                                                                                                       aGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGAAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTTCCAGCAGTATTAAGCCGCGAACCACGTTCTGCGCGaa       >  2:527810/1‑150 (MQ=255)
                                                                                                                                        gCCAACCGTGATGGCGGTTGGCTATGGGGGGAAAAAACGCGAAATTGTGGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGCAGATCCCCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa      >  2:959421/1‑150 (MQ=255)
                                                                                                                                         ccAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg     <  2:361506/150‑1 (MQ=255)
                                                                                                                                          cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGATGATCACCGCTGCCAGCAGGTTTAAGCCGCGAACCCCGTAGTGCGCGAAAgg    >  2:99676/1‑150 (MQ=255)
                                                                                                                                          cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAgg    >  2:328945/1‑150 (MQ=255)
                                                                                                                                          cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTCCGCGAAAgg    >  2:826097/1‑150 (MQ=255)
                                                                                                                                            aCCGTTTTGGCGGTTGGCTATGGTGGGAAAAAAAGCTAATTTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag  <  2:84095/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAG  >  NC_000913/1982356‑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: