Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,679 A→T 5.1% intergenic (+33/‑333) 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,6790AT 5.1% 341.0 / 6.1 97intergenic (+33/‑333)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base T (2/3);  total (46/51)
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.31e-01

ACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462531‑2462826
                                                                                                                                                    |                                                                                                                                                    
aCGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATa                                                                                                                                                      <  3:459947/149‑1 (MQ=255)
   ttGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                                  <  1:197314/149‑1 (MQ=255)
    tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                                 <  1:395560/149‑1 (MQ=255)
    tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                                 <  3:614365/149‑1 (MQ=255)
      gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                               <  1:110694/149‑1 (MQ=255)
       tgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTg                                                                                                                                              >  2:19378/1‑149 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                             <  1:878354/149‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                             <  1:669461/149‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                             <  1:19378/149‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                             <  2:687574/149‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                             <  4:181394/149‑1 (MQ=255)
         tttCTTATATGCACGGTCAGAGCGTGTGAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTAGGTACTAACTTTAACTACGCGTTCTGATAACTCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGAc                                                                                                                                            >  1:142654/1‑149 (MQ=255)
          ttCCTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACt                                                                                                                                           <  4:484341/149‑1 (MQ=255)
           tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACtt                                                                                                                                          <  3:476726/149‑1 (MQ=255)
            cTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACttt                                                                                                                                         >  3:447026/1‑149 (MQ=255)
             ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                                  <  4:367886/139‑1 (MQ=255)
             ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                                  >  3:367886/1‑139 (MQ=255)
             ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTa                                                                                                                                        >  4:249483/1‑149 (MQ=255)
              tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTATTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTaa                                                                                                                                       <  2:178538/149‑1 (MQ=255)
               atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACTTTatc                                                                                                                                      <  4:498061/149‑3 (MQ=255)
               atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAAc                                                                                                                                      <  1:458819/149‑1 (MQ=255)
                    cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                                 <  2:424900/149‑1 (MQ=255)
                    cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                                 <  1:414307/149‑1 (MQ=255)
                                gTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                     >  3:38155/1‑149 (MQ=255)
                                 tGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                                    <  1:581082/149‑1 (MQ=255)
                                  gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAag                                                                                                                   >  2:378314/1‑149 (MQ=255)
                                  gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAag                                                                                                                   >  3:797841/1‑149 (MQ=255)
                                      aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                               <  1:658266/149‑1 (MQ=255)
                                      aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                               <  3:653361/149‑1 (MQ=255)
                                      aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                               >  1:8654/1‑149 (MQ=255)
                                          aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTcc                                                                                                           >  4:136977/1‑149 (MQ=255)
                                          aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTcc                                                                                                           >  4:441371/1‑149 (MQ=255)
                                                  ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCt                                                                                                   >  3:297428/1‑149 (MQ=255)
                                                   cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                                  >  3:168503/1‑149 (MQ=255)
                                                   cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                                  <  4:11750/149‑1 (MQ=255)
                                                     tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                      >  1:125375/1‑143 (MQ=255)
                                                     tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                      <  2:125375/143‑1 (MQ=255)
                                                     tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAc                                                                                                <  3:442737/149‑1 (MQ=255)
                                                     tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGCCACCTGCATA‑GCAGGTGACTTTAACTCCACCACTTCACCGAAAGTAGTGTCCTCATTGCTTAc                                                                                                <  2:613998/149‑1 (MQ=255)
                                                      aCCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCTTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcc                                                                                               <  2:227943/149‑1 (MQ=255)
                                                          gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATTACGTGTGCGCCGGGATAAAGTCACCAGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCTAAAGTAGTGTCCTCATTGCTTAcctcct                                                                                           <  2:638658/149‑1 (MQ=255)
                                                            tCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                         <  3:615048/149‑1 (MQ=255)
                                                            tCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                         <  3:819310/149‑1 (MQ=255)
                                                              gAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGt                                                                                       >  1:16099/1‑149 (MQ=255)
                                                                      aGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               >  1:504083/1‑149 (MQ=255)
                                                                         tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAg                                                                            >  2:303870/1‑149 (MQ=255)
                                                                          aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                           <  1:340778/149‑1 (MQ=255)
                                                                          aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                           <  1:869403/149‑1 (MQ=255)
                                                                            aGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCt                                                                         <  4:447026/149‑1 (MQ=255)
                                                                               cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                      <  1:378314/149‑1 (MQ=255)
                                                                                 ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                    >  2:367291/1‑149 (MQ=255)
                                                                                 ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                    >  3:593102/1‑149 (MQ=255)
                                                                                       tCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCtt                                                              <  1:252474/149‑1 (MQ=255)
                                                                                          gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                           >  4:524101/1‑149 (MQ=255)
                                                                                          gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                           <  1:110110/149‑1 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               <  3:18168/126‑1 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               >  4:18168/1‑126 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                        >  1:87542/1‑149 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                        >  3:595337/1‑149 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                        >  1:510754/1‑149 (MQ=255)
                                                                                              tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                       >  3:842864/1‑149 (MQ=255)
                                                                                                 cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                                    >  3:267211/1‑149 (MQ=255)
                                                                                                  tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   >  3:397935/1‑149 (MQ=255)
                                                                                                  tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                   >  1:730418/1‑149 (MQ=255)
                                                                                                      aCTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCaa                                               <  3:249483/149‑1 (MQ=255)
                                                                                                       cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                              <  4:297428/149‑1 (MQ=255)
                                                                                                       cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                              <  2:470190/149‑1 (MQ=255)
                                                                                                        tACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCGGCTTGCCCTTTCAAATCAATa                                             >  1:844497/1‑149 (MQ=255)
                                                                                                             gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCgg                                        >  4:61736/1‑149 (MQ=255)
                                                                                                             gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCgg                                        >  1:362429/1‑149 (MQ=255)
                                                                                                                 tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                                 >  1:209377/1‑136 (MQ=255)
                                                                                                                 tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                                 <  2:209377/136‑1 (MQ=255)
                                                                                                                     aaCGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                                >  1:297366/1‑149 (MQ=255)
                                                                                                                            tCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAc                         <  3:461006/149‑1 (MQ=255)
                                                                                                                                cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGGAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               >  2:767310/1‑91 (MQ=255)
                                                                                                                                cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                     >  3:252624/1‑149 (MQ=255)
                                                                                                                                cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                     >  2:549722/1‑149 (MQ=255)
                                                                                                                                cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTGAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                               <  1:767310/91‑1 (MQ=255)
                                                                                                                                 tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    >  1:525345/1‑149 (MQ=255)
                                                                                                                                 tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                    <  2:586108/149‑1 (MQ=255)
                                                                                                                                  ggATAAAGTCACCTGCATA‑GCAGGTGCCTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                   <  2:142654/149‑1 (MQ=255)
                                                                                                                                  ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                   <  1:303870/149‑1 (MQ=255)
                                                                                                                                   gATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                    >  3:596072/1‑131 (MQ=255)
                                                                                                                                   gATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                    <  4:596072/131‑1 (MQ=255)
                                                                                                                                    aTAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAggg                 >  1:175226/1‑149 (MQ=255)
                                                                                                                                     tAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGa                <  4:397935/149‑1 (MQ=255)
                                                                                                                                      aaaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa               >  3:225062/1‑149 (MQ=255)
                                                                                                                                       aaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa              >  4:701763/1‑149 (MQ=255)
                                                                                                                                       aaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa              <  1:175229/149‑1 (MQ=255)
                                                                                                                                        aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg             <  1:232745/149‑1 (MQ=255)
                                                                                                                                        aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg             <  2:643126/149‑1 (MQ=255)
                                                                                                                                        aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg             <  3:45234/149‑1 (MQ=255)
                                                                                                                                          tCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGc           >  4:736711/1‑149 (MQ=255)
                                                                                                                                            aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATcc                                   >  2:269385/1‑123 (MQ=255)
                                                                                                                                            aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATcc                                   <  1:269385/123‑1 (MQ=255)
                                                                                                                                            aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat         >  4:499962/1‑149 (MQ=255)
                                                                                                                                                gCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTACCCTTTCAAATGAATAAGCGGTATCCACAGAGAAACCTGAAGGGAATGGCATAttt     >  2:117598/1‑149 (MQ=255)
                                                                                                                                                   tA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTATCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCa  >  4:37948/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCA  >  NC_000913/2462531‑2462826

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