Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,678 T→A 5.2% intergenic (+32/‑334) fadL → / → yfdF long‑chain fatty acid outer membrane transporter/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6780TA 5.2% 810.8 / 25.6 232intergenic (+32/‑334)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (113/107);  new base A (6/6);  total (119/113)
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.39e-01

Only 100 of 237 total aligned reads displayed.

GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462530‑2462826
                                                                                                                                                    |                                                                                                                                                     
gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCAt                                                                                                                                                       >  1:932574/1‑149 (MQ=255)
   gTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑Gc                                                                                                                                                   >  4:377380/1‑149 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATCCGAGTTCGAGTCTGAAGGTAAAGCCTGGCGGTTCGGTACTAACTTTAACTACGCGTGCTGATAGCGCGTTCGCCTGGATAAAGTCACCTGCATA‑GTAGGTGa                                                                                                                                             <  4:901610/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                             <  1:1712271/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                             <  3:1006719/149‑1 (MQ=255)
          tttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGAc                                                                                                                                            >  3:383061/1‑149 (MQ=255)
          tttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGAc                                                                                                                                            >  2:797985/1‑149 (MQ=255)
                atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAAc                                                                                                                                      >  1:296301/1‑149 (MQ=255)
                   tGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcc                                                                                                                                   <  3:726751/149‑1 (MQ=255)
                   tGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcc                                                                                                                                   <  1:692859/149‑1 (MQ=255)
                    gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTccc                                                                                                                                  <  3:444542/149‑1 (MQ=255)
                    gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTccc                                                                                                                                  <  2:485201/149‑1 (MQ=255)
                     cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                                 <  4:1351070/149‑1 (MQ=255)
                     cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                                 <  2:932574/149‑1 (MQ=255)
                     cACGGTCAGAGCATGAAAATTAACGAAGTCCCCTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTTGGTACTAACTTTAACTACGCGTTCTGTTAACGCGTGCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGAATTTAACTcccc                                                                                                                                 <  2:450757/149‑1 (MQ=255)
                      aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTccccc                                                                                                                                <  4:598100/149‑1 (MQ=255)
                        ggTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCAc                                                                                                                              <  2:1021197/149‑1 (MQ=255)
                             gagCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGAc                                                                                                                                            >  2:432070/1‑130 (MQ=255)
                             gagCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGAc                                                                                                                                            <  1:432070/130‑1 (MQ=255)
                                 gTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     >  1:879213/1‑149 (MQ=255)
                                 gTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     >  4:1526191/1‑149 (MQ=255)
                                  tGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                                    <  4:233845/149‑1 (MQ=255)
                                     aaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACt                                                                                                                                           >  1:106960/1‑123 (MQ=255)
                                      aaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCa                                                                                                                          >  1:1145015/1‑139 (MQ=255)
                                      aaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCa                                                                                                                          <  2:1145015/139‑1 (MQ=255)
                                        ttAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTGTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCAACGAAAGTAgtg                                                                                                              >  2:1174056/1‑149 (MQ=255)
                                        ttAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAAATTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAgtg                                                                                                              >  1:3899/1‑149 (MQ=255)
                                        tcagcGAAGGCCCATACCAGTTCGAGGCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTGACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAgtg                                                                                                              >  4:1469589/5‑149 (MQ=255)
                                          aaCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTc                                                                                                            >  4:929362/1‑149 (MQ=255)
                                          aaCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTc                                                                                                            >  2:429135/1‑149 (MQ=255)
                                             gAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTc                                                                                                         >  1:1313730/1‑149 (MQ=255)
                                              aaGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCa                                                                                                        >  1:1415851/1‑149 (MQ=255)
                                                ggCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                      <  3:679152/149‑1 (MQ=255)
                                                   ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTaa                                                                                                                                       >  1:1015049/1‑113 (MQ=255)
                                                    cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                                  >  1:1257020/1‑149 (MQ=255)
                                                    cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                                  <  1:1116653/149‑1 (MQ=255)
                                                      tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAc                                                                                                <  2:1689168/149‑1 (MQ=255)
                                                               gAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                               <  2:1369687/93‑1 (MQ=255)
                                                               gAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                               >  1:1369687/1‑93 (MQ=255)
                                                                  tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagta                                                                                                                 <  1:1053572/120‑1 (MQ=255)
                                                                  tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                                    <  4:382484/149‑1 (MQ=255)
                                                                     gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                                 >  3:1327888/1‑85 (MQ=255)
                                                                     gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCa                                                                                 >  2:1260265/1‑149 (MQ=255)
                                                                     gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCa                                                                                 >  1:971907/1‑149 (MQ=255)
                                                                        ggTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAAc                                                                              <  2:1710290/149‑1 (MQ=255)
                                                                         gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACa                                                                             >  3:1473635/1‑149 (MQ=255)
                                                                            aaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACttt                                                                                                                                         <  4:575441/86‑1 (MQ=38)
                                                                               ccTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGt                                                                       <  2:1313730/149‑1 (MQ=255)
                                                                                 tGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTg                                                                     <  3:884314/149‑1 (MQ=255)
                                                                                  ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGACAGTAGTGTCCTCATTGCTTACCTCCTGTGTTTTGCAAACAGCCTGTTgg                                                                    >  4:24829/1‑149 (MQ=255)
                                                                                           gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                           >  4:169775/1‑149 (MQ=255)
                                                                                               tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                       <  3:1536188/149‑1 (MQ=255)
                                                                                                  cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                                    >  2:1096460/1‑149 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTg                                                                     >  3:1279134/1‑131 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTg                                                                     <  4:1279134/131‑1 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   >  4:764436/1‑149 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   <  4:1245682/149‑1 (MQ=255)
                                                                                                    ttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAAt                                                  >  3:579625/1‑149 (MQ=255)
                                                                                                        cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                           >  4:68198/1‑136 (MQ=255)
                                                                                                        cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                              <  2:362503/149‑1 (MQ=255)
                                                                                                         tccgcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATa                                             >  3:1102599/3‑149 (MQ=255)
                                                                                                            gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                          >  3:1365131/1‑149 (MQ=255)
                                                                                                            gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                          <  3:942732/149‑1 (MQ=255)
                                                                                                               ttCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGt                                       <  3:205327/149‑1 (MQ=255)
                                                                                                                 cTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   <  4:476282/135‑1 (MQ=255)
                                                                                                                 cTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   >  3:476282/1‑135 (MQ=255)
                                                                                                                    aTAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                          >  2:1459628/1‑141 (MQ=255)
                                                                                                                    aTAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCca                                  <  3:1714949/149‑1 (MQ=255)
                                                                                                                      aaCGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                                >  3:1493963/1‑149 (MQ=255)
                                                                                                                         gcgTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCGGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGAc                             >  1:755259/1‑149 (MQ=255)
                                                                                                                         gcgTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGAc                             >  4:804040/1‑149 (MQ=255)
                                                                                                                               gCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCt                       >  1:1721932/1‑149 (MQ=255)
                                                                                                                               gCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCt                       >  2:649223/1‑149 (MQ=255)
                                                                                                                                 cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                     >  3:1185792/1‑149 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  3:688927/149‑1 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  2:1721932/149‑1 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    >  2:735446/1‑149 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  1:1342042/149‑1 (MQ=255)
                                                                                                                                   ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                   >  4:1223825/1‑149 (MQ=255)
                                                                                                                                     aTAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAggg                 >  2:560940/1‑149 (MQ=255)
                                                                                                                                      tAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGa                >  4:560817/1‑149 (MQ=255)
                                                                                                                                       aaaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa               >  3:25006/1‑149 (MQ=255)
                                                                                                                                         aGTCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                         >  4:157783/1‑73 (MQ=38)
                                                                                                                                         aGTCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                         <  3:157783/73‑1 (MQ=37)
                                                                                                                                         aGTCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg             >  2:299402/1‑149 (MQ=255)
                                                                                                                                          gTCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAgg            >  1:1194398/1‑149 (MQ=255)
                                                                                                                                          gTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCt                                                                         >  4:1427171/1‑88 (MQ=255)
                                                                                                                                          gTCACCTGC‑TATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAgg            >  3:379998/1‑149 (MQ=255)
                                                                                                                                           tCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAc                         <  2:1539472/135‑1 (MQ=255)
                                                                                                                                           tCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAc                         >  1:1539472/1‑135 (MQ=255)
                                                                                                                                           tCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                   >  3:82802/1‑141 (MQ=255)
                                                                                                                                           tCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGc           <  2:68200/149‑1 (MQ=255)
                                                                                                                                             aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat         >  4:821229/1‑149 (MQ=255)
                                                                                                                                              ccTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCg                                         <  1:408761/116‑1 (MQ=255)
                                                                                                                                              ccTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCg                                         >  2:408761/1‑116 (MQ=255)
                                                                                                                                               cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGa                              <  4:1027164/126‑1 (MQ=255)
                                                                                                                                                tGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAtt      >  4:59185/1‑149 (MQ=255)
                                                                                                                                                 gCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAttt     >  1:518707/1‑149 (MQ=255)
                                                                                                                                                 gCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGCTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAttt     >  2:876337/1‑149 (MQ=255)
                                                                                                                                                    tA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCa  <  3:764436/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462530‑2462826

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